> !!!本程序适用于求解形如f(x)=1/2*x'Ax+bx+c二次函数的稳定点;, Q4 }* W1 L/ N& R4 J
,allocatable::x,x1,gradtf,gradts,dirf,dirs,b
,allocatable::hessin
> !!!子程序,返回函数值
,intent(in)::x,b
,intent(in)::A' u u) y9 t' b, z9 ^) b5 j
> !!!精确线搜索0.618法子程序,返回迭代步长% ]' ^$ E$ c; V3 p5 P2 R* A% Q
,intent(in)::x,d
,intent(in)::b
,intent(in)::A
>本程序由Fortran90编写,在Vistual Fortran 5上调试通过!希望大家批评指正!</P>
>谢谢</P>| 欢迎光临 数学建模社区-数学中国 (http://www.madio.net/) | Powered by Discuz! X2.5 |