DiBa := proc (phi) options operator, arrow; ((-1+1/h*H)*k/Pi*So+Ba(phi)^2)^(1/ 2) end proc;