ich möchte die Anisotropie eines Einkristalls mittels des E-Moduls im dreidimensionalen Raum in Abhängigkeit der Kristallrichtungen darstellen.
Dazu soll zunächst die Elastizitätsmatrix von Silizium genutzt werden:
$\boldsymbol{C} =\begin{bmatrix}C_{11} & C_{12} & C_{12} &0&0&0 \\C_{12} & C_{11} & C_{12} &0&0&0 \\C_{12} & C_{12} & C_{11} &0&0&0 \\0&0&0&C_{44}&0&0\\0&0&0&0&C_{44}&0\\0&0&0&0&0&C_{44}\\\end{bmatrix}$
$C_{11} = 165,7\,\mathrm{GPa}, \; C_{12} = 63,9\,\mathrm{GPa}, \; C_{44} = 79,6\,\mathrm{GPa}$

Zur Umsetzung bin ich auf das Paket 3dplot gestoßen. (http://www.ctan.org/tex-archive/graphic ... kz-3dplot/) Dazu hangele ich mich an folgendem Minimalbeispiel entlang.
\documentclass{minimal}\usepackage{verbatim}\usepackage{tikz}\usepackage{3dplot}\begin{document}\tdplotsetmaincoords{70}{135}\begin{tikzpicture}[scale=1,line join=bevel,tdplot_main_coords, fill opacity=.7]\pgfsetlinewidth{.1pt}\tdplotsphericalsurfaceplot[parametricfill]{72}{36}%{sqrt(15/2)*sin(\tdplottheta)*cos(\tdplottheta)}{black}{\tdplotphi}%{\draw[color=black,-latex] (0,0,0) -- (2,0,0) node[anchor=north east]{$x$};}%{\draw[color=red,-latex] (0,0,0) -- (0,2,0) node[anchor=north west]{$y$};}%{\draw[color=black,-latex] (0,0,0) -- (0,0,2) node[anchor=south]{$z$};}%\end{tikzpicture}\end{document}
Kann hier jemand einen Tip geben bzw. gibt es alternative Lösungsmöglichkeiten?
Viele Grüße
Dirk