偏微分方程数值解实验报告.docxVIP

偏微分方程数值解 上机实验报告 （一）实验一 一、 上机题目： 用线性元求解下列边值问题的数值解： -??′′ + ??2 4 ?? = ??2sin?? ??，0x1 2 2 y(0)=0，??′(1) = 0 二、 实验程序： function S=bz x=fzero(@zfun,1); [t y]=ode45(@odefun,[0 1],[0 x]); S.t=t; S.y=y; plot(t,y) xlabel(x:′从0一直到1) ylabel(y) title(线性元求解边值问题的数值解) function dy=odefun(x,y) dy=[0 0]; dy(1)=y(2