IDL编程学习之解非线性方程(组)

15195775117

IDL中解非线性方程(组)有4个函数，解非线性方程组有Broden法和牛顿法，解一元非线性方程有Fx_root函数，Fz_roots用于复系数。

案例：解f(x)=x^2-4
系数：
0次项：-4；
1次项：0
2次项：1
IDL> fz_roots([-4,0,1])
(      -2.0000000,      0.00000000)(       2.0000000,      0.00000000)
案例：解p(x)=exp[sin(x^2)+cos(x^2)-1]-1
function func,x
return,exp(sin(x)^2+cos(x)^2-1)-1
end
pro linshi
x=[0.0,-!pi/2,!pi]
root=fx_root(x,'func',/double)
print,root
end
fx_root函数中的x是这么说的：
A 3-element real or complex initial guess vector.
也就是说使用了“先验值”，算法原理写在了Numerical Recipes in C: The Art of Scientific Computing (Second Edition)一书9.5中(中文版见附件)。