1

odeplot を使用したいので、後で結果をプロットするのではなく、段階的に結果を取得します。こんな感じで書いてみたのですが上手くいかないので、どなたか教えていただければ幸いです。

%Parameters

s = 1; 

q = 1; 

w = 0.1610;

y0 = [30 1 30]; % Initial values

tspan = [0 10]; % Time 0<t<10

plot=odeset('OutputFcn','odeplot');

[t, y] = ode45(@(t,y) concentration(t, y, s, q, w), plot, tspan, y0);
4

1 に答える 1