271.html

Manual Maple 2.7.1.mws

Capítulo 2: Comandos Úteis e Simples

Seção 2.7: Plotando no Maple

Session 2.7.1: Plots bi-dimensionais

> restart;

> plot( Zeta(t), t=-3..3, y=-3..3, discont=true, colour=black );

> alias( W=LambertW );

> plot( {W(x),W(-1,x)}, x=-0.5..1.5, y=-4..1, colour=black );

> plot( [t*exp(t), t, t=-4..1], view=[-0.5..1.5,-4..1], colour=black );

> plot( {seq(orthopoly[T](k,x),k=0..20)}, x=-1..1, y=-1..1, axes=BOXED, colour=black, scaling=CONSTRAINED, numpoints=101 );

> with( orthopoly ):

>    plot( {seq( [T(2,t),T(k,t),t=-1..1], k=1..5) ,
       seq( [T(2,t),-T(k,t),t=-1..1],k=1..5)},            view=[-1..1,-1..1], colour=black,
      scaling=CONSTRAINED, axes=BOXED );

>    plot( {seq( [T(2,t),T(k,t),t=-1..1], k=1..20) ,
       seq( [T(2,t),-T(k,t),t=-1..1],k=1..20)},            view=[-1..1,-1..1], colour=black,
      scaling=CONSTRAINED, axes=BOXED );

> restart;

> with(student):

> leftbox(1/(1+t), t=0..1, 6, colour=black );

>

> restart;

> f := 1/(2+sin(theta));

> a := n -> Int( f*cos(n*theta), theta=0..2*Pi )/Pi;

> b := n -> Int( f*sin(n*theta), theta=0..2*Pi )/Pi;

> Approx := N -> add( b(k)*sin(k*theta), k=1..N )
+add( a(k)*cos(k*theta), k=1..N ) + a(0)/2;

> err := evalf( f - Approx(5) ):

> plot( err, theta=0..2*Pi );

> restart;

> G := t -> frac(1/t);

$G := proc (t) options operator, arrow; frac(1/t) end proc$

> plot( G, 0..1, 0..1, numpoints=101, colour=black );

> y := [ seq( k/10., k=0..10 ) ];

> y[11] := y[11] - Float(1,3-Digits);

> t := Vector[row]( 1..1100 ):

>    for i to 100 do
  for j to 11 do
    t[(i-1)*11+j] := 1./(i+y[j]);
  end do;
end do;

> gt := map( G, t ):

> gt_0 := map( G, t[1100-7..1100] );

> pts := zip( (x,y)->[x,y], t, gt ):

> whattype(pts);

> plot( convert(pts,list), view=[0..1,0..1], colour=black );

> pieces := array(1..100):

>    for i to 100 do
  pieces[i] := NULL;
  for j to 11 do
    pieces[i] := pieces[i], 1./(1+y[j]), y[j];
  end do;
  pieces[i] := [pieces[i]];
end do:

> pieces := array(1..100):

> ys := [seq(y[j],j=1..11)]:

> for i to 100 do
pieces[i] := zip((i,j)->[i,j],[seq(1./(i+y[j]),j=1..11)],ys);
end do:

> plot( {seq(pieces[k],k=1..100)}, view=[0..1,0..1], colour=black);

> restart;

> with( plots, polarplot );

> polarplot( {seq(cos(5*theta)+n*cos(theta),n=-5..5)},
theta=0..2*Pi, colour=black, axes=BOXED, scaling=CONSTRAINED);

> r := (n,theta) -> cos(5*theta) + n*cos(theta);

> plots[display]( [seq(plot( [r(n,t)*cos(t),r(n,t)*sin(t),t=0..2*Pi]) ,n=-5..5)],insequence=true);

> polarplot( (4*cos(3*theta)+cos(13*theta))/cos(theta), theta=0..2*Pi, colour=black, axes=BOXED );

> polarplot( (4*cos(theta)+cos(9*theta))/cos(theta), theta=0..2*Pi, colour=black, axes=BOXED );

> polarplot( exp(cos(theta))-2*cos(4*theta)+sin(theta/12)^5, theta=0..24*Pi, colour=black, axes=BOXED );

>