> !!!本程序适用于求解形如f(x)=1/2*x'Ax+bx+c二次函数的稳定点;
,allocatable::x,gradt,dir,b ,s,y,gradt1,x1
,allocatable::hessin ,B1 ,G,G1
> do i=1,n, t& ~+ r2 ^; K5 t
> !!!子程序,返回函数值 5 R& P: i' e$ C3 r: [; }8 S
,intent(in)::x,b* w f) ]/ d% [" K2 O6 x
,intent(in)::A! m3 s; _$ j" ^3 L6 z3 j2 B4 ~
,intent(in)::p
,intent(out)::G
=p(i)*p
,intent(in)::x,d
,intent(in)::b9 M$ X" _% S! H4 R6 o: t. v
,intent(in)::A
>
>本程序用Fortran 90编写,在Virual Fortran 5上编译通过,本程序由沙沙提供!</P>| 欢迎光临 数学建模社区-数学中国 (http://www.madio.net/) | Powered by Discuz! X2.5 |