2017年3月16日
Matlab:Crank Nicolson方法求解线性抛物方程
摘要： 1 tic; 2 clear 3 clc 4 M=[10,20,40,80,160,320,640];%x的步数 5 K=M; %时间t的步数 6 for p=1:length(M) 7 hx=1/M(p); 8 ht=1/K(p); 9 r=ht/hx^2; %网格比 10 x=0:hx:1; 11 t=0:ht:1; 12 numerical=zeros(M(p)+1,K... 阅读全文
posted @ 2017-03-16 17:37 胡冬冬 阅读(6495) 评论(0) 推荐(1) 编辑