Plan Omega to Alpha - POA

Energy calculation function[]=energy(E) m=9.11*10e-19; haba=1.054e-34; k=0.021; dt=1e-11; x=[-5e-9:dt:5e-9]; a=1; C(1)=1; B(1)=1; for x=[-5e-9:dt:5e-9]         A(a)=(-2*m/haba^2)*(E-0.5*k*x^2)*C(a);     B(a+1)=B(a)+A(a)*dt;     C(a+1)=C(a)+B(a)*dt;     a=a+1; end EE=C(a) function[]=energy1(E0,E1) dE=1e-24; q=1; for E=E0:dE:E1     E3(q)=E;     Final(q)=energy(E);     q=q+1; end plot(E3,Final);

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)

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   ...