作业帮 > 综合 > 作业

matlab du/dt=d(du)/dx^2 x属于(0,1),t属于(0,T]u(0,t)=u(1,t)=0u(x,

来源:学生作业帮 编辑:拍题作业网作业帮 分类:综合作业 时间:2024/05/01 15:43:30
matlab
du/dt=d(du)/dx^2 x属于(0,1),t属于(0,T]
u(0,t)=u(1,t)=0
u(x,0)=0
若在u(1/2,0)处有一个扰动1/(2^10),取h=1/16,r=1/2
分别用古典显格式和richardson格式计算8层
1)打印出第8层上各结点处的计算值
2)预测继续算下去计算值的变化趋势
3)分析上述趋势产生的原因
要用matlab编程
望告知具体程序~55555自己做不来啊………= =|||都没有人会做?分数已经到上限了
#include
#include
#include
void main()
{ double u[16][16],x[16];
double h=0.0625,r=0.5,y;
int a=1,i,j;
y=r*h*h/a;
for(i=0;i