IDL的积分运算与MATLAB有相似之处,都是需要先定义一个函数。 [part_1:二重积分例子] ;积分函数:
function f1,x,y
return,y*cos(x^5)
end
;内部积分限
function f2,x
return,[0.0,x^2]
end
pro int_example
print,int_2d('f1',[0.0,2.0],'f2',48);这个48是固定的
end 输出:0.0551419
二重积分练习: 极坐标求圆的面积 代码:(用t表示角度) function f1,t,r→这里t与r的顺序不能错,t是外部积分变量,r是内部积分变量,顺序变了结果不对
return,r
end
function f2,r
return,[0.0,1.0]
end
pro linshi
print,int_2d('f1',[0.0,2.0*!pi],'f2',6)
end 输出:3.14159