017 MATLAB reasonable energy

reasonable energy

syms me haba pusi a V E

p=0;
final(1)=0;

for E=linspace(0,10^-20,100)
  
    p=p+1;
  
    me=9.11*10^-31;
    haba=1.054*10^-34;
    a=10*10^-9;
    k=0;
    C(1)=0;
    B(1)=1;
    A(1)=0;
    x=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
  
    finalx(p)=C(101);
    E1(p)=E;  
  
end
       
plot(E1,finalx)
grid on

for o=1:101;
    if finalx(o)<=0;
        E2(1)=E1(o);
        break
    end
end

for m=o:101;
    if finalx(m)>=0;
        E2(2)=E1(m);
        break
    end
end

for n=m:101;
    if finalx(n)<=0;
        E2(3)=E1(n);
        break
    end
end
A=E2