jacobidiedai.m

来自「在matlab环境下对偏微分方程的Jacobi迭代解法的算法实现」· M 代码 · 共 14 行

function [u,u0]=jacobidiedai(f,m,n,h1,h2,u,x,y)
for i=2:m
    for j=2:n
        u0(i,j)=(f(x(i),y(j))+(u(i,j-1)+u(i,j+1))/(h2*h2)+(u(i-1,j)+u(i+1,j))/(h1*h1))/(2*(1/(h1*h1)+1/(h2*h2)));
    end
end
for j=1:n+1
    u0(1,j)=sin(y(j))+cos(y(j));
    u0(m+1,j)=exp(1)*(sin(y(j))+cos(y(j)));
end
for i=1:m+1
    u0(i,1)=exp(x(i));
    u0(i,n+1)=exp(x(i))*(sin(1)+cos(1));
end