parabcyl.m

来自「用于计算脉冲激光照射下物体的温度分布」· M 代码 · 共 39 行

function f=Parabcyl(x,r,N)
% computer parabolic cylinder funcion
%    input:x,r,N
%    output the value of f
%
if nargin<3 N=100;end

n=length(x);

a0=1;
a(1)=1;
ac=-r-0.5;
a(2)=ac;
a(3)=ac;

i=4;
a(4)=ac*a(i-2)+(1/4)*(i-2)*(i-3)*a(1);
for i=5:N,
    a(i)=ac*a(i-2)+(1/4)*(i-2)*(i-3)*a(i-4);
end

%for i=1;n,
i=1;
m=0;
ya(m/2+1)=a0*x(i)^(m);
for m=2:2:N,
    ya(m/2+1)=a(m)*x(i)^(m)/prod([1:m]);
end
y1=sum(ya);

for m=1:2:N,
    yb((m+1)/2)=a(m)*x(i)^m/prod([1:m]);
end
y2=sum(yb);

f=cos(pi*r/2)/(pi)^0.5*gamma(0.5+0.5*r)/(2^(-r/2))*y1+sin(pi*r/2)/(pi)^0.5*gamma(1+0.5*r)/(2^(-r/2-0.5))*y2;

%end
%f=F;