Robert M. Corless
Department of Applied Mathematics
University of Western Ontario
London, Canada

Copyright 2001 by Robert M. Corless
All rights reserved

Useful one-word commands 

Использование монокоманд

 Plotting in Maple 

Графика в Maple

Getting hardcopy of your plots 

Получение жёстких копий Ваших графиков

> restart;
> B := proc(v) local p,v0,u0;
if not type(v, numeric) then
'B(v)'
elif v < -1 then
(1.+v)^2/(1.+2*v)*B(v*(v+1.))
elif v=-1 then
0.
elif -1 < v and v < 0 then
v0 := v;
p := 1.;
while v0 < -0.1 do
p := p*(1.+v0)^2/(1+2.*v0);
v0 := v0*(v0+1.);
od:
p*Bseries(v0)
elif v = 0 then
1.
else
u0 := v;
p := 1.;
while u0 > 0.1 do
u0 := 2.*u0/(1.+(1.+4.*u0)^(1/2));
p := p*(1.+2.*u0)/(1+u0)^2;
od:
p*Bseries(u0)
fi
end:
> Bseries := v -> evalf(Sum('c(n)*v^(n-1)', n=1..infinity));

Bseries := proc (v) options operator, arrow; evalf(...
> c := proc( n ) option remember;
-1/(n-1)*add( binomial(n-i+1, i+1)*c(n-i), i=1..n-1 )
end proc;

> c(1) := 1;

c(1) := 1
> Bseries(0.05);

.9534459937
> B(0.05);

.9534459937
> B(1.);

.5665285893
> plotsetup( ps, plotoutput=`Beyn.ps`, plotoptions=`portrait,noborder`);
> plot( B, -2..2, -2..2, discont=true, colour=black );

С официального разрешения                    © 2002 Waterloo Maple, Inc

 
Hosted by uCoz