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

西西

<>本人紧急求助，希望各位大虾们伸出你们友情之手，帮我一下。</P>
<>1。牛顿迭代法</P>
<>1。全主消元法</P>
<>3。改进平方根法</P>
1 C0 m& S0 T6 H# B( K4 r! f- H" I! D<>4。牛顿向前，向后插值法</P>
+ p0 E) K  G0 N: P3 k<>5。加权最小二乘法</P>
<>6。龙贝格积分法</P>
6 b( [  T5 \& v<>拜托了。。。。</P>

loveyuan

<>1。牛顿迭代法</P><>
  #include "stdio.h"
  #include "math.h"
  int dnewt(x,eps,js)
  int js;
. _- ]  ?* O) c" n- J& j; G# J  double *x,eps;
2 V5 d. D" w" T% X  { extern void dnewtf();
7 s  d' I2 {- p    int k,l;
, e/ ~" r  D3 S3 [5 \- g! M    double y[2],d,p,x0,x1;
- i( j2 \1 o; [" U* i) x    l=js; k=1; x0=*x;
# v- j5 G& z, e( h( _+ _9 D% n* v, ]    dnewtf(x0,y);
* U3 A/ `- }9 E  V. g6 }8 G    d=eps+1.0;
    while ((d&gt;=eps)&amp;&amp;(l!=0))
      { if (fabs(y[1])+1.0==1.0)
          { printf("err\n"); return(-1);}
0 |9 n: V6 A: X( R' m        x1=x0-y[0]/y[1];
& \, _0 N" \  i# Q3 X* ?, r        dnewtf(x1,y);
        d=fabs(x1-x0); p=fabs(y[0]);
        if (p&gt;d) d=p;
5 c( K( I6 A" W" C( M        x0=x1; l=l-1;
3 n& n0 X8 T0 \# r6 Q& c      }
2 b5 Q& k2 I0 A# R    *x=x1;
/ q% k$ u" T% Y( Z$ [    k=js-l;
    return(k);
  }</P><>全主消元法</P><>#include "stdlib.h"
8 i: i0 U/ q2 k7 P7 u  #include "stdio.h"
2 c) S  W1 S2 U( h- V. U( q  int acgas(ar,ai,n,br,bi)
  int n;
& r( z+ ]: x7 }1 E: e8 i  double ar[],ai[],br[],bi[];
  { int *js,l,k,i,j,is,u,v;
. w- X. u' E. J! }" B6 J* l    double p,q,s,d;
. Q) h& o# I% @8 l' n. _: O    js=malloc(n*sizeof(int));
1 Z( S9 f  x- p) F) R9 A. d) C    for (k=0;k&lt;=n-2;k++)
' L- |4 g& G/ ~* l      { d=0.0;
. E% V) K5 F% _, R/ Y; y        for (i=k;i&lt;=n-1;i++)
        for (j=k;j&lt;=n-1;j++)
          { u=i*n+j;
. X$ l' U/ C0 T1 y1 D            p=ar*ar+ai*ai;
            if (p&gt;d) {d=p;js[k]=j;is=i;}
          }
        if (d+1.0==1.0)
& X5 k) M0 c3 H9 s7 B( X/ z5 P          { free(js); printf("err**fail\n");
            return(0);
          }
& T4 B  X8 L! S. R! ?; G# E( `        if (is!=k)
" S2 m( [/ n3 c& A9 Z, G6 q          { for (j=k;j&lt;=n-1;j++)
1 ?  c! Z; b9 F( F              { u=k*n+j; v=is*n+j;
/ S) b; U* c! K6 g                p=ar; ar=ar[v]; ar[v]=p;
                p=ai; ai=ai[v]; ai[v]=p;
              }
% n7 S- j6 r" |$ \% |6 Z! ]2 ]% t            p=br[k]; br[k]=br[is]; br[is]=p;
( y# ~; d" D( a            p=bi[k]; bi[k]=bi[is]; bi[is]=p;
" m: \6 F& q7 s, @8 L. X( Q$ Z0 F          }
* ^7 d; w0 A5 l' E7 K        if (js[k]!=k)
4 f9 c7 A. \" }2 E( r2 Y8 S          for (i=0;i&lt;=n-1;i++)
            { u=i*n+k; v=i*n+js[k];
2 M) \  O5 b5 l3 p6 y' w$ W2 p; f              p=ar; ar=ar[v]; ar[v]=p;
              p=ai; ai=ai[v]; ai[v]=p;
- ~+ e! K+ w% s) ~, G6 n: O            }
        v=k*n+k;
        for (j=k+1;j&lt;=n-1;j++)
          { u=k*n+j;
            p=ar*ar[v]; q=-ai*ai[v];
& [; ?- ^1 V4 w  h) t            s=(ar[v]-ai[v])*(ar+ai);
9 i1 u* s% ^- [3 p" y            ar=(p-q)/d; ai=(s-p-q)/d;
& b3 t9 S3 G: J  j' Q          }
, m$ T) O" p# M4 d- a8 _        p=br[k]*ar[v]; q=-bi[k]*ai[v];
# g) _+ U: \* N' \        s=(ar[v]-ai[v])*(br[k]+bi[k]);
. t! e5 \9 N, P+ F& {* U        br[k]=(p-q)/d; bi[k]=(s-p-q)/d;
        for (i=k+1;i&lt;=n-1;i++)
% S8 C4 G+ [/ ?          { u=i*n+k;
            for (j=k+1;j&lt;=n-1;j++)
              { v=k*n+j; l=i*n+j;
+ T; ]/ {" Y) x$ {" r                p=ar*ar[v]; q=ai*ai[v];
$ k; x  m, K0 K7 x: ^! T                s=(ar+ai)*(ar[v]+ai[v]);
  F5 b& y4 l) t0 e7 j                ar[l]=ar[l]-p+q;
) b# J. H6 q; l1 p6 v; T                ai[l]=ai[l]-s+p+q;
              }
1 k. E* s' r. m! q& K% l* E            p=ar*br[k]; q=ai*bi[k];
" Y/ B  f  q9 ~$ O0 d" T6 {/ {0 K            s=(ar+ai)*(br[k]+bi[k]);
' F2 d; _! m- q+ x2 v8 z: K! B            br=br-p+q; bi=bi-s+p+q;
          }
  h- l2 j" P8 x7 R& Z7 B; v      }
* \% x2 l9 h5 ~  G  q( ?    u=(n-1)*n+n-1;
- |$ b0 f% q+ ?* P& n  J" p+ u    d=ar*ar+ai*ai;
    if (d+1.0==1.0)
) p# C. w$ k4 m; {4 a# a      { free(js); printf("err**fail\n");
        return(0);
      }
    p=ar*br[n-1]; q=-ai*bi[n-1];
    s=(ar-ai)*(br[n-1]+bi[n-1]);
7 P, A: J: O: B6 s  I- }+ c+ i    br[n-1]=(p-q)/d; bi[n-1]=(s-p-q)/d;
    for (i=n-2;i&gt;=0;i--)
' ?: a2 J1 ~* E+ c    for (j=i+1;j&lt;=n-1;j++)
      { u=i*n+j;
6 x+ t, Z8 C! S" w5 F        p=ar*br[j]; q=ai*bi[j];
3 t8 x: l% s3 F8 q; e        s=(ar+ai)*(br[j]+bi[j]);
        br=br-p+q;
$ |& T: R5 {) r8 a  f8 s        bi=bi-s+p+q;
' g( Y" S% h2 G6 v, W      }
    js[n-1]=n-1;
    for (k=n-1;k&gt;=0;k--)
      if (js[k]!=k)
* K' l1 l1 O6 Q/ Z+ R" j+ S        { p=br[k]; br[k]=br[js[k]]; br[js[k]]=p;
  A8 {0 x& B$ u/ x: ?, T+ |          p=bi[k]; bi[k]=bi[js[k]]; bi[js[k]]=p;
% i$ z7 K1 Z# h        }
: z2 A. m7 n9 w1 e* T( f: [  K# I+ n    free(js);
4 R7 v+ W, k( ]) z+ ?+ h    return(1);
  }</P><>平方根法</P><>#include "math.h"
  #include "stdio.h"
5 [  z. ~) G( ]2 |9 a! \& u  int achol(a,n,m,d)
  int n,m;
  double a[],d[];
  { int i,j,k,u,v;
    if ((a[0]+1.0==1.0)||(a[0]&lt;0.0))
8 }( G0 e% _; j7 e1 `, s) E      { printf("fail\n"); return(-2);}
. @. m; G6 A$ E    a[0]=sqrt(a[0]);
    for (j=1; j&lt;=n-1; j++) a[j]=a[j]/a[0];
9 M) O" v0 M& d    for (i=1; i&lt;=n-1; i++)
      { u=i*n+i;
        for (j=1; j&lt;=i; j++)
, l/ n! ~: L8 e* S' [          { v=(j-1)*n+i;
6 `9 P0 \* D% H9 W* W7 {/ F            a=a-a[v]*a[v];
          }
        if ((a+1.0==1.0)||(a&lt;0.0))
! l0 G" {8 e$ H          { printf("fail\n"); return(-2);}
( w6 O& I, K1 z, h        a=sqrt(a);
        if (i!=(n-1))
1 `  z0 c* x  X3 P          { for (j=i+1; j&lt;=n-1; j++)
/ [' f( w4 U" `4 h. W  J4 c              { v=i*n+j;
9 {1 M/ r5 O  a" k5 k                for (k=1; k&lt;=i; k++)
6 h+ e5 P3 y" {" ]% w                  a[v]=a[v]-a[(k-1)*n+i]*a[(k-1)*n+j];
                a[v]=a[v]/a;
              }
          }
5 b; w# }; P8 b; h" h% F& y" m      }
0 \; n" K, M  w$ `2 V    for (j=0; j&lt;=m-1; j++)
4 h  f/ Y& M7 R" C! d9 d      { d[j]=d[j]/a[0];
% a: P5 x) `) B) [( J        for (i=1; i&lt;=n-1; i++)
" p( _# }9 K6 h3 ~          { u=i*n+i; v=i*m+j;
. \( b6 j! Q- c& k  f            for (k=1; k&lt;=i; k++)
              d[v]=d[v]-a[(k-1)*n+i]*d[(k-1)*m+j];
0 @) y  A) ]3 l8 h+ U            d[v]=d[v]/a;
          }
      }
$ J; I% @- G$ S" G8 ?4 r    for (j=0; j&lt;=m-1; j++)
3 M& G: }7 q! ]% z      { u=(n-1)*m+j;
        d=d/a[n*n-1];
        for (k=n-1; k&gt;=1; k--)
          { u=(k-1)*m+j;
/ {! ]' u" R) Y5 r; d0 z            for (i=k; i&lt;=n-1; i++)
4 |; [- ]% o$ @; Q& k2 i              { v=(k-1)*n+i;
                d=d-a[v]*d[i*m+j];
1 f) _! m! |9 X7 ?              }
            v=(k-1)*n+k-1;
  X6 y; D, z- b5 c& l: f            d=d/a[v];
          }
* G3 A8 l& Y2 ]" n: O2 @" n      }
  ]$ ]5 R5 Z3 J/ k    return(2);
+ C" ^) r4 q3 t" [" C6 A+ B( q2 f& e  }</P><>牛顿向前，向后插值法</P><>double eelgr(x0,h,n,y,t)
5 [' c  k; |5 b8 d$ M1 j& D6 `1 b  int n;
6 q1 g' z" O7 c  double x0,h,t,y[];
  { int i,j,k,m;
    double z,s,xi,xj;
0 i- I4 {* G9 m% g$ \& |    float p,q;
    z=0.0;
    if (n&lt;1) return(z);
4 R; c# |) E. m3 X' x0 Y- ~7 F    if (n==1) { z=y[0]; return(z);}
8 y3 e4 O) X* n! d+ r% s* B    if (n==2)
$ U& ~, N3 f5 f/ M      { z=(y[1]*(t-x0)-y[0]*(t-x0-h))/h;
; L. z4 D* R. s7 _8 M6 v/ U0 o        return(z);
      }
# r! x7 `2 a. P5 ?8 S% K    if (t&gt;x0)
      { p=(t-x0)/h; i=(int)p; q=(float)i;
& a, n4 f6 o% ^, q) y( A        if (p&gt;q) i=i+1;
& M5 @& X, C$ m. ]/ g6 l, K      }
1 Y" _4 i' v% y6 n) ^% m    else i=0;
    k=i-4;
    if (k&lt;0) k=0;
0 ^( w/ C; n9 K    m=i+3;
: F$ s% [/ D9 L8 _9 ~    if (m&gt;n-1) m=n-1;
1 {8 X8 }/ {" E6 u/ i    for (i=k;i&lt;=m;i++)
      { s=1.0; xi=x0+i*h;
        for (j=k; j&lt;=m; j++)
          if (j!=i)
            { xj=x0+j*h;
- U3 B2 @: d( T& H3 H& f              s=s*(t-xj)/(xi-xj);
            }
        z=z+s*y;
      }
# |' @  t' s6 i- M9 m    return(z);
8 b2 e5 g5 l, S1 E  }
向前，向后是一样的思想！！</P><>加权最小二乘法</P><>#include "math.h"
- M5 ^8 `/ q" z; _) C  void hpir1(x,y,n,a,m,dt)
  `$ e; l/ F* C: u2 ]  int n,m;
  double x[],y[],a[],dt[];
  { int i,j,k;
2 H7 S# ~. v4 p7 k, q) y    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;
) ^0 Y! D: }* a8 w9 M6 Q    if (m&gt;n) m=n;
    if (m&gt;20) m=20;
8 V* d3 y- |3 u) q& `" U    z=0.0;
: Q; b5 U7 f2 x# o  @8 r, k    for (i=0; i&lt;=n-1; i++) z=z+x/(1.0*n);
0 H6 k) q3 g6 T' x: \    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;}
    c=c/d1; p=p/d1;
3 J9 w8 J  [( k1 ]0 d    a[0]=c*b[0];
    if (m&gt;1)
- |0 m( Y. c( y- k- |5 a1 S7 o1 i      { t[1]=1.0; t[0]=-p;
0 n! o& K5 q- T5 j8 h4 ^0 D        d2=0.0; c=0.0; g=0.0;
        for (i=0; i&lt;=n-1; i++)
          { q=x-z-p; d2=d2+q*q;
            c=c+y*q;
' v3 {, l. I3 x# ?, y* b- p" M4 M            g=g+(x-z)*q*q;
          }
! z+ i+ Q1 t0 N; I        c=c/d2; p=g/d2; q=d2/d1;
2 Y' Z0 h, j1 X) |9 I        d1=d2;
$ N* p9 \" w' K! H        a[1]=c*t[1]; a[0]=c*t[0]+a[0];
8 j+ K) b& P0 B7 }% c4 b      }
! e( J0 P3 B1 m+ t  F- ]    for (j=2; j&lt;=m-1; j++)
      { s[j]=t[j-1];
' v" L) m3 ?7 Q, n/ g2 `0 B        s[j-1]=-p*t[j-1]+t[j-2];
        if (j&gt;=3)
4 m- m$ B2 h( P8 B          for (k=j-2; k&gt;=1; k--)
            s[k]=-p*t[k]+t[k-1]-q*b[k];
        s[0]=-p*t[0]-q*b[0];
        d2=0.0; c=0.0; g=0.0;
3 ^6 p* O. {. \0 P3 ?; I        for (i=0; i&lt;=n-1; i++)
          { q=s[j];
            for (k=j-1; k&gt;=0; k--)
              q=q*(x-z)+s[k];
( W6 O: @0 i; N! r            d2=d2+q*q; c=c+y*q;
            g=g+(x-z)*q*q;
9 _/ c0 B! I4 X+ Y8 H8 Z  j$ F9 Z9 @          }
        c=c/d2; p=g/d2; q=d2/d1;
        d1=d2;
8 C9 _9 G+ r+ T# `& i4 y; n0 ]5 N+ u        a[j]=c*s[j]; t[j]=s[j];
        for (k=j-1; k&gt;=0; k--)
+ f5 U- g+ M" ]$ L          { a[k]=c*s[k]+a[k];
            b[k]=t[k]; t[k]=s[k];
! ?8 k6 S- q3 j6 a( L0 B- Q) i, Z          }
      }
    dt[0]=0.0; dt[1]=0.0; dt[2]=0.0;
+ N4 _7 j, Y# I% Z8 ]" f$ z/ p( U    for (i=0; i&lt;=n-1; i++)
  [) j5 Z( U. X4 t: o: s      { q=a[m-1];
        for (k=m-2; k&gt;=0; k--)
          q=a[k]+q*(x-z);
9 l4 |2 l# l6 _8 b8 f0 t        p=q-y;
4 c' b! m; ^5 B  L: F        if (fabs(p)&gt;dt[2]) dt[2]=fabs(p);
        dt[0]=dt[0]+p*p;
        dt[1]=dt[1]+fabs(p);
# [# ]  [# T+ K4 d/ m8 D      }
# H: T3 ~4 c1 i2 b# q3 L6 ?    return;
  }</P><>龙贝格积分法</P><>#include "math.h"
  double fromb(a,b,eps)
' ~" {; o$ x8 n! L0 E  double a,b,eps;
+ D9 ^8 D% t) J1 o& j6 T, ^% e  { extern double frombf();
    int m,n,i,k;
$ r+ C+ m& H" M* Z    double y[10],h,ep,p,x,s,q;
+ P' _! O- i7 ]. z# R    h=b-a;
    y[0]=h*(frombf(a)+frombf(b))/2.0;
) G2 o2 U' o) @5 M1 u4 G, K    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++)
! _' v+ [0 N8 I* |" s          { x=a+(i+0.5)*h;
            p=p+frombf(x);
          }
        p=(y[0]+h*p)/2.0;
5 v/ z2 h* p0 z) E6 E( y        s=1.0;
/ s5 ]% P3 _* B/ v) |5 u. S% d" z" P; I        for (k=1;k&lt;=m;k++)
" {, ?6 x+ S8 D/ Z1 v          { s=4.0*s;
7 `6 |3 H6 j; F) L  |. c' v! F( T# p3 k            q=(s*p-y[k-1])/(s-1.0);
            y[k-1]=p; p=q;
, N! V% _: t; j* n8 l          }
        ep=fabs(q-y[m-1]);
        m=m+1; y[m-1]=q; n=n+n; h=h/2.0;
      }
    return(q);
0 Y' B: J  D0 [2 h3 f" C9 Y  }</P><>呵呵 希望对你有用！！</P>

西西

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

陈了

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