Solve the s-wave Schrodinger equation for ground state hydrogen atom using scilab(5th sem practical)

 Read all Books:  CLICK

CODING:

clc;
clear();
a=6,8,10
h=627.783
m=0.521e6
n=200
r0=1e-15
rm=10
e=3.805
r=linspace(r0,rm,n)
d=(rm-r0)/n
A=zeros(n,n)
V=zeros(n,n)
A(1,[1:2])=[-2,1]
A(n,[n-1:n])=[1,-2]
for i=2:n-1
    A(i,[i-1:i+1])=[1,-2,1];
end
for i=1:n
    V(i,i)=((e.^2)/r(i))*(exp(-r(i)/a));
    V1(i)=(-(e.^2)/r(i))*(exp(-r(i)/a));
end
H=((h.^2)/(2*m*d*d)*A)+V
[y,eig]=spec(H)
subplot(2,2,1)
for i=10:n
    plot(r(i),V1(i),'y+')
    xlabel("X-Axis")
    ylabel("Y-Axis")
end
disp(eig(n-1,n-1))
disp(eig(n-2,n-2))
disp(eig(n-3,n-3))
subplot(2,2,2)
plot(r,abs(y(:,n-1)),'r')
xlabel("X-Axis")
ylabel("Y-Axis")
subplot(2,2,3)

plot(r,abs(y(:,n-2)),'r')
xlabel("X-Axis")
ylabel("Y-Axis")
subplot(2,2,4)

plot(r,abs(y(:,n-3)),'g')
xlabel("X-Axis")
ylabel("Y-Axis")
subplot(2,2,5)
legend('y, Ground State n=1, y2 First Excited State n=2,1,0')

OUTPUT:


Scilab Programming
C ++ Programming
Model syllabus (+3): Downoad
List of all practicals:   Download
Get all Books:  Download
Physics Notes: CLICK
For online classes & notes, Contact: quantaphysics3@gmail.com
                                          Join Telegram: quantaphysics.com
                                    Subscribe Youtube: BS quantaphysics 
NOTES:
1. Electric Field Intensity 👉 Click
2. Coulomb's Law 👉 Click
3. IIT-JAM Question Sets 👉 Click
4. Blackbody Radiation 👉 Click
5. Inadequacies of Classical Mechanics 👉  Click
6. Constraints  👉 Click

For online classes & notes, Contact: quantaphysics3@gmail.com
                                          Join Telegram: quantaphysics.com
                                    Subscribe Youtube: BS quantaphysics 
                                                           
How to Prepare for IIT-JAM  ⇠ ⇠ Click Here
How to prepare for JEST ⇠ ⇠ Click Here
Get everything about your B.Sc (+3) Physics Practical⇠ ⇠Click