page 1 page 2 page 3 page 4 page 5 page 6 page 7 page 8

> curlA:=ContractIndices(DA,epsilon,[[1,2],[2,3]]);


2.3 Obtaining the graphics

> restart;

> addcoords(conic,[r,mu,lambda],[r*JacobiDN(lambda,sin(3))*JacobiSN(mu,cos(3))
,r*JacobiSN(lambda,sin(3))*JacobiDN(mu,cos(3)),r*JacobiCN(lambda,sin(3))*JacobiCN(mu,cos(3))]);
> fsolve(jacobiSN(lambda,sin(3))=0,lambda, lambda=-4..-3 );#plane y=0

-3.157411419

> fsolve(JacobiCN(mu,cos(3))=0,mu, mu=-4..-3);#plane z=0

-3.356234394

> plot3d(2,mu=-3.356234394..3.356234394,lambda=-3.157411419..3.157411419,
coords=conic,axes=boxed,labels=[x, y, z]);
 
 

7