matlab计算抛物型偏微分方程的数值解

更新时间:2023-06-08 22:33:45 阅读: 评论:0

function uxy = Ini2U(x,y)
format long;
孙红娟灭鼠工作总结uxy = exp(-10*x*x-10*y*y);
爱我中华手抄报function u = peParabImp(c,dt,n,minx,maxx,lbu,rbu,M)
format long;
maxx=100;
minx=1;
n=10
h = (maxx-minx)/(n-1);
新生儿笑是为什么%u0(1) = lbu;
u0(1) =0;
%u0(n) = rbu;
u0(n) =0;
for j=2:n-1
    u0(j) = PrIniU(minx+(j-1)*h);
end
it经理怎么治疗老花眼u1 = u0;
M=100
dt=0.01
c=0.1
lbu=0
rbu=0
中国僵尸片
for k=1:M
    A = zeros(n-2,n-2);
    cb =  - transpo(u0(2:(n-1)));
    cb(1) = cb(1) - dt*c*lbu/h/h;
    cb(n-2) = cb(n-2) - dt*c*rbu/h/h;
    A(1,1) = -2*dt*c/h/h -1;
    A(1,2) = dt*c/h/h ;
    for i=2:n-3
        A(i,i-1) = dt*c/h/h ;
        A(i,i) = - 2*dt*c/h/h -1 ;
        A(i,i+1) = dt*c/h/h  ;
    end
    A(n-2,n-2) = -2*dt*c/h/h -1;
    A(n-2,n-3) = dt*c/h/h;
    u1(2:(n-1)) = A\cb;
养阴清肺膏
    u0 = u1;
end
u = u1;
format short;

本文发布于:2023-06-08 22:33:45,感谢您对本站的认可!

本文链接:https://www.wtabcd.cn/fanwen/fan/82/906116.html

版权声明:本站内容均来自互联网,仅供演示用,请勿用于商业和其他非法用途。如果侵犯了您的权益请与我们联系,我们将在24小时内删除。

标签:灭鼠   老花眼   僵尸   总结   治疗
相关文章
留言与评论(共有 0 条评论)
   
验证码:
推荐文章
排行榜
Copyright ©2019-2022 Comsenz Inc.Powered by © 专利检索| 网站地图