非线性规划程序

dx999

#include<iostream.h> 0 }3 l' r7 I3 }" _2 q#include<math.h> double fx(int n,double p[],double x) { double z=0,y=1; 0 R8 @' I+ I7 i5 y* }7 D int i=0,j=0; for(;i<n+1;i++) { for(;j<i;j++) ; _, c1 L6 l$ o$ Z! F9 w y*=x; 3 M& R( z# I z5 m. v, }; E z+=*(p+i)*y; } : s* ~# v ^8 X/ f& K5 ~* | return z; }

* o6 i# i3 p3 i- o. v

int fi(int k) 6 Z2 t- E1 G/ ~" y{ if(k==0) , o% Z* w# T) E1 M return 1; else { if(k==1) return 1; 5 o6 p! N8 C( S' K# ~% Y9 e else { int i=0,i1=1,i2=1,w=0; 1 N5 P# n) x. L* ]5 e0 J for(;i<k-1;i++) ) W; u9 s/ Q4 o4 @) X { w=i1+i2; i1=i2; i2=w; : {$ w/ ]4 K, ?( T1 q8 Z/ Y7 w' x& M% y } 6 \" a/ X+ ]3 y% b. k return w; # g3 l0 H; @7 V h } 6 {% z+ j$ d) H# e+ e7 Q4 `} }

# s. F6 |1 V& k. ^. ^( O) \3 f

int r(double w,double a,double b) { int i=0; ( d4 e, H1 O e' L for(;fi(i)<=(b-a)/w;i++) {} return i; } : ~" c' @. D7 S7 }7 O, _0 @void main() , H/ r8 h0 o& b6 S6 l{ cout<<"输入函数的最高次幂"; int n; cin>>n; double *p=new double[n+1]; double a,b,g=0,j=0,c=0; ( g' y3 h$ k$ Q" h& b& s. |- X cout<<"依升降幂输入x项的系数"<<endl; T) b4 D- {& }' F$ b# E* \ for(int i=0;i<n+1;i++) cin>>*(p+i); " ?# T. j7 K# A) W# i cout<<"输入误差值"<<endl; ! C; d5 i5 t. g$ K* A1 O double w; P$ d* o$ a& }: B6 N2 g; A) o6 O7 | cin>>w; cout<<"输入初始区间"<<endl; Y: K' `& `( Y4 S9 m2 ? cin>>a>>b; cout<<"输入辨别常数"<<endl; cin>>c; ! `3 ^! x6 Y' a. r0 x9 P n=r(w,a,b); i=0; 2 N! s& [/ h l& A& E6 f cout<<fi(n-i-2)<<endl; cout<<a<<" "<<b; 7 Y2 n1 e1 E$ _$ Q$ k/ v g=a+fi(n-2-i)*(b-a)/fi(n-i); $ ]8 U/ P- H3 f, Z" m& w j=a+fi(n-1-i)*(b-a)/fi(n-i); cout<<g<<" "<<j; for(i=1;i<n-1;i++) * u! @8 P4 P6 v6 }/ J { 9 t0 f) Q# S4 p8 v/ l 9 W5 u/ Z. f7 w% ^ F9 B8 _ if(fx(n,p,g)>fx(n,p,j)) - F, k$ G1 V: B9 K6 r { - H* @3 H: ^- n+ y% N a=g; ( Y. ~$ \' ^& l4 z. Y g=j; ( }9 P4 s5 V# ~/ ]; N- k j=a+fi(n-1-i)*(b-a)/fi(n-i); } - }( ~( S# _/ E( S6 B m9 o else " p3 s s1 J/ ^* W3 i9 b0 e9 P { b=j; j=g; g=a+fi(n-2-i)*(b-a)/fi(n-i); } & B4 A7 R( n$ t/ ?5 {4 } } - C9 H! M P& p. \+ ` j=g+c; if(fx(n,p,g)>fx(n,p,j)) . B4 c1 w1 r: X2 T6 n a=g; + _* U1 W1 y6 u/ [ else b=g; 5 ?- o5 q$ k' }" {5 W5 R cout<<"函数值在 ("<<g<<","<<j<<")之间；"; - I3 D8 ?; g! u( f }

无语的昭昭

c啊，，，，历害，，我一般用lingo的说

weigao1234

lsyzzhl

城榜公_DNO

沙发，板凳都没偶的份了。
! L8 i9 Y! k* y9 e) Y

9 j# l3 U# P) u+ y3 e& }7 ]
" o7 R' x& ?( h: m0 j$ s, j

  o  \  x: U5 ~  s- _2 J

7 X- p! u) g" w
% M  a0 `1 d; a
5 C- j7 s. O" L8 i
1 N% c9 {" m9 T3 L3 C1 \9 x
- @/ V/ w# {1 C) N" t+ J

淘宝专享地址：zhfqw.taobao.com/

6 L5 ?- B2 k; B8 S如何将宝贝卖向全世界           索马里论坛**598专卖

2 H/ C( N5 ~0 E2 z7 n. Q0 x6 Q如何将宝贝卖向全世界           索马里论坛**598专卖

shuxuezaozhuang

有些乱！！

水木年华zzu

楼主要是能把格式整理下就更好了

515genius

慢慢看好了

家家开水

头都昏了···

999mmg

hehe&nbsp;&nbsp; 绝对是好东西啊