Solve the s-wave radial Schrodinger equation for a particle of mass, with the anharmonic potential v. find the GS & FES using scilab

 Read all Books:  CLICK

CODING:

clc;
clear();
h=197.3;
m=941;
d=0.755501;
k=110;
b=0,15,30;
a=1.44;
q0=3.795;
r0=0.131349;
rm=10;
n=500;
r=linspace(r0,rm,n)
d=(r-r0)/r
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
A(i,[i-1:i+1])=[1,-2,1]
for i=1:n
    V(i)=d*((exp(-2*a*(r(i)-r0)/r(i))))
    -exp(-a*(r(i)-r0)/r(i))
end
H=-(((h^2)/(2*m*d*d))*A)+V;
[R,y]=spec(H)
w=find(H>0.0);
disp(y(3,3))
disp("Energy of first excited state is:")
disp(H(w(2)))
R(1,:)=R(:,1)
R(2,:)=R(:,2)
plot(r,(R(1,:)),'red')
plot(r,(R(2,:)),'green')
xgrid()

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