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:
Model syllabus (+3): Downoad
Get everything about your B.Sc (+3) Physics Practical⇠ ⇠Click
Post a Comment