aias.m

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

function f=aias(X);
%> This function gives the asymptotics of Airy functions 
%> for positive  arguments, without the exponential factor.
%> Call: aias(X)
%> Input: X a row vector of positive arguments. 
%> For accuracy min(X) > 2
%> Output: a matrix of dimension 4 x dim(X) =
%> [ exp(-zet)*airy(0,X); exp(-zet)*airy(1,X);
%>   exp(+zet)*airy(2,X); exp(+zet)*airy(3.X)]
%> 	where zet = (2/3)*(X)^1.5
%> Reference Handbook of Mathematical Functions ch 10
%>
%>