{VERSION 3 0 "IBM INTEL NT" "3.0" } {USTYLETAB {CSTYLE "Maple Input" -1 0 "Courier" 0 1 255 0 0 1 0 1 0 0 1 0 0 0 0 }{PSTYLE "Normal" -1 0 1 {CSTYLE "" -1 -1 "" 0 1 0 0 0 0 0 0 0 0 0 0 0 0 0 }0 0 0 -1 -1 -1 0 0 0 0 0 0 -1 0 }} {SECT 0 {EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 8 "restart:" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 92 "with(plots):\nclo:=seq(COLOR(RGB,i/ 9,0,(9-i)/9),i=0..9),seq(COLOR(RGB,(9-i)/9,i/9,0),i=1..8):" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 701 "p:=2:\nn:=60:\nni:=12:\ng:=NULL:\n \nfor j from 5 to 4 do b:=j/4:\n#Calcul de AM en fonction de alpha\nL: =NULL:\nfor k from 0 to n do alpha:=k*3.14/n+0.001:f:=unapply((t**2+4* t*cos(alpha)+4)**(p/2)-(2*b**p-t**p),t);\naa:=0.01:bb:=1:\nfor i to ni do\n m:=(aa+bb)/2:\n if f(aa)*f(m)<0 then bb:=m else aa:=m fi:\nod:L: =L,m od: \nL:=[L]:\nf1:=plot([[seq([1+op(i,L)*cos((i-1)*Pi/n+0.001),o p(i,L)*sin((i-1)*Pi/n+0.001)],i=1..n+1)],[seq([(1+op(i,L)*cos((i-1)*Pi /n+0.001)),-op(i,L)*sin((i-1)*Pi/n+0.001)],i=1..n+1)],[seq([-(1+op(i,L )*cos((i-1)*Pi/n+0.001)),op(i,L)*sin((i-1)*Pi/n+0.001)],i=1..n+1)],[se q([-(1+op(i,L)*cos((i-1)*Pi/n+0.001)),-op(i,L)*sin((i-1)*Pi/n+0.001)], i=1..n+1)]],colour=clo[j]):\ng:=g,display([f1]):\nod:\n" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 468 "for j from 6 to 13 do b:=j/4:\nL:= NULL:\nfor k from 0 to n do alpha:=k*3.14/n+0.001:f:=unapply((t**2+4*t *cos(alpha)+4)**(p/2)-(2*b**p-t**p),t);\naa:=0.1:bb:=10:\nfor i to ni \+ do\n m:=(aa+bb)/2:\n if f(aa)*f(m)>0 then aa:=m else bb:=m fi:\nod:L:= L,m od: \nL:=[L]:\nf1:=plot([[seq([1+op(i,L)*cos((i-1)*Pi/n+0.001),op (i,L)*sin((i-1)*Pi/n+0.001)],i=1..n+1)],[seq([(1+op(i,L)*cos((i-1)*Pi/ n+0.001)),-op(i,L)*sin((i-1)*Pi/n+0.001)],i=1..n+1)]],colour=clo[j]): \ng:=g,display([f1]):\nod:" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 101 "display([g,plot([[1,0],[-1,0]],style=point,symbol=circle,color=bl ue)],axes=none,scaling=constrained);" }{TEXT -1 0 "" }}}}{MARK "4 0 1 " 0 }{VIEWOPTS 1 1 0 1 1 1803 }