Differential equations - recursive anonymous functions

[genehacker]

So,

    I have to tell you this forum is great for octave users who seeks help
and advise. thanks very much guys.

So i have this differential equations setup and a couple of recursive
anonymous functions in there.

***********************
g = @(t,x2,x4);
x4 = @(t,x2,g);
x5 = @(x1)   k1max- ((k1max- k1min)/(1 + (f1 * x1/IC2)));
x6 = @(x1)   k2max- ((k2max- k2min)/(1 + (f1 * x1/IC2)));
x7 = @(x1)   k23min- ((k23min- k23max)/(1 + (f1 * x1/IC2)));

xdot = @(x,t) [ ( g * x4 / pre_OC) - x5(x(1)) * x(1);
                        x6(x(1)) * pre_OB - x7(x(1)) * x(2);
                        x7(x(1)) * x(2) - k3 * x(3)];

x = [0,0,0; 0,0,0];
time = linspace(1,50,51);
for i=1:length(time)-1
        seedtime = time(i:i+1);
        [x,istate,msg] = lsode(xdot,[0 0 0],seedtime);
endfor

**************
so that I can run differential equations step by step to control the
time-dependent parameterrs g,x4,x5,x6 and x7. The problem Im having is
because g(x4) and x4(g), how should I write the xdot(1) equation. simply
doing this gives error
xdot(1) = g * x4 / pre_OC) - x5(x(1)) * x(1);
Note that g(0)=0; so system can compute other parameters from it and start
running.

How should I write xdot? Also if I want to call these parameters from within
the for loop at each run, how should I call them?

thanks a bunch,
k3.
-- 
View this message in context: 
http://www.nabble.com/Differential-equations---recursive-anonymous-functions-tp20080482p20080482.html
Sent from the Octave - General mailing list archive at Nabble.com.