Solve Laplace Equation Using FDM (In Cartesian Coordinate)

소스코드

clear all;

d=1;
Nx=100;
dx=d/Nx;
V0=0;
V1=2;
N=5000;

V=zeros(Nx+1);
V(1)=V0;
V(Nx+1)=V1;

for n=1:1:N
    for i=2:1:Nx
        V(i) = (V(i+1)+V(i-1))/2;
    end
end

x=[0:dx:d];

figure();
plot(x,V);

이미지

(준비중...)