Vektorprodukt

$$\vec{u} \times \vec{v} = \left[ \begin{array}{c} u_{y} v_{z} ... ..._{x} v_{z} \\ u_{x} v_{y} - u_{y} v_{x} \\ \end{array}\right]$$ (1.13)

void   VPR(double *v0,double *v1,double *erg){ //Vektorprodukt v0 x v1
   *erg     =    *(v0+1) * *(v1+2) - *(v0+2) * *(v1+1);
   *(erg+1) =    *(v0+2) * *(v1)   - *(v0)   * *(v1+2);
   *(erg+2) =    *(v0)   * *(v1+1) - *(v0+1) * *(v1)  ;
   return;
}