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

西西

<>本人紧急求助，希望各位大虾们伸出你们友情之手，帮我一下。</P>
; X7 ~4 V9 E/ t# b3 E8 P" u9 ~, l<>1。牛顿迭代法</P>
<>1。全主消元法</P>
<>3。改进平方根法</P>
<>4。牛顿向前，向后插值法</P>
<>5。加权最小二乘法</P>
<>6。龙贝格积分法</P>
<>拜托了。。。。</P>

loveyuan

<>1。牛顿迭代法</P><>
. J2 I( u+ x7 R7 A9 D  #include "stdio.h"
7 L, B3 ?3 \) l' m, \8 Q  #include "math.h"
6 M% u( V* ]5 f$ T4 `/ m6 x  int dnewt(x,eps,js)
  int js;
3 G! D+ ]$ `4 m& w% ?9 k+ U  double *x,eps;
  { extern void dnewtf();
3 v& H+ ^; ]8 a4 I4 f4 U4 r    int k,l;
: H) a8 D; I4 W" P. b    double y[2],d,p,x0,x1;
" N, n1 p2 O8 T# O: c+ w* D  [5 @. Y    l=js; k=1; x0=*x;
; N& E4 [; b( F6 u5 s' Z    dnewtf(x0,y);
    d=eps+1.0;
- x" O* G  j: t+ C. ~+ J    while ((d&gt;=eps)&amp;&amp;(l!=0))
      { if (fabs(y[1])+1.0==1.0)
          { printf("err\n"); return(-1);}
4 Q! J  _1 l1 c% ~- F        x1=x0-y[0]/y[1];
! ^' m1 v  p/ D- O+ @/ G9 m5 @; u( J        dnewtf(x1,y);
0 p$ x! E8 A( f        d=fabs(x1-x0); p=fabs(y[0]);
. ]! l% T7 k6 c( k$ |        if (p&gt;d) d=p;
        x0=x1; l=l-1;
      }
2 I1 p5 U- m# o    *x=x1;
    k=js-l;
* D: k' }8 D4 P  ^    return(k);
  }</P><>全主消元法</P><>#include "stdlib.h"
0 o/ _' ^, M+ a" W6 J& C  #include "stdio.h"
6 w% V' U: S: C  q: c  int acgas(ar,ai,n,br,bi)
6 Q% k, x3 K4 h  int n;
  double ar[],ai[],br[],bi[];
  { int *js,l,k,i,j,is,u,v;
    double p,q,s,d;
  T; u" M; a2 {8 r/ Q3 b1 ?    js=malloc(n*sizeof(int));
& v' h  N. o1 Q' Q" b  @    for (k=0;k&lt;=n-2;k++)
1 r9 R% V+ i4 J/ C3 R. `/ v) J      { d=0.0;
        for (i=k;i&lt;=n-1;i++)
        for (j=k;j&lt;=n-1;j++)
          { u=i*n+j;
) i0 o" Y& `, v4 U            p=ar*ar+ai*ai;
            if (p&gt;d) {d=p;js[k]=j;is=i;}
& s' g) {) d$ |3 `/ d) K- F- N0 a          }
& J9 \, z$ K, v6 U4 d        if (d+1.0==1.0)
2 l# v4 v. s5 n+ D          { free(js); printf("err**fail\n");
            return(0);
          }
5 m% i& C+ S+ r: S: ?        if (is!=k)
          { for (j=k;j&lt;=n-1;j++)
& a/ J. @8 a7 I+ n6 L2 w              { u=k*n+j; v=is*n+j;
                p=ar; ar=ar[v]; ar[v]=p;
                p=ai; ai=ai[v]; ai[v]=p;
              }
            p=br[k]; br[k]=br[is]; br[is]=p;
            p=bi[k]; bi[k]=bi[is]; bi[is]=p;
# J' F* i- i& V/ K9 g          }
        if (js[k]!=k)
! j! P" C  d8 t          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;
              p=ai; ai=ai[v]; ai[v]=p;
5 G0 D3 O0 K7 G8 T6 j2 u            }
) t; X3 c2 k7 U3 T. C- D  Z* s        v=k*n+k;
        for (j=k+1;j&lt;=n-1;j++)
          { u=k*n+j;
3 X2 e7 o! _0 S: W) i            p=ar*ar[v]; q=-ai*ai[v];
            s=(ar[v]-ai[v])*(ar+ai);
            ar=(p-q)/d; ai=(s-p-q)/d;
/ Q/ D8 ?' Q* z# I$ {5 X5 O          }
        p=br[k]*ar[v]; q=-bi[k]*ai[v];
3 H/ e' Y7 h/ }0 t. \        s=(ar[v]-ai[v])*(br[k]+bi[k]);
+ [- |0 Y: b9 C* S" U. q, S7 @% X        br[k]=(p-q)/d; bi[k]=(s-p-q)/d;
' X' E, v) L# d        for (i=k+1;i&lt;=n-1;i++)
/ v% U$ U5 I8 f; @# f5 t6 u          { u=i*n+k;
6 L1 j( {) G( K/ L            for (j=k+1;j&lt;=n-1;j++)
              { v=k*n+j; l=i*n+j;
3 `# d  W: r4 ~9 p2 e, S; f4 z                p=ar*ar[v]; q=ai*ai[v];
                s=(ar+ai)*(ar[v]+ai[v]);
4 s6 a& g2 _' F& ^! k+ u! Q5 e                ar[l]=ar[l]-p+q;
                ai[l]=ai[l]-s+p+q;
+ P, r9 b, \& T5 ^% J% J6 I              }
            p=ar*br[k]; q=ai*bi[k];
            s=(ar+ai)*(br[k]+bi[k]);
  }$ x& G1 @" }: ~7 `2 r0 S, b4 S            br=br-p+q; bi=bi-s+p+q;
          }
3 r" L# g% Q1 y! M  [1 ?% O% \1 M* H8 K      }
' V" J) b( C3 [% p  Q. [    u=(n-1)*n+n-1;
  G9 v$ S6 ~, h$ |  B    d=ar*ar+ai*ai;
& r. t- C! G& _( v    if (d+1.0==1.0)
9 H; O/ t4 A3 i# p      { free(js); printf("err**fail\n");
        return(0);
* A6 Y/ ?* m8 L8 Q! o. ?$ z      }
3 p& S( F9 ^! E- _8 [; k1 D    p=ar*br[n-1]; q=-ai*bi[n-1];
9 a$ \2 X* f' K# x/ P# J+ }  ?4 j+ k    s=(ar-ai)*(br[n-1]+bi[n-1]);
) {4 W5 P1 n' S/ S    br[n-1]=(p-q)/d; bi[n-1]=(s-p-q)/d;
    for (i=n-2;i&gt;=0;i--)
6 E1 R" \9 `& k0 ]7 j  J    for (j=i+1;j&lt;=n-1;j++)
, z; T+ z+ M" G      { u=i*n+j;
7 S1 P+ d' d9 P5 v        p=ar*br[j]; q=ai*bi[j];
        s=(ar+ai)*(br[j]+bi[j]);
        br=br-p+q;
        bi=bi-s+p+q;
      }
    js[n-1]=n-1;
    for (k=n-1;k&gt;=0;k--)
9 L8 ]# u8 R$ f* _, I, h4 w$ ^7 K      if (js[k]!=k)
        { p=br[k]; br[k]=br[js[k]]; br[js[k]]=p;
          p=bi[k]; bi[k]=bi[js[k]]; bi[js[k]]=p;
# }* E5 g+ y2 u9 a* x) M& N        }
    free(js);
    return(1);
. g: d3 p2 ?3 q- e( x: S) N% t  }</P><>平方根法</P><>#include "math.h"
  #include "stdio.h"
$ W- P$ I) P2 j( Q8 ?  int achol(a,n,m,d)
  int n,m;
  double a[],d[];
. `7 `5 A2 E! G  { int i,j,k,u,v;
    if ((a[0]+1.0==1.0)||(a[0]&lt;0.0))
3 K  m  j) F, N- h* T' L8 D) X      { printf("fail\n"); return(-2);}
9 T9 [5 w: A% J" o4 C# h2 w    a[0]=sqrt(a[0]);
% A- m1 C: e" k' M6 J3 p& [    for (j=1; j&lt;=n-1; j++) a[j]=a[j]/a[0];
; ^8 z4 E7 O; C' \5 K3 \    for (i=1; i&lt;=n-1; i++)
      { u=i*n+i;
; P( z7 Z! y- \1 h/ e0 X" i, U' J# U        for (j=1; j&lt;=i; j++)
          { v=(j-1)*n+i;
            a=a-a[v]*a[v];
- J4 b) m1 @7 G1 i          }
        if ((a+1.0==1.0)||(a&lt;0.0))
- k/ R- S) a$ `2 n, S          { printf("fail\n"); return(-2);}
) v  o4 t6 N0 u) r+ {        a=sqrt(a);
: k1 D2 @8 Z/ \& y        if (i!=(n-1))
* j2 O; [% g$ z6 F* J          { for (j=i+1; j&lt;=n-1; j++)
) K5 P  q6 M: ?              { v=i*n+j;
, V# T* m, Y, ?1 L) m1 s6 E; d                for (k=1; k&lt;=i; k++)
" @4 l: B# V% X$ [3 P' Q                  a[v]=a[v]-a[(k-1)*n+i]*a[(k-1)*n+j];
, U) ~$ Y! }' @) r% L5 ~8 v                a[v]=a[v]/a;
              }
& q& i; J$ p# g/ b; v          }
      }
    for (j=0; j&lt;=m-1; j++)
      { d[j]=d[j]/a[0];
  n5 d0 H* t. q* e8 r1 R, t        for (i=1; i&lt;=n-1; i++)
          { u=i*n+i; v=i*m+j;
" J; e/ X! h, Z0 ?) F            for (k=1; k&lt;=i; k++)
              d[v]=d[v]-a[(k-1)*n+i]*d[(k-1)*m+j];
1 @: p/ x8 N# [: f+ S. B  l( z/ j) e            d[v]=d[v]/a;
8 S" b, E& |" {; F          }
      }
    for (j=0; j&lt;=m-1; j++)
- `% P" [% E* ?$ H      { u=(n-1)*m+j;
        d=d/a[n*n-1];
        for (k=n-1; k&gt;=1; k--)
          { u=(k-1)*m+j;
3 i5 Q/ T4 H$ G- K5 U, W            for (i=k; i&lt;=n-1; i++)
; A( c2 R6 {. ^4 @6 H) l, k              { v=(k-1)*n+i;
                d=d-a[v]*d[i*m+j];
. h% B; R2 P# D4 D& z              }
            v=(k-1)*n+k-1;
0 [  g* Z7 V& ~% d. ~# h  V            d=d/a[v];
" f; v' d. O0 `8 X          }
$ T# @, z" w4 P% Y* f      }
* X, M) i5 ?8 A! {' N+ G8 |) g    return(2);
. c* W. k4 a. ]+ y3 n  }</P><>牛顿向前，向后插值法</P><>double eelgr(x0,h,n,y,t)
  int n;
  double x0,h,t,y[];
  x6 R8 Z* D2 O% M! a' ?  { int i,j,k,m;
    double z,s,xi,xj;
    float p,q;
    z=0.0;
    if (n&lt;1) return(z);
4 e% k0 j% }2 j( X1 ]    if (n==1) { z=y[0]; return(z);}
  |& ~+ ]! R( g3 m    if (n==2)
9 R1 W+ L+ K8 w- r' G& T- }      { 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;
2 o+ k  t& Z) @5 P5 W" Y; b& I2 c! r& ~) W        if (p&gt;q) i=i+1;
% z2 n* l4 z  ~& F* L! z2 l      }
    else i=0;
' a% w0 n, ~+ p6 R& A6 H$ S5 K    k=i-4;
4 S( u) p7 W+ G- i% |; L    if (k&lt;0) k=0;
    m=i+3;
    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++)
6 K4 L: G6 R, p3 ]2 p7 {' U          if (j!=i)
8 j. D, G" v% {0 \" W6 T7 I9 R            { xj=x0+j*h;
              s=s*(t-xj)/(xi-xj);
            }
        z=z+s*y;
: ^1 q9 D5 Q7 H  G- r      }
2 d+ q& l- L- ^2 c3 u$ Y* d    return(z);
4 B9 n$ j- I/ h0 t9 m4 ?% C: D  }
7 W+ u  p+ z' y" c8 ]7 K: Y$ D向前，向后是一样的思想！！</P><>加权最小二乘法</P><>#include "math.h"
3 |0 r( Y* d& M2 {4 z* @, @  void hpir1(x,y,n,a,m,dt)
  int n,m;
  double x[],y[],a[],dt[];
  { int i,j,k;
    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;
. f0 C" I% d, G    if (m&gt;n) m=n;
    if (m&gt;20) m=20;
    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++)
! A+ ^- H( w+ y! M; S      { p=p+(x-z); c=c+y;}
    c=c/d1; p=p/d1;
9 z2 C  X/ Z5 G    a[0]=c*b[0];
    if (m&gt;1)
; }$ F& [" N' g$ Z      { t[1]=1.0; t[0]=-p;
        d2=0.0; c=0.0; g=0.0;
& ~7 m0 s% S! s& X+ }- A  A* l        for (i=0; i&lt;=n-1; i++)
          { q=x-z-p; d2=d2+q*q;
1 J; W" n9 o( f$ Z6 H) I( a            c=c+y*q;
            g=g+(x-z)*q*q;
          }
        c=c/d2; p=g/d2; q=d2/d1;
! h4 x* U7 n: W5 P7 w7 z$ z        d1=d2;
3 d2 [" \0 l6 u) e* M7 B        a[1]=c*t[1]; a[0]=c*t[0]+a[0];
& M$ _  j$ w& J4 o4 i/ n; g      }
    for (j=2; j&lt;=m-1; j++)
3 d  g- i2 B6 N& Y6 H. @6 L# b" w+ C      { s[j]=t[j-1];
        s[j-1]=-p*t[j-1]+t[j-2];
        if (j&gt;=3)
          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];
- D" h0 h8 ^. f' h        d2=0.0; c=0.0; g=0.0;
        for (i=0; i&lt;=n-1; i++)
- Z8 O+ k' F7 X6 u+ E) [8 A7 F- j          { q=s[j];
            for (k=j-1; k&gt;=0; k--)
* c- e: _' ~6 W2 x. F              q=q*(x-z)+s[k];
# u& v6 \% n! V( P# u5 _* [            d2=d2+q*q; c=c+y*q;
$ U& @! l7 A. }9 ?  B- j0 s7 q            g=g+(x-z)*q*q;
, N. W3 O- {# L9 e1 r; ?) z3 C8 K) J          }
        c=c/d2; p=g/d2; q=d2/d1;
        d1=d2;
3 v5 e9 y3 u$ J& ]* ^5 B1 P! h        a[j]=c*s[j]; t[j]=s[j];
7 P6 r$ `& r4 r2 W+ V# K+ T; @        for (k=j-1; k&gt;=0; k--)
          { a[k]=c*s[k]+a[k];
/ Z2 A5 L8 m9 S6 x. A  `            b[k]=t[k]; t[k]=s[k];
+ X7 K" @4 p9 B& X& u& h  Y3 {          }
4 j  [% b; b, _5 ~$ Q      }
    dt[0]=0.0; dt[1]=0.0; dt[2]=0.0;
    for (i=0; i&lt;=n-1; i++)
$ w0 `. n" F" i" M) i/ l. r+ l      { q=a[m-1];
# L0 o' n9 k, ?6 M$ o8 Y        for (k=m-2; k&gt;=0; k--)
4 @" ]/ C& b! {$ `7 c" |- j' _          q=a[k]+q*(x-z);
/ W* F% s0 v6 x! C2 T" N- l1 K4 P        p=q-y;
, e( Z3 }8 G: w* {        if (fabs(p)&gt;dt[2]) dt[2]=fabs(p);
        dt[0]=dt[0]+p*p;
        dt[1]=dt[1]+fabs(p);
& H1 y. l. s/ m$ e7 Q% t( t5 ?      }
  N$ x# e9 Y: g! f0 d    return;
  }</P><>龙贝格积分法</P><>#include "math.h"
  double fromb(a,b,eps)
  double a,b,eps;
, O) e+ M5 T) w. i  { extern double frombf();
    int m,n,i,k;
9 B! M+ l! f- ?2 `. V- {! I    double y[10],h,ep,p,x,s,q;
    h=b-a;
    y[0]=h*(frombf(a)+frombf(b))/2.0;
( n+ ?& j3 E% \: U# P2 l; j9 i    m=1; n=1; ep=eps+1.0;
    while ((ep&gt;=eps)&amp;&amp;(m&lt;=9))
8 B/ r; f0 B9 p- z7 ^      { p=0.0;
        for (i=0;i&lt;=n-1;i++)
          { x=a+(i+0.5)*h;
            p=p+frombf(x);
          }
# h) B0 E+ I  F, B" A# W: g        p=(y[0]+h*p)/2.0;
8 h. r4 \2 _1 j) s: ^6 |  F4 L        s=1.0;
& c5 ?, e% G& [6 ?. v        for (k=1;k&lt;=m;k++)
          { s=4.0*s;
' U) m, }+ o, M( u            q=(s*p-y[k-1])/(s-1.0);
- k# x, b4 n, V* ^, j            y[k-1]=p; p=q;
" f/ }* C' C7 e2 ]$ N0 d. {5 }          }
8 X8 \. F" D) C: h7 m! j        ep=fabs(q-y[m-1]);
; d4 c9 |1 @) _( x3 T  s# c        m=m+1; y[m-1]=q; n=n+n; h=h/2.0;
$ M7 p2 ~  i  O8 Q      }
    return(q);
  }</P><>呵呵 希望对你有用！！</P>

西西

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

陈了

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