exemples/polaires/pol01.mp

   1 input gdd;
   2 input gdd-rep;
   3 input gdd-lbl;
   4 input gdd-tra;
   5 input svgnames;
   6
   7 lblPreambule := "\usepackage[frenchb]{babel}"       &
   8                 "\usepackage[garamond]{mathdesign}" &
   9                 "\usepackage{amsmath}";
  10
  11 input gdd-fct;
  12 input gdd-plt;
  13
  14 labeloffset := 6;
  15 gddU:=1.2cm;
  16
  17 vardef r(expr t) = 1/(cos(t)-cos(2*t)) enddef;
  18 vardef rp(expr t) = (r(t)*cos(t),r(t)*sin(t)) enddef;
  19 vardef f(expr t) = 3/2*t*t-7/18 enddef;
  20 def DroiteN(expr r, t) = Droite((r*cos(t),r*sin(t)), (r*cos(t)-2*sin(t),r*sin(t)+2*cos(t))) enddef;
  21 def traceDoubleVecteur(expr o,d)= drawdblarrow ((o-d)--(o+d)) gddEnPlace enddef;
  22
  23 def Crayon(expr t,c) = drawoptions(withpen pencircle scaled t withcolor c) enddef;
  24
  25
  26
  27 beginfig(1);
  28
  29 Repere(10,10,2,5,3.5,3.5);
  30 Axes;
  31 Debut;
  32 Graduations; Unites(1);
  33 trace Cercle(origine,1);
  34 Crayon(0.75,1.1*LightSlateGrey);
  35 trace rp(-arccos(1/4))--(0,0)--rp(arccos(1/4)) dashed evenly;
  36 Crayon(0.75,1.1*FireBrick);
  37 trace DroiteN(-2/3/sqrt(3),7*Pi/6);
  38 trace DroiteN(-2/3/sqrt(3),5*Pi/6);
  39 trace Representation(f,-3,3,100) rotated -90;
  40 Crayon(1.5,0.8*LightSlateGrey);
  41 trace CourbeEnPolaires(r,0.1,2*Pi/3-0.1,100);
  42 trace CourbeEnPolaires(r,-0.1,-2*Pi/3+0.1,100);
  43 trace CourbeEnPolaires(r,2*Pi/3+0.1,4*Pi/3-0.1,100);
  44 Crayon(1,DodgerBlue);
  45 traceDoubleVecteur(rp(Pi),(0,0.2));
  46 traceDoubleVecteur(rp(arccos(1/4)),0.05*(-sqrt(15),1));
  47 traceDoubleVecteur(rp(-arccos(1/4)),0.05*(sqrt(15),1));
  48 pointe rp(Pi);
  49 pointe rp(arccos(1/4));
  50 pointe rp(-arccos(1/4));
  51 Fin;
  52
  53 endfig;
  54
  55 end