function cal_chalf(n,delta,r)
Sqrt(4.0*Power(-1.0 + n,2.0)*Power(r,2.0) + Power(delta,2.0)*Power(-1.0 + Power(r,2.0),2.0))/Sqrt(4.0*Power(n,2.0)*Power(r,2.0) + Power(delta,2.0)*Power(-1.0 + Power(r,2.0),2.0))
end
