> !!!本程序适用于求解形如f(x)=1/2*x'Ax+bx+c二次函数的稳定点;; ]' D Z: X/ i% N. t5 c6 D1 H
,allocatable::x,gradt,dir,b ,s,y,p ,gradt1,x1* l) m9 i o9 s, H( o% V. u
,allocatable::hessin ,H ,G ,U
> do i=1,n
> !!!子程序,返回函数值
,intent(in)::x,b) u: V9 n* D$ @; n) x5 d
,intent(in)::A1 d- ^* r& [, w" r3 i" h+ n# f
,intent(in)::p3 `0 Y2 b j1 G ^3 ]- [3 I/ t8 [8 Q. m
,intent(out)::G( }/ L, J+ p8 t! |. d
,intent(in)::x,d
,intent(in)::b
,intent(in)::A. t6 D2 D) x) M* @" ?! z$ w C/ R
>!!!本程序适用于求解形如f(x)=1/2*x'Ax+bx+c二次函数的稳定点;2 H8 l# Y }" o2 [
,allocatable::x,gradt,dir,b ,s,y,p ,gradt1,x1
,allocatable::hessin ,H ,G ,U
> do i=1,n. G8 H+ u9 q- }9 t2 w
> !!!子程序,返回函数值 / Q4 |7 M/ e, M# a2 ]4 q( N6 _
,intent(in)::x,b
,intent(in)::A$ |( N7 f7 v4 S/ h, T
,intent(in)::p
,intent(out)::G0 N0 f" x( k6 H1 W6 _
,intent(in)::x,d
,intent(in)::b7 W9 ^; a( A* \& i& `. {
,intent(in)::A
>本程序由Fortran 90编写,在Visual Fortran 5上编译通过,本程序由沙沙提供!| 欢迎光临 数学建模社区-数学中国 (http://www.madio.net/) | Powered by Discuz! X2.5 |