> !!!本程序适用于求解形如f(x)=1/2*x'Ax+bx+c二次函数的稳定点;
,allocatable::x,gradt,dir,b ,s,y,gradt1,x1
,allocatable::hessin ,B1 ,G,G1! x$ p+ L9 H" T8 B) I1 X' t1 ]
> do i=1,n
> !!!子程序,返回函数值 V4 o* s) ]' C- P% X* @& }
,intent(in)::x,b. S8 y# L4 N% @5 [) {
,intent(in)::A$ r& O6 a4 ^; O# c
,intent(in)::p( z* }/ o+ }- T5 r1 I, ^, j, {. U
,intent(out)::G% Y9 t9 G* }# x. N5 K
=p(i)*p5 ]* J5 r( l" L& G i
,intent(in)::x,d
,intent(in)::b, |7 U$ H; j+ ^ L% r- E7 N- W$ m" O4 X/ v
,intent(in)::A) r E! {8 v( h# X
> 4 g* J* m( C v+ |" l
>本程序用Fortran 90编写,在Virual Fortran 5上编译通过,本程序由沙沙提供!</P>| 欢迎光临 数学建模社区-数学中国 (http://www.madio.net/) | Powered by Discuz! X2.5 |