Dear friends,
i am facing a problem. here is an equation,
dE/dZ+c*dE/dT=-i*k*y(3);--------------------------A
here y(3) is the output from the system of 9 coupled differential equations solved by Ode45. the ode system is an pasted as under;
dydt(1)=-k1*y(3)+k1*y(7)+k5*y(9);
dydt(5)=-k2*y(6)+k2*y(
;%Pop of the metastable state
dydt(9)=k1*y(3)+k2*y(6)-k1*y(7)-k2*y(
-k5*y(9);
dydt(2)=-k3*y(2)-k2*y(3)+k1*y(
;
dydt(3)=-k1*y(1)-k2*y(2)-k4*y(3)+k1*y(9);
dydt(4)=-k3*y(4)-k1*y(6)+k2*y(7);
dydt(6)=-k1*y(4)-k2*y(5)+k2*y(9)-k6*y(6);
dydt(7)=k1*y(1)+k2*y(4)-k4*y(7)-k1*y(9);
dydt(
=k1*y(2)+k2*y(5)-k6*y(
;
now here, y(3) is a vector output from the call of ode45(@fname,t,s) where 'fname' file contains the code above. now what i want to do is;
to get the plot of E(z,t) as a function of space and time from the above equation A. let;
t=linspace(0,10,100);
z=linspace(0,10,100);
how would i obtain the plot of E(z,t) from the above equation?
??
i would be grateful for your kind help.
thanks alot
