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);