Fitzhugh Nagumo code question

Question

0 个投票

I wish to solve the equations v'=v^3/3-w+I, w'=g(v+a-bw) with a=0.8, b=0.7, g=0.08, I=0.5 using ode23/45

I googled this but for one unfamiliar with the code, it is hard to fathom what they are solving

I also would like the code to show the plots of both variables changing with time and the phase plots of both variables.

Thanks

0 个评论
显示 -2更早的评论隐藏 -2更早的评论

Answer 1

Azzi Abdelmalek 2014-7-20

编辑：Azzi Abdelmalek 2014-7-20

0 个投票

Write this function and save it as myode.m

function dy=myode(t,y)
dy=zeros(2,1)
a=0.8 
b=0.7
g=0.08 
I=0.5
dy(1)=y(1)^3/3-y(2)+I
dy(2)=g*(y(1)+a-b*y(2))

Call the above function

tspan=[0 10];
v0=0; 
w0=0;
y0=[v0;w0]; % Initials consitions
[t,y]=ode45(@myode,tspan,y0)

Andrew 2014-7-20

Hi, thanks. It says the variable y is undefined

Azzi Abdelmalek 2014-7-20

Did you do the two steps