حل معادله پواسون یک بعدی با متلب

clc
clear all
n=200;
l=5E-6;
ro=2E-12;
e0=3.4;
x=linspace(0,l,n); dx=x(2)-x(1);

v1=10;
v2=0;

V=zeros(1,n)+(v1+v2)/2;
V(1)=v1;V(n)=v2;
VV=V;

error=2;
j=1;
while(error>1E-5)
    VV=V;
    for i=2:n-1
        V(i)=(-2*dx^2*ro/e0+V(i+1)+V(i-1))/2.;
    end
    error=max(abs(V-VV));
    j=j+1
end

plot(x/1E-6,V,'k','LineWidth',2)
xlabel(' X (*10^{-6}m) ','FontSize',14,'FontWeight','bold',...
    'color','r','FontName','Times New Roman')
ylabel('V','FontSize',14,'FontWeight','bold',...
    'color','r','FontName','Times New Roman')