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