📄 newu.m
字号:
function [u]=Newu(u,v,d,n)
A=500;B=500;D=500;C=200;
u0=0.01;lan=1.0000e-005;delta=0.1;
tmp=0;
for k=1:n
for k1=1:n
tmp=tmp+v(k,k1);
end% k1
end % k
for xk=1:n
for k=1:n
u(xk,k)=(1-lan)*u(xk,k);
for yk=1:n
if yk~=k
u(xk,k)=u(xk,k)-v(xk,yk)*A*lan;
end
if yk~=xk
u(xk,k)=u(xk,k)-v(yk,k)*B*lan;
if k==1
u(xk,k)=u(xk,k)-d(xk,yk)*(v(yk,2)+v(yk,n))*D*lan;
elseif k==n
u(xk,k)=u(xk,k)-d(xk,yk)*(v(yk,1)+v(yk,n-1))*D*lan;
else
u(xk,k)=u(xk,k)-d(xk,yk)*(v(yk,k+1)+v(yk,k-1))*D*lan;
end% k==1
end% yk~=xk
end% yk=1:n
u(xk,k)=u(xk,k)-(tmp-n)*C*lan;
end% k=1:n
end% xk=1:n
⌨️ 快捷键说明
复制代码
Ctrl + C
搜索代码
Ctrl + F
全屏模式
F11
切换主题
Ctrl + Shift + D
显示快捷键
?
增大字号
Ctrl + =
减小字号
Ctrl + -