function f=ody(t,y)
f=zeros(2,1);
f(1)=y(2);
f(2)=100.*(1-y(1).^2).*y(2)-y(1);
%[T,Y]=ode45(@ody,[0 300],[2 0]);
%plot(T, Y(:,1));
%plot(Y(:,1),Y(:,2));