04-54.txt

MATLABR2006a基础教程

function dydt = vdp1(t,y)
dydt = [y(2); (1-y(1)^2)*y(2)-y(1)];


>> [t,y] = ode45(@vdp1,[0 20],[2; 0]);

>> plot(t,y(:,1),'-',t,y(:,2),'--')
>> title('Solution of van der Pol Equation, \mu = 1');
>> xlabel('time t');
>> ylabel('solution y');
>> legend('y_1','y_2')