strfunc ='100*(1-y^2)*yprime - y';
[T,Y1,Y2]=ode45_2(strfunc, [0, 30], [2; 0]);

