Comandos Úteis e Simples

Seção 2.7: Plotando no Maple

Session 2.7.5: Fazendo cópias dos seus 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));

>	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;

>	Bseries(0.05);

>

B(0.05);

>

B(1.);

>	plotsetup( ps, plotoutput=`Beyn.ps`, plotoptions=`portrait,noborder`);

>	plot( B, -2..2, -2..2, discont=true, colour=black );

>

>