200# {version 1.95-beta-10} Cmd [Fenetre(-4+4*i,4-4*i,1.125+1.125*i), Marges(0,0,0,0), Border(0)]; [OriginalCoord(1),IdMatrix()]; [theta:=0.1745, phi:=0.9599, IdMatrix3D(), ModelView(ortho)]; [GrayScale(0), ComptGraph:=5]; Mac grandCercle = [ SaveAttr(), $O:=%1, $R:=%2, $A:=%3, $B:=%4, $U:=Prodvec(A-O,B-O), $N:=normalize(Prodvec(\n,U)), $A':=O+R*N, $B':=O-R*N, $V:=Prodvec(A'-O,\n), LineStyle:=solid, Arc3D(A',O,B',R,1,U), LineStyle:=dotted, Arc3D(A',O,B',R,-1,U), if StrComp(%5,"")=0 then $v:=R*normalize((A+B)/2), LabelDot(Proj3D(O+v),%5,%6) fi, RestoreAttr() ]; dsphere = [ SaveAttr(), FillStyle:=none, LineStyle:=solid, $np:=%3, $R:=%2, $O:=%1, $nm:=%4, $C:=O+R*vecK, $h:=2*R/(1+np), $H:=R, $R2:=R^2, for $k from 1 to np do Inc(C,-h*vecK), Inc(H,-h), $r:=sqrt(R2-H^2), Cercle3D(C,r,vecK), od, $alpha:=0, $dalpha:=2*pi/(nm-1), for k from 1 to nm do Inc(alpha,dalpha), Cercle3D(O, R, Prodvec([R*exp(i*alpha),0],vecK)) od, RestoreAttr() ]; Cmd [tMin:=-3.7566,tMax:=3.7566,DotScale:=1+i]; Graph objet1 = [ view(-4,4,-4,4), Marges(0,0,0,0), size(9), O:=[0,0], N:=3*vecK, S:=-3*vecK, A:=3*normalize(M(1,1,0)), B:=3*normalize(M(1,3,0)), Color:=lightgray, Width:=2, dsphere([0,0],3,20,20), Color:=black, Width:=8, grandCercle( O, 3, A, B), LineStyle:=dotted, Ligne3D( [N,S], 0), LabelDot3D(N,"$P_1$","N",1), LabelDot3D(S,"$P_2$","S",1), LabelDot3D(O,"$O$","NO",1) ];