function B = angle(a,b,c)
B = asin((b-c)/a);