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