exam2_solve.m

解偏微分方程

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;