{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 704 "p:=-0.3:\nn:=60:\nni:=12:\ng:=NULL :\n\nfor j from 1 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 t o ni do\n m:=(aa+bb)/2:\n if f(aa)*f(m)<0 then bb:=m else aa:=m fi:\no d:L:=L,m od: \nL:=[L]:\nf1:=plot([[seq([1+op(i,L)*cos((i-1)*Pi/n+0.00 1),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)],[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.00 1)],i=1..n+1)]],colour=clo[j]):\ng:=g,display([f1]):\nod:\n" }}} {EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 468 "for j from 5 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=c lo[j]):\ng:=g,display([f1]):\nod:" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 101 "display([g,plot([[1,0],[-1,0]],style=point,symbol=ci rcle,color=blue)],axes=none,scaling=constrained);" }{TEXT -1 0 "" }}}} {MARK "4 0 1" 0 }{VIEWOPTS 1 1 0 1 1 1803 }