IDL编程学习之三重积分

15195775117

(积分区域好像是个半球)

积分函数:
function fxyz,x,y,z
return,z*(x^2+y^2+z^2)^1.5
end
中间积分限：
function pq_limits,x
return,[-sqrt(4.0-x^2),sqrt(4.0-x^2)]
end
最内积分限：
function uv_limits,x,y
return,[0.0,sqrt(4.0-x^2-y^2)]
end
pro excise
print,int_3d('fxyz',[-2.0,2.0],'pq_limits','uv_limits',6)
end
输出： 57.4177
注意int_3d中积分限的排列，是从外到内哒：
最外：[-2.0,2.0]
中间：'pq_limits'
最内：'uv_limits'