Code:
\documentclass[10pt,a4paper,twoside]{book}
\usepackage[latin1]{inputenc}
\usepackage{ngerman,a4wide}
\usepackage{amsmath,amsfonts,amssymb,textcomp}
\usepackage[amsmath,framed,thmmarks,thref]{ntheorem}
\usepackage{thmbox}
\usepackage{vector}
\usepackage{etex}
\usepackage{morefloats}
\usepackage{graphicx}
\usepackage{xkeyval}% Important to have it before
\usepackage{pst-pdf}
\usepackage{pst-all,pst-func,pst-xkey,pst-eucl,pst-math,pst-map3d,
pst-map3dII,pst-vue3d,pst-solides3d,pst-light3d,pst-slpe,infix-RPN,
pst-lens,pst-tree,pst-ob3d}
\usepackage{pst-3dplot}% Important to have it at the end!
\usepackage{epsfig}
\usepackage{floatflt} %Text around figures
\usepackage[german]{varioref}
\usepackage[margin=8pt,font=footnotesize,labelfont=bf]{caption}
\usepackage{ifthen}
\usepackage{pifont} %dinglist
\usepackage{fancyhdr}
\usepackage{fancybox}
\usepackage{multirow} %Array
\usepackage{array} %Array
\usepackage{enumerate}
\usepackage{ifsym}
\usepackage[normalem]{ulem} %Undeline Mode in Text
\usepackage{txfonts}
\usepackage{cancel}
\usepackage{caption}
\usepackage{datatool} %Data Tool / Arithmetic
\usepackage{intcalc} %Simple Arithmetic
\usepackage{tensor} %Tensors
\usepackage{braket} %Quantenmechanik
\usepackage{siunitx} %Units
%\usepackage{esvect} %Vectors
%
%
%
\begin{document}
\def\RadA{1}\def\RadB{2}
\psset[pst-solides3d]{lightsrc=100 20 40,
lightintensity=2,Decran=17}
\begin{pspicture}(-2,-1.8)(2.1,3)
%\psgrid[gridcolor=gray,subgridcolor=lightgray]
%Enable Memory
\psset{solidmemory}
%Ellipsoid as streched sphere
\psSolid[object=sphere,r=1,fillcolor=gray,grid=false,
ngrid=50 50,transform={1 1 2 scaleOpoint3d},
name=MySolid]
%Plane of Tangent (reset of strech before required)
\psSolid[object=plan,definition=solidface,args=MySolid 614,
base=-1.25 1.25 -0.3 0.3,transform={1 1 1 scaleOpoint3d},
fillcolor=gray,
name=MyPlan]%,showBase
%Definition of 3D Point in Center of Ellipsoid
\psSolid[object=point,args=0 0 0,action=none,
name=MyPO]
%Definition of 3D Point in Center of Plane of Tangent
\psSolid[object=point,definition=solidcentreface,
args=MySolid 614,action=none,
name=MyPF]
%Vector Omega and Line to center of ellipsoid
\psSolid[object=vecteur,args=MyPF,
transform={MyPF translatepoint3d},name=MyTMP]
\psSolid[object=line,args=MyPO MyPF,
linestyle=dashed,linecolor=darkgray]
%LABEL ?????????????????????????????????????????
\psSolid[object=point,definition=translatepoint3d,
args=MyPF MyTMP,action=none,name=A,text={$\omega$}]
%LABEL ?????????????????????????????????????????
%Definition of 3D Point as Projection of Center of
%Ellipsoid to Plane of the Tangent
\psSolid[object=point,definition=orthoprojplane3d,
args=MyPO MyPlan,action=none,
name=MyPS]
%Vector of Angualer Momentum and Line to center of ellipsoid
\psSolid[object=vecteur,args=MyPS,
transform={MyPS translatepoint3d},
name=MyVecI]
\psSolid[object=line,args=MyPO MyPS,
linestyle=dashed,linecolor=darkgray]
%Construct of vector in Plane and Vector Product (Normal Vector)
\psSolid[object=vecteur,definition=vecteur3d,
args=MyPF MyPS,action=none,
name=MyVecII]
\psSolid[object=vecteur,definition=vectprod3d,
args=MyVecII MyVecI,action=none,
name=MyNVecI]
%Define Plane
\psSolid[object=plan,definition=normalpoint,
args={MyPS [MyVecII MyNVecI]},action=none,
%base=-0.25 0 0 0.25,fillcolor=gray,
name=MyPlanI]
%Test
\psProjection[object=point,plan=MyPlanI,
args=-0.2 0,action=none,
name=P1]
\psProjection[object=point,plan=MyPlanI,
args=0 0,action=none,
name=P2]
\psProjection[object=point,plan=MyPlanI,
args=0 0.2,action=none,
name=P3]
\psProjection[object=rightangle,
linecolor=lightgray,args=P1 P2 P3]
%Draw AllComponents
\composeSolid
%Core System
\DTLadd{\valA}{\RadA}{0.8}
\DTLadd{\valB}{\RadB}{1.0}
\axesIIID[axisnames={\xi,\eta,\zeta},
arrowinset=0,arrowscale=1.5
](\RadA,\RadA,\RadB)(\valA,\valA,\valB)
\end{pspicture}
\end{document}
Lesezeichen