gegenb.m

求解量子力学的薛定谔方程

function f=gegenb(a,n);
%> This file calculates the coefficients of the Gegenbauer 
%> polynomial C(a,n), HMF Chapter 22
%> Call gegenb(a,n),  a > 0 , n >= 0 integer
%> Output: (n+1)x(n+1) matrix with integer entries =
%> the coefficients of the n+1 first Gegenbauer polynomials,
%> i.e. the polynomials up to order n inclusive,
%> the k:th row = coefficients of the k:th polynomial
%> To get the polynomials call evalpol(f,x);

%>