> !!!本程序适用于求解形如f(x)=1/2*x'Ax+bx+c二次函数的稳定点;. s1 Y* g# H$ w* u9 W4 X
,allocatable::x,gradt,dir,b ,s,y,p ,gradt1,x1
,allocatable::hessin ,H ,G ,U
> do i=1,n
> !!!子程序,返回函数值
,intent(in)::x,b
,intent(in)::A. `7 i+ T( F2 {8 [* i, W
,intent(in)::p! C# y: s- E. G7 o: m2 @3 w
,intent(out)::G+ j" c4 y6 M/ q. `8 d& V! r( d
,intent(in)::x,d
,intent(in)::b
,intent(in)::A
>!!!本程序适用于求解形如f(x)=1/2*x'Ax+bx+c二次函数的稳定点;
,allocatable::x,gradt,dir,b ,s,y,p ,gradt1,x1
,allocatable::hessin ,H ,G ,U6 Z4 n/ r; q/ l& g2 F
> do i=1,n
> !!!子程序,返回函数值 ) y8 A7 h. y2 s
,intent(in)::x,b) L4 ? ^: T$ D: s8 m
,intent(in)::A1 t7 w0 u; t" d+ [+ u% d+ D
,intent(in)::p4 T8 x! [2 T% o$ F8 R- w" G+ T( U
,intent(out)::G
,intent(in)::x,d
,intent(in)::b
,intent(in)::A
>本程序由Fortran 90编写,在Visual Fortran 5上编译通过,本程序由沙沙提供!| 欢迎光临 数学建模社区-数学中国 (http://www.madio.net/) | Powered by Discuz! X2.5 |