exam2_solve.m

来自「解偏微分方程」· M 代码 · 共 18 行

function exam2_solve

t   = [0 10];
y0  = [4;0];
figure
ode45(@f,t,y0)


function dydt = f(t,y)

a = y(2);
b = (t/5-1)*y(2)-t^2*y(1)+1;
dydt=[a;b];
% dydt(1)=a;
% dydt(2)=b;