let
mhalfpi = asin (-1.)