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


