一个轴对称原子核，采用柱坐标系，形变时保持体积不变，即始终为4/3（半径为1，圆周率已经考虑），并且已知ρ的平方与z的函数关系，满足体积守恒（用fortran，高斯勒让德积分，比较容易验证，轴对称，所以几何体垂直对称轴的截面都是圆，对应面积就是πρ²，积分就可以得到体积，程序也是如此，变换形状，始终有1.3333的值）
但是，在某个z值处算ρ²，在核形状变换过程中，常出现负值，始终查不出到底哪里出错，求助大佬，必要时我可以发具体代码，供大佬们验看