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 ProgrammingC ++ ProgrammingModel syllabus (+3): DownoadList of all practicals:   DownloadGet all Books:  DownloadPhysics Notes: CLICKFor online classes & notes, Contact: quantaphysics3@gmail.com                                          Join Telegram: quantaphysics.com                                    Subscribe Youtube: BS quantaphysics NOTES:1. Electric Field Intensity ðŸ‘‰ Click2. Coulomb's Law ðŸ‘‰ Click3. IIT-JAM Question Sets ðŸ‘‰ Click4. Blackbody Radiation ðŸ‘‰ Click5. Inadequacies of Classical Mechanics ðŸ‘‰  Click6. Constraints  ðŸ‘‰ ClickFor online classes & notes, Contact: quantaphysics3@gmail.com                                          Join Telegram: quantaphysics.com                                    Subscribe Youtube: BS quantaphysics                                                            How to Prepare for IIT-JAM  ⇠ ⇠ Click HereHow to prepare for JEST ⇠ ⇠ Click HereGet everything about your B.Sc (+3) Physics Practical⇠ ⇠Click`