double VWI(double *v0,double *v1){ //Winkel zwischen v0 und v1 double erg; erg = SPR(v0,v1)/(VBG(v0)*VBG(v1)); if(erg<-1.0){ erg=-1.0; } if(erg> 1.0){ erg= 1.0; } return acos(erg); };