编了个程序,得到了预期的多条x-y曲线,但是想得到各时间点的具体值。研究了好几天,还是做不出。 哪位高手能演示下? 随便输出哪一条曲线都可以。 程序如下:
function IV
tspan=[0, 30];
x0 = [0 0 0 0 0 100 0 0 0 0 0 0 0 0 0 0 100];
[T, X] = ode45(@lun_fun, tspan, x0);
plot(T, X)
function dx = lun_fun(t, x)
dx = zeros(17, 1);
k12=0.69;
k13=4.6;
k21=0.35;
k26=0.33;
k27=0.67;
k34=0.66;
k30=6.6;
k43=0.33;
k46=0.17;
k47=0.33;
k56=0.01;
k62=0.67;
k64=0.33;
k65=0.005;
k68=0.17;
k69=1.7;
k610=10;
k611=1.7;
k612=2.2;
k76=0.023;
k86=0.083;
k96=0.023;
k100=0.0046;
k103=0.023;
k106=0.023;
k110&#