rem88

HI, i am using ODE45 to solve this:

My M file is:

function dx=p(t,x)

dx=zeros(size(x));

u = (0.177*x(2)/(0.005 + x(2)));

f =((u*3443*20*exp(u*t))/(28*0.45));

v = 3443 + f;

dx(1)=((0.177*x(1)*x(2))/(0.005+x(2))) - ((f*x(1))/v);

dx(2)=(28*f/v) - (0.177*x(1)*x(2)/(0.45*((0.005+x(2)))));


and the code in the work space is:

[t,x]=ode45(@p, [0,4],[28;60]);
plot(t,x)


I get a plot of x1 and x2 against t.

I was just wondering if there is some way of plottin u or f against time?

Apprecaite any help.

Thanks,

rem