質問: odeset('refine',10) sol=ode45(@(t,y) myfun(t,y),[0 50],[0 10]) % グラフ化 plot(sol.x, sol.y(1,:)) xlabel('t') ylabel('振れ角θ') title('時間,s')
回答: odeset: https://jp.mathworks.com/help/matlab/ref/odeset.html の refine の項目を見ると、以下の記載がありました。 refine は length(tspan) > 2 の場合や ode ソルバーが解を構造体として返す場合には適用されません。 ですので、構造体ではなく数値を返すよう ode45 を実行して見てください。 例: [t,y] = ode45(@(t,y) myfun(t,y),[0 50],[0 10])