018 MATLAB Permossible Energy

Permossible Energy

syms me haba pusi a V E

n=3;
me=9.11*10^-31;
haba=1.054*10^-34;
a=10*10^-9;
E=(haba^2*pi^2*n^2)/(2*me*a^2);

f1=dsolve('D2y+(pi^2*n^2)/(a^2)*y=0','y(0)=0','x');

C(1)=0;
B(1)=1;
A(1)=0;
x=0;

k=0;

for i=linspace(0,10*10^-9,100);
 
    dx=(10*10^-9)/100;
 
    k=k+1;
      
    A(k)=((-2*me*E)/haba^2)*C(k);
 
    B(k+1)=B(k)+A(k)*dx;
 
    C(k+1)=C(k)+B(k)*dx;
 
    x(k+1)=x(k)+dx;
  
end

plot(x,C)