//Bogenmaß (rad) zwischen punkt a und punkt b im 3d koordinatensystem?
double top = a[0]*b[0]+a[1]*b[1]+a[2]*b[2]; //Skalarprodukt beider vektoren
double bot_left = sqrt((pow_d(a[0])+pow_d(a[1])+pow_d(a[2]))); //Betrag vektor (Raumdiagonale ?) 1
double bot_right = sqrt((pow_d(b[0])+pow_d(b[1])+pow_d(b[2]))); //Betrag vektor (Raumdiagonale ?) 2
return acos (top / (bot_left * bot_right)) * 180.0 / PI; //in Bogenmaß umrechnen