Hi , I have tried solving the following system of ODE's (eq1 attached) using Matlab. first i reduced it into a system of four 1st order ODE's (eq2 attached). thani tried to solve it using ode45() in the following manner: function xprime = Mirage(t,x); k=2; xprime=[x(1); k*(1-x(1).^2)./(1+k*x(2)); x(3) ; -(k*x(1)*x(3))./(1+k*x(2))]; end x0=[1 1 1 1]; tspan=[0,20]; [t,x]=ode45(@Mirage,tspan,x0) the program seem to be in infinite loop(stays Busy). where did i wrong? i will appreciate your help.