comunque beta si dovrebbe calcolare così:
beta=acos((a*a+c*c-b*b)/2*a*c);