Can someone help me with my code ? I can not see error.

function dy = vdp1000(t,y)

dy = zeros(2,1);

dy(1) = y(2);

dy(2) = 1000*(1 - y(1)^2)*y(2) - y(1);

[T,Y] = ode15s(@vdp1000,[0 3000],[2 0]);

plot(T,Y(:,1),'-o')

[Harry]

The code runs fine for me. (I called the function with input 'y' as 1:2). The input variable 't' is never used, so you can delete that.

Answer 1

[Peng Liu]

There is nothing wrong with your code. You just need to save the function vdp1000 as an independent m file with name vdp1000.m and add it to the MATLAB current folder. Then you should be able to run: [T,Y] = ode15s(@vdp1000,[0 3000],[2 0]); plot(T,Y(:,1),'-o')