求助！！ 数值分析程序（c语言）

西西

<>本人紧急求助，希望各位大虾们伸出你们友情之手，帮我一下。</P>
<>1。牛顿迭代法</P>
<>1。全主消元法</P>
0 Q0 }7 s9 \8 M4 U9 S$ z, S$ o<>3。改进平方根法</P>
( k  U. S+ L) Y, p: h<>4。牛顿向前，向后插值法</P>
<>5。加权最小二乘法</P>
<>6。龙贝格积分法</P>
<>拜托了。。。。</P>

loveyuan

<>1。牛顿迭代法</P><>
9 K6 g: r- U5 P; r1 I3 n  #include "stdio.h"
9 w$ `5 B* d1 o' D$ V0 w  #include "math.h"
  int dnewt(x,eps,js)
+ D' @- Y6 x/ {0 b: U) m  int js;
6 n9 I8 t: c' T: }  double *x,eps;
- B+ l* i$ W, F# Z+ t0 \  { extern void dnewtf();
. J) E4 V! U# |5 D    int k,l;
7 U- ]6 Y8 F/ M% u$ Y    double y[2],d,p,x0,x1;
9 n0 B+ l) R* G% y) A5 S+ D+ q    l=js; k=1; x0=*x;
0 d6 B7 _2 J' f    dnewtf(x0,y);
$ A1 ^& F7 q% ?5 ]$ g1 H, S0 e    d=eps+1.0;
" J, b6 O. ~& w! x    while ((d&gt;=eps)&amp;&amp;(l!=0))
! g/ \2 n; j% `% j$ X+ z      { if (fabs(y[1])+1.0==1.0)
          { printf("err\n"); return(-1);}
        x1=x0-y[0]/y[1];
        dnewtf(x1,y);
        d=fabs(x1-x0); p=fabs(y[0]);
" g& H6 b  n- n0 @% A% G5 `% @        if (p&gt;d) d=p;
        x0=x1; l=l-1;
; w2 p2 ?# L9 ~* ?* ~      }
5 b7 |. D1 o" Z    *x=x1;
    k=js-l;
! f& ~" Z) s; R0 i) R6 e. L    return(k);
  }</P><>全主消元法</P><>#include "stdlib.h"
  #include "stdio.h"
  int acgas(ar,ai,n,br,bi)
* q# o2 ]- f. Q" Q* z( q  int n;
6 x5 ~* Q: T( d' F/ U+ i* p. ]  double ar[],ai[],br[],bi[];
  { int *js,l,k,i,j,is,u,v;
    double p,q,s,d;
! m* p8 u: I" d7 s& k    js=malloc(n*sizeof(int));
- V; c1 H+ S7 S$ Q$ y/ n$ B, K  ~  e# ]    for (k=0;k&lt;=n-2;k++)
' i/ S) C& V' u! }1 k      { d=0.0;
        for (i=k;i&lt;=n-1;i++)
        for (j=k;j&lt;=n-1;j++)
- q, B  U! j5 e0 o  K9 r          { u=i*n+j;
; C/ A/ V; l' @) B9 K            p=ar*ar+ai*ai;
3 e  t- j3 ]5 r9 T  }            if (p&gt;d) {d=p;js[k]=j;is=i;}
1 k9 P+ X3 c1 ?4 u! j% Z          }
        if (d+1.0==1.0)
: C  R' Z1 q2 i& S. w# s+ ]) {          { free(js); printf("err**fail\n");
$ _1 ~9 W/ n( N+ R' \1 M            return(0);
          }
        if (is!=k)
          { for (j=k;j&lt;=n-1;j++)
# ]$ u! T0 _7 z0 P- [7 }/ F              { u=k*n+j; v=is*n+j;
                p=ar; ar=ar[v]; ar[v]=p;
6 P4 }% m$ p8 C$ m) x3 R3 F                p=ai; ai=ai[v]; ai[v]=p;
              }
0 w- _+ ^: U2 ~# b+ H2 o+ m            p=br[k]; br[k]=br[is]; br[is]=p;
            p=bi[k]; bi[k]=bi[is]; bi[is]=p;
( Y- ^9 o4 }" j; Q& B& m          }
        if (js[k]!=k)
) G  G) C- c" [1 _; l          for (i=0;i&lt;=n-1;i++)
            { u=i*n+k; v=i*n+js[k];
              p=ar; ar=ar[v]; ar[v]=p;
) m1 _9 [5 A! z              p=ai; ai=ai[v]; ai[v]=p;
3 K: |, y( e# j/ l8 j            }
        v=k*n+k;
- S9 I( u) P+ |' w4 j  _        for (j=k+1;j&lt;=n-1;j++)
  M3 C6 p7 m( f+ R          { u=k*n+j;
! |6 x% K; \- x            p=ar*ar[v]; q=-ai*ai[v];
! F# u' i- k% _  a            s=(ar[v]-ai[v])*(ar+ai);
! `1 c  a- ?4 Y            ar=(p-q)/d; ai=(s-p-q)/d;
          }
        p=br[k]*ar[v]; q=-bi[k]*ai[v];
, K: g1 Z  Q: \1 q2 h. C0 }( e. O        s=(ar[v]-ai[v])*(br[k]+bi[k]);
        br[k]=(p-q)/d; bi[k]=(s-p-q)/d;
( y1 a$ X2 y. o5 K        for (i=k+1;i&lt;=n-1;i++)
6 Q% t/ K: G; z. a. B  \          { u=i*n+k;
            for (j=k+1;j&lt;=n-1;j++)
              { v=k*n+j; l=i*n+j;
                p=ar*ar[v]; q=ai*ai[v];
                s=(ar+ai)*(ar[v]+ai[v]);
                ar[l]=ar[l]-p+q;
1 b" r! S" {2 t: W( ^& d" ]) T                ai[l]=ai[l]-s+p+q;
; I% n. W" X" ^3 N              }
9 _1 B+ O4 ?6 u; j& x! T, y            p=ar*br[k]; q=ai*bi[k];
            s=(ar+ai)*(br[k]+bi[k]);
            br=br-p+q; bi=bi-s+p+q;
+ C- O" s9 |% H          }
      }
" r" @; J" m3 \( T3 {    u=(n-1)*n+n-1;
1 T# g' E6 M7 o8 e: D    d=ar*ar+ai*ai;
    if (d+1.0==1.0)
3 {* p% x8 ^/ {. {; V+ k      { free(js); printf("err**fail\n");
9 C2 [' o8 M4 v: E) V        return(0);
% R/ R' e7 L9 j3 w- e: F      }
    p=ar*br[n-1]; q=-ai*bi[n-1];
, k4 e5 p( N2 L* _; W  e    s=(ar-ai)*(br[n-1]+bi[n-1]);
6 K/ A+ A/ L% E$ [/ J    br[n-1]=(p-q)/d; bi[n-1]=(s-p-q)/d;
    for (i=n-2;i&gt;=0;i--)
& {1 w3 L- o/ O/ [4 f* W    for (j=i+1;j&lt;=n-1;j++)
      { u=i*n+j;
        p=ar*br[j]; q=ai*bi[j];
! Y; n/ e5 ^% Z  f. |        s=(ar+ai)*(br[j]+bi[j]);
        br=br-p+q;
        bi=bi-s+p+q;
* u; ]' F+ x3 N' p2 u4 X      }
    js[n-1]=n-1;
    for (k=n-1;k&gt;=0;k--)
* @. f" S  n& e2 x* z      if (js[k]!=k)
' o( p; Y- I/ ]- G% [        { p=br[k]; br[k]=br[js[k]]; br[js[k]]=p;
7 J, c9 \! i$ M. O( k          p=bi[k]; bi[k]=bi[js[k]]; bi[js[k]]=p;
        }
    free(js);
& `  p9 ?" m5 U5 x; ~    return(1);
1 Z, q* w' J" x8 I0 q  }</P><>平方根法</P><>#include "math.h"
  #include "stdio.h"
; i% B8 L7 e, @0 ], n1 E& [9 Q: `  int achol(a,n,m,d)
  int n,m;
  double a[],d[];
  { int i,j,k,u,v;
! t) Z$ ], C$ F7 F) d% J( u    if ((a[0]+1.0==1.0)||(a[0]&lt;0.0))
      { printf("fail\n"); return(-2);}
) N- H& Q" \2 b. s* Z/ G% c) a8 B* h    a[0]=sqrt(a[0]);
3 I- o) D9 D" I% ~5 n4 c! J    for (j=1; j&lt;=n-1; j++) a[j]=a[j]/a[0];
    for (i=1; i&lt;=n-1; i++)
2 L) f- ?, }, d& H8 H6 R/ t* H      { u=i*n+i;
        for (j=1; j&lt;=i; j++)
          { v=(j-1)*n+i;
: \& Z, K' O2 G% ^* C/ c4 r3 P            a=a-a[v]*a[v];
          }
        if ((a+1.0==1.0)||(a&lt;0.0))
- n% g) g/ a' y: K4 X          { printf("fail\n"); return(-2);}
$ P, \% M6 ?, @- y+ D0 ^& L5 v# i        a=sqrt(a);
        if (i!=(n-1))
          { for (j=i+1; j&lt;=n-1; j++)
              { v=i*n+j;
                for (k=1; k&lt;=i; k++)
& I! n+ u3 W" o4 r: Y4 [$ u                  a[v]=a[v]-a[(k-1)*n+i]*a[(k-1)*n+j];
                a[v]=a[v]/a;
- F5 e5 W& m/ L5 o- z8 d( H              }
# o; y% K. T6 z. B! Y          }
! ?/ Q5 ?5 f; v! ?( Y9 F- M. A      }
) ]' e5 z7 K, O9 o    for (j=0; j&lt;=m-1; j++)
      { d[j]=d[j]/a[0];
- i& y+ v2 ?4 b5 ]% e        for (i=1; i&lt;=n-1; i++)
! `# P( \4 T2 L( \" W7 d          { u=i*n+i; v=i*m+j;
& n8 |+ p4 A& w& a            for (k=1; k&lt;=i; k++)
5 x) N! z: w3 [+ u4 k/ @7 }3 B3 z              d[v]=d[v]-a[(k-1)*n+i]*d[(k-1)*m+j];
1 P1 o5 H$ S+ d) {            d[v]=d[v]/a;
% n* U- a# J) q6 w7 M- _          }
      }
    for (j=0; j&lt;=m-1; j++)
( P- A4 ^0 R% {! F6 A      { u=(n-1)*m+j;
2 D8 F! E( g. S* V, }        d=d/a[n*n-1];
        for (k=n-1; k&gt;=1; k--)
$ `! k5 n: Q! l0 b. G& e          { u=(k-1)*m+j;
            for (i=k; i&lt;=n-1; i++)
              { v=(k-1)*n+i;
/ g7 r$ N3 |) \' b                d=d-a[v]*d[i*m+j];
& s  K! |% `; M! ]% t              }
            v=(k-1)*n+k-1;
5 `, }+ w* N* o            d=d/a[v];
          }
8 r- r6 x+ f- H: L$ y" s$ E3 Y5 j      }
; S, g6 [8 n" T( G6 Y    return(2);
, c* I" b& ]5 q" R$ j  }</P><>牛顿向前，向后插值法</P><>double eelgr(x0,h,n,y,t)
  int n;
/ m! P5 m  x& k) g' s; V% `5 Q  double x0,h,t,y[];
; L# V7 j8 J/ R2 L" e  ?; q, u  { int i,j,k,m;
7 X. x& i( o3 _% m1 M    double z,s,xi,xj;
    float p,q;
    z=0.0;
    if (n&lt;1) return(z);
  a7 I$ o( J) ]* l( o8 Z0 I    if (n==1) { z=y[0]; return(z);}
    if (n==2)
5 T; S% [2 I8 \; m      { z=(y[1]*(t-x0)-y[0]*(t-x0-h))/h;
        return(z);
      }
    if (t&gt;x0)
      { p=(t-x0)/h; i=(int)p; q=(float)i;
        if (p&gt;q) i=i+1;
      }
    else i=0;
    k=i-4;
    if (k&lt;0) k=0;
* s3 _4 }: ?' o' ^    m=i+3;
; |7 l! C$ V: R( x% O9 t2 t    if (m&gt;n-1) m=n-1;
    for (i=k;i&lt;=m;i++)
      { s=1.0; xi=x0+i*h;
        for (j=k; j&lt;=m; j++)
) d! y+ [2 x1 F          if (j!=i)
; {/ X* Y. g! o) v; t! s            { xj=x0+j*h;
# t) A- t$ ?( e              s=s*(t-xj)/(xi-xj);
4 m' M# R7 @, z: z3 G            }
  K( n( e. K3 d( w1 @$ Q  P; }        z=z+s*y;
% n3 j! u' y/ ~* I) s. e3 p      }
; q& Q+ j+ E5 Q5 i, }    return(z);
  }
2 p  Q8 l9 b0 {$ ?向前，向后是一样的思想！！</P><>加权最小二乘法</P><>#include "math.h"
  void hpir1(x,y,n,a,m,dt)
  int n,m;
  double x[],y[],a[],dt[];
  { int i,j,k;
- Z0 C; Y# X" J4 [2 u    double z,p,c,g,q,d1,d2,s[20],t[20],b[20];
    for (i=0; i&lt;=m-1; i++) a=0.0;
    if (m&gt;n) m=n;
    if (m&gt;20) m=20;
4 T" j3 B7 d( d  ^+ c    z=0.0;
    for (i=0; i&lt;=n-1; i++) z=z+x/(1.0*n);
    b[0]=1.0; d1=1.0*n; p=0.0; c=0.0;
    for (i=0; i&lt;=n-1; i++)
      { p=p+(x-z); c=c+y;}
0 V, ]; z( e& g  \    c=c/d1; p=p/d1;
    a[0]=c*b[0];
/ M8 J5 p+ }8 K7 S% [    if (m&gt;1)
2 M* L% f  L; T. w. {, Y      { t[1]=1.0; t[0]=-p;
; U( I5 B& h7 N4 a2 h0 B        d2=0.0; c=0.0; g=0.0;
, c  I6 |# R% ~3 h+ ?2 n        for (i=0; i&lt;=n-1; i++)
          { q=x-z-p; d2=d2+q*q;
5 ~5 y+ T. E1 T/ b            c=c+y*q;
; q7 w7 k* X, ~8 e1 o' I( Q+ O) c            g=g+(x-z)*q*q;
4 B; d- G- S7 E0 a' j8 H: r& s          }
6 G! V: p# _( {" D7 t) d- `' p7 X        c=c/d2; p=g/d2; q=d2/d1;
        d1=d2;
: d* n2 t: s2 h        a[1]=c*t[1]; a[0]=c*t[0]+a[0];
      }
; b3 U- k8 |0 v& `: U4 v" {    for (j=2; j&lt;=m-1; j++)
- E: Y' }4 @' x5 y' U      { s[j]=t[j-1];
* U3 D/ k- }( {" K        s[j-1]=-p*t[j-1]+t[j-2];
& U' x2 S( _; b        if (j&gt;=3)
5 e5 Q. E- W" v6 x7 b! e          for (k=j-2; k&gt;=1; k--)
" D3 C; ~4 Q" C6 f" r/ l% `            s[k]=-p*t[k]+t[k-1]-q*b[k];
& c5 ~1 n! u5 j1 ~3 s        s[0]=-p*t[0]-q*b[0];
/ t6 u, o; ]3 l        d2=0.0; c=0.0; g=0.0;
        for (i=0; i&lt;=n-1; i++)
          { q=s[j];
            for (k=j-1; k&gt;=0; k--)
8 n, K$ R) s3 j! q% c7 T              q=q*(x-z)+s[k];
1 s2 g# a* W: X: t6 }            d2=d2+q*q; c=c+y*q;
            g=g+(x-z)*q*q;
          }
5 K# E9 V! D( l! j8 ?" K3 C9 A, n        c=c/d2; p=g/d2; q=d2/d1;
        d1=d2;
        a[j]=c*s[j]; t[j]=s[j];
        for (k=j-1; k&gt;=0; k--)
; y9 |$ |8 J4 J. j& F, d, j$ J          { a[k]=c*s[k]+a[k];
( v6 N$ P+ b! z( K, [            b[k]=t[k]; t[k]=s[k];
$ D' A/ X) |! j% X          }
, v( V/ ]' o1 N, z  c$ W      }
, u' h  G; H2 C: ?    dt[0]=0.0; dt[1]=0.0; dt[2]=0.0;
    for (i=0; i&lt;=n-1; i++)
      { q=a[m-1];
( E& C1 c2 B( [# A8 _+ Y) W        for (k=m-2; k&gt;=0; k--)
          q=a[k]+q*(x-z);
        p=q-y;
' L5 [# M3 y  [/ U7 _% w        if (fabs(p)&gt;dt[2]) dt[2]=fabs(p);
* d& T- V! B8 A8 |' J2 M        dt[0]=dt[0]+p*p;
- U5 j3 j( g; y1 ?        dt[1]=dt[1]+fabs(p);
; U1 G  B1 e! C4 e( J/ V7 }; g! a      }
    return;
  }</P><>龙贝格积分法</P><>#include "math.h"
  double fromb(a,b,eps)
8 v  n, w1 b! h# J  double a,b,eps;
+ D/ g+ ~4 O) I* m  { extern double frombf();
" J9 U! }! X3 u  _( z( {3 o7 f7 u1 _    int m,n,i,k;
1 t$ k, v9 \8 `* h0 X! V    double y[10],h,ep,p,x,s,q;
! W4 T5 }; {  X6 T' l% w    h=b-a;
0 v2 M5 o! }2 ^1 V3 e1 s    y[0]=h*(frombf(a)+frombf(b))/2.0;
6 E+ d0 i6 Y! x8 K) L  J3 L    m=1; n=1; ep=eps+1.0;
    while ((ep&gt;=eps)&amp;&amp;(m&lt;=9))
      { p=0.0;
        for (i=0;i&lt;=n-1;i++)
; |6 v4 L- s) t, i. d4 W# J% l1 ~          { x=a+(i+0.5)*h;
            p=p+frombf(x);
          }
4 [6 H) {0 t" [6 k; p        p=(y[0]+h*p)/2.0;
$ S6 k1 w) I2 J; N) R        s=1.0;
        for (k=1;k&lt;=m;k++)
          { s=4.0*s;
            q=(s*p-y[k-1])/(s-1.0);
            y[k-1]=p; p=q;
/ T' x' `  R  U          }
        ep=fabs(q-y[m-1]);
        m=m+1; y[m-1]=q; n=n+n; h=h/2.0;
      }
    return(q);
9 F1 ~. Y" ]4 p: ]  w9 F" |& C) h& v& k  }</P><>呵呵 希望对你有用！！</P>

西西

<>谢谢你哦，呵呵</P>

陈了

哇哇哇哇哇哇。。。好复杂呀！看不懂