029 Matlab Energy calculation

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);