LaguerreEig: siyyam1.m - File Exchange

Code covered by the BSD License

Highlights from
LaguerreEig

D=deriv(n)
P=param physiological values of the parameters
X=mult(n)
[Ps,Pu,L,V]=proj(A) calculates the projection operators, the eigenvalues
[x,A,f,kont]=jacob(x) calculates the critical points x and the corresponding Jacobian A
[x,coef,sol]=leuclag1(n,bas,P...
[x,w,T,iT,X,D]=lagtool(n,b) Tools for Laguerre spectral method
[x,w]=pd(n) Nodes and weights for Gauss-Laguerre quadrature
events_leuc( t,u )
f=leucemia3d(t,u)
laguer(a, x) http://dip.sun.ac.za/~weideman/research/mfiles/laguer.m
main_fermi
t=x2t(n,x)
airy_LagEig.m
anharmonic.m
hydrogen.m
laguerrestartup.m
main_leuc.m
marletta.m
pryce.m
siyyam.m
siyyam1.m
work2.m
View all files

from LaguerreEig by Damian Trif
Calculates numerical eigenvalues and eigenfunctions of Schroedinger problems on (0,infinity).

siyyam1.m

n=128;
[x,w]=pd(n);t=x2t(n,x);D=deriv(n);X=mult(n);
A=D^3+3*D^2+(3*eye(n)+2*X)*D+(2*eye(n)+2*X);
A(n-2,:)=ones(1,n);
A(n-1,:)=sum(D);A(n,:)=sum(D^2);
B=t'\(exp(-x).*((1-2*x).*sin(x)+cos(x)));
B(n-2)=1;B(n-1)=-1;B(n)=0;
c=A\B;
rez=(t'*c).*exp(x)-cos(x);
ind=find(x<10);
plot(x(ind),rez(ind));

Highlights from LaguerreEig

Highlights from
LaguerreEig