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)