matlab中怎么编这个程序啊：-u"=4-6x(很急啊)

来源：百度知道 编辑：UC知道 时间：2024/06/25 01:30:54
-u"=4-6x in[0，1]
u(0)=u'(1)=0
分10，20，40，80，160个点五种情况求数值解，并画图和精确解比较，在k=1，2，3时算阶
就算只指导一部分也好……

%归一化：
%令u=y1,u'=y2,x=y3
%dy1/dt=y2
%dy2/dt=-4+6*y3
%dy3/dt=1
%函数文件：
function dy=fun(t,y)
dy=[y(2);4-6*y(3);1];
%脚本程序文件：
clear
ts=1/10*(0:10);%此处可改20，40，80，160
y0=[0;?;0];
[t,y]=ode45('fun',ts,y0);
u=y(:,1)
du=y(:,2)
x=t