legendre7_demi_entier.mws

>   

>    restart;with(plots):

Warning, the name changecoords has been redefined

>    bool:=proc(expression)
piecewise(expression,1,0);
end:

     

[Maple OLE 2.0 Object]

>    LegendreQ_demi:=(n,z)->
(sqrt(Pi)/exp((n+1/2)*(arccosh(z)+ln(2)))/GAMMA(n+1))*
hypergeom([n+1/2,1/2],[n+1],exp(-2*arccosh(z)));

LegendreQ_demi := proc (n, z) options operator, arrow; sqrt(Pi)/exp((n+1/2)*(arccosh(z)+ln(2)))/GAMMA(n+1)*hypergeom([n+1/2, 1/2],[n+1],exp(-2*arccosh(z))) end proc

>    f:=(v,x,y)->(Re(LegendreQ_demi(v,x+I*y)),Im(LegendreQ_demi(v,x+I*y))):
g:=(v,p,x,y)->(Re(LegendreP(v,p,x+I*y)),Im(LegendreP(v,p,x+I*y))):

>    m:=(F,v,p,x,y)->sqrt(F(v,p,x,y)[1]^2+F(v,p,x,y)[2]^2):
phi:=(F,v,p,x,y)->arctan(F(v,p,x,y)[2]/F(v,p,x,y)[1])+bool(F(v,p,x,y)[2]<0)*Pi:
Kx:=-10..10; Ky:=-10..10;Kz:=-200..200;

Kx := -10 .. 10

Ky := -10 .. 10

Kz := -200 .. 200

>    animate3d([x,y,m(f,p+0.5,0,x,y)],x=Kx,y=Ky,p=0..8,frames=10,
                              axes=BOXED,labels=["x","y","mod(f)"]);

Animation : (x, y, module) des fonctions de Legendre de 2nde espèce ( n  = 5, 0 <  m < 8)

[Maple Plot]

>    nu:=40:

>    plot3d([x,y,m(g,0.5,0,x,y)],x=Kx,y=Ky,
                              axes=BOXED,labels=["x","y","mod(f)"],grid=[nu,nu]);

[Maple Plot]

>   

>    nu:=40:

>    plot3d([x,y,m(f,1,0,x,y)],x=Kx,y=Ky,                              axes=BOXED,labels=["x","y","mod(f)"],grid=[nu,nu]);

[Maple Plot]