Scribd
Upload
282 views3 pages

Orthogonal Eigenfunction Hermitian

1) The document contains information about a student named Rohit enrolled in a BSc Physics course. It provides his roll number, subject and university roll number. 2) The problem is about orthogonal polynomials as eigenfunctions of Hermitian differential operators. 3) The solution involves calculating the third order Hermite polynomial and using it to verify the orthogonality relation between two Hermite polynomials of the same order.

Uploaded by

Alpha
Copyright
© © All Rights Reserved
We take content rights seriously. If you suspect this is your content, claim it here.
Available Formats
Download as PDF, TXT or read online on Scribd
282 views3 pages

Orthogonal Eigenfunction Hermitian

1) The document contains information about a student named Rohit enrolled in a BSc Physics course. It provides his roll number, subject and university roll number. 2) The problem is about orthogonal polynomials as eigenfunctions of Hermitian differential operators. 3) The solution involves calculating the third order Hermite polynomial and using it to verify the orthogonality relation between two Hermite polynomials of the same order.

Uploaded by

Alpha
Copyright
© © All Rights Reserved
We take content rights seriously. If you suspect this is your content, claim it here.
Available Formats
Download as PDF, TXT or read online on Scribd
You are on page 1/ 3

Name - Rohit

Course - BSc (H) Physics


Roll No - PHH19032
Subject - AMP Practical
Univ Roll No - 19048567003

Problem :-Orthogonal Polynomials as eigenfunctions of


Hermitian differential operators.
Solution:

INPUT:

clc;
funcprot(0)
function Y=H(n, X)
Y=0
for i=0:floor(n/2)
Y=Y+((factorial(n)*((-1)^i)*(2*X)^(n-2*i))/(factorial(i)*factorial(n-2*i)))
end
endfunction
n=input("Order of Hermite polynomial(n):")
X=poly(0,"X");
disp(H(n,X),"Hermite polynomial H(n,X)")
S= -2:0.0001:2
plot2d(S,H(n,S),5)
m=input("order of 2nd Hermite Polynomial for orthogonality
relation(m):")
z=integrate ('(exp(-(t)^2))*H(n,t)*H(m,t)','t',-100,100)
disp(z,"orthogonality relation of Hermite polynomials")
disp("Theotrical Value")
if n==m then
disp(factorial(n)*(2^n)*sqrt(%pi))
else
disp(0)
end
Output:

Order of Hermite polynomial(n):3


"Hermite polynomial H(n,X)"
-12X +8X³
order of 2nd Hermite Polynomial for orthogonality relation(m):3
85.077785
"orthogonality relation of Hermite polynomials"
"Theotrical Value"
85.077785

You might also like