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

西西

<>本人紧急求助，希望各位大虾们伸出你们友情之手，帮我一下。</P>
<>1。牛顿迭代法</P>
<>1。全主消元法</P>
<>3。改进平方根法</P>
8 c! J& b% N/ E6 f9 ]8 Z<>4。牛顿向前，向后插值法</P>
<>5。加权最小二乘法</P>
<>6。龙贝格积分法</P>
<>拜托了。。。。</P>

loveyuan

<>1。牛顿迭代法</P><>
  #include "stdio.h"
  #include "math.h"
4 l0 z. S5 M7 a; O" z1 U  int dnewt(x,eps,js)
  int js;
4 P! B# X5 u' S) ?  double *x,eps;
7 Z1 }/ x" m: B4 Z7 y( d  { extern void dnewtf();
( p. F) F* f3 u- y    int k,l;
    double y[2],d,p,x0,x1;
. u( a* R# i/ R. `! s    l=js; k=1; x0=*x;
    dnewtf(x0,y);
    d=eps+1.0;
5 Z& f# ]0 K" B8 q5 e    while ((d&gt;=eps)&amp;&amp;(l!=0))
# n* G) I0 b& o0 e$ {; g      { if (fabs(y[1])+1.0==1.0)
          { printf("err\n"); return(-1);}
        x1=x0-y[0]/y[1];
        dnewtf(x1,y);
4 |/ s5 `: y, |6 |+ b2 n% X        d=fabs(x1-x0); p=fabs(y[0]);
        if (p&gt;d) d=p;
) y- B+ }& }$ f% \% y* k        x0=x1; l=l-1;
      }
( m! U: J; V% ]8 m    *x=x1;
    k=js-l;
9 [% z' l8 s! R: y6 w. c, R    return(k);
$ J" l2 m, {3 L4 e3 R4 m' Z0 |- K  }</P><>全主消元法</P><>#include "stdlib.h"
  #include "stdio.h"
* O% E- _- ?- e. Q- |/ W  int acgas(ar,ai,n,br,bi)
  int n;
1 m/ F4 `7 ?4 Z8 W' c  double ar[],ai[],br[],bi[];
; [: e5 |8 l5 D! Y" b  { int *js,l,k,i,j,is,u,v;
    double p,q,s,d;
* N; i# T( U  Q/ {( b2 f$ u7 m6 o    js=malloc(n*sizeof(int));
    for (k=0;k&lt;=n-2;k++)
; Q2 F4 s: M- Z' R% y9 f( ?      { d=0.0;
, i1 A$ ~6 l* I7 g' X        for (i=k;i&lt;=n-1;i++)
        for (j=k;j&lt;=n-1;j++)
" b: s& ?3 q1 J- u6 G& t1 Q, K' J          { u=i*n+j;
1 C: w; u$ c' F9 v) d            p=ar*ar+ai*ai;
            if (p&gt;d) {d=p;js[k]=j;is=i;}
+ T0 r! n. i" K          }
        if (d+1.0==1.0)
          { free(js); printf("err**fail\n");
/ J! }" X- m1 g( [4 e( j* r/ {8 X            return(0);
          }
        if (is!=k)
          { for (j=k;j&lt;=n-1;j++)
              { u=k*n+j; v=is*n+j;
! r( U3 f( \7 }: q: \" d                p=ar; ar=ar[v]; ar[v]=p;
, M! D2 d$ n" U9 D$ [3 r                p=ai; ai=ai[v]; ai[v]=p;
; I7 r$ x9 T- x3 y9 H/ d              }
            p=br[k]; br[k]=br[is]; br[is]=p;
            p=bi[k]; bi[k]=bi[is]; bi[is]=p;
          }
$ E" Z- V, p1 I! l5 {9 D8 f        if (js[k]!=k)
1 b6 b$ T$ i3 j! p          for (i=0;i&lt;=n-1;i++)
  w% L) u6 r$ z. o            { u=i*n+k; v=i*n+js[k];
              p=ar; ar=ar[v]; ar[v]=p;
8 @: [# S. i( o( ?" ^              p=ai; ai=ai[v]; ai[v]=p;
            }
        v=k*n+k;
        for (j=k+1;j&lt;=n-1;j++)
& c* L0 Z3 I8 z8 U. h1 K; }          { u=k*n+j;
            p=ar*ar[v]; q=-ai*ai[v];
            s=(ar[v]-ai[v])*(ar+ai);
" V* y7 f% C/ x4 f" \3 M( i            ar=(p-q)/d; ai=(s-p-q)/d;
8 r3 |# _; A% t7 b          }
        p=br[k]*ar[v]; q=-bi[k]*ai[v];
' X1 K; G: S8 O- o        s=(ar[v]-ai[v])*(br[k]+bi[k]);
0 g2 ~/ x  H: m) l. I. m        br[k]=(p-q)/d; bi[k]=(s-p-q)/d;
, H7 T# D8 \! \. t: \% ?3 A, c        for (i=k+1;i&lt;=n-1;i++)
          { u=i*n+k;
2 {; ^, j( a2 n+ j- |            for (j=k+1;j&lt;=n-1;j++)
& ^4 C$ V4 g. ?/ y* B8 k              { v=k*n+j; l=i*n+j;
                p=ar*ar[v]; q=ai*ai[v];
                s=(ar+ai)*(ar[v]+ai[v]);
- y' P! O2 H% S* l0 m% d                ar[l]=ar[l]-p+q;
                ai[l]=ai[l]-s+p+q;
              }
            p=ar*br[k]; q=ai*bi[k];
! f) b8 |' o' K( F, n/ c* n            s=(ar+ai)*(br[k]+bi[k]);
* ^8 I% ^6 K& \5 |1 J            br=br-p+q; bi=bi-s+p+q;
          }
; K% q! K4 e$ C7 e& ]      }
* I) _" u" N! U    u=(n-1)*n+n-1;
    d=ar*ar+ai*ai;
# ?7 @7 E) H- \# f; i# I    if (d+1.0==1.0)
      { free(js); printf("err**fail\n");
        return(0);
+ G$ l! ?8 u# ]; |1 _. N      }
4 i& Z% @% d# h7 f; c    p=ar*br[n-1]; q=-ai*bi[n-1];
    s=(ar-ai)*(br[n-1]+bi[n-1]);
; g% M/ X/ T' c4 d' k" i- K/ ]5 s, z, B    br[n-1]=(p-q)/d; bi[n-1]=(s-p-q)/d;
( T: l8 U3 V2 \) w; Y- U    for (i=n-2;i&gt;=0;i--)
( R) O- v) h& M! @( A2 g    for (j=i+1;j&lt;=n-1;j++)
$ O3 m2 p+ h% G      { u=i*n+j;
        p=ar*br[j]; q=ai*bi[j];
0 D6 `1 r* w- X        s=(ar+ai)*(br[j]+bi[j]);
4 u6 D6 f% g, M0 H        br=br-p+q;
" a/ S2 ]6 A! r( T4 Y        bi=bi-s+p+q;
( N) p. Q! c. I; d) @4 S5 u% d      }
    js[n-1]=n-1;
    for (k=n-1;k&gt;=0;k--)
      if (js[k]!=k)
        { p=br[k]; br[k]=br[js[k]]; br[js[k]]=p;
- L* M# j# @, o          p=bi[k]; bi[k]=bi[js[k]]; bi[js[k]]=p;
( K' j9 Q! k" R+ O: ^2 |" |        }
5 ?* M) n: }6 D# V6 S1 T    free(js);
    return(1);
7 }6 {' X5 ]4 j7 U  }</P><>平方根法</P><>#include "math.h"
) k. u$ F  J- Z+ I  E# e  #include "stdio.h"
  int achol(a,n,m,d)
  int n,m;
. Q( ^' b5 ~( E" B( g: a  double a[],d[];
- [" H2 B. g2 X9 x$ _( n  { int i,j,k,u,v;
& F  a( ~  S! y    if ((a[0]+1.0==1.0)||(a[0]&lt;0.0))
% P* C# t. q" m% J" O  Q      { printf("fail\n"); return(-2);}
    a[0]=sqrt(a[0]);
    for (j=1; j&lt;=n-1; j++) a[j]=a[j]/a[0];
    for (i=1; i&lt;=n-1; i++)
4 Z) T; e) }, N9 W      { u=i*n+i;
  j' g! y/ f. m1 F/ u5 t( p$ [0 `        for (j=1; j&lt;=i; j++)
          { v=(j-1)*n+i;
            a=a-a[v]*a[v];
          }
$ ]7 ?  j- E. k        if ((a+1.0==1.0)||(a&lt;0.0))
          { printf("fail\n"); return(-2);}
        a=sqrt(a);
        if (i!=(n-1))
          { for (j=i+1; j&lt;=n-1; j++)
9 Y  T' s- r3 [. |              { v=i*n+j;
  Q# N0 K5 u' v! `' x                for (k=1; k&lt;=i; k++)
                  a[v]=a[v]-a[(k-1)*n+i]*a[(k-1)*n+j];
                a[v]=a[v]/a;
              }
8 ]5 {' n" @' P8 ]0 n. u  _          }
      }
9 }; y" S+ z% R1 W3 G3 v    for (j=0; j&lt;=m-1; j++)
      { d[j]=d[j]/a[0];
9 }1 g9 E: N# y: ?& b        for (i=1; i&lt;=n-1; i++)
; v. A/ g) y3 i$ R( U7 q# N3 ~          { u=i*n+i; v=i*m+j;
            for (k=1; k&lt;=i; k++)
8 u1 `8 }( `, ~; e+ P              d[v]=d[v]-a[(k-1)*n+i]*d[(k-1)*m+j];
            d[v]=d[v]/a;
          }
; {" P7 ^) Z5 G' n      }
+ ]' s3 C) _/ r) B/ K3 O$ \6 q! Q    for (j=0; j&lt;=m-1; j++)
' b* k* I# H' K      { u=(n-1)*m+j;
        d=d/a[n*n-1];
        for (k=n-1; k&gt;=1; k--)
          { u=(k-1)*m+j;
            for (i=k; i&lt;=n-1; i++)
( M9 t. {" i( O3 ~. S7 E% A- }              { v=(k-1)*n+i;
. |9 H" u! }; a/ P9 Y2 T4 J6 E. X                d=d-a[v]*d[i*m+j];
              }
, a  ?! _3 O; S5 V! d6 @! t            v=(k-1)*n+k-1;
            d=d/a[v];
          }
      }
0 g. T. `2 L" X    return(2);
$ f' f  g7 I4 I1 s, {' ~  }</P><>牛顿向前，向后插值法</P><>double eelgr(x0,h,n,y,t)
  int n;
$ q8 v" K; ^7 C+ w  double x0,h,t,y[];
5 L2 k. X" @  L$ l  { int i,j,k,m;
    double z,s,xi,xj;
6 m- e$ [& E4 h. c/ j: |) M2 T    float p,q;
    z=0.0;
8 S' J& y( P" @! o    if (n&lt;1) return(z);
    if (n==1) { z=y[0]; return(z);}
    if (n==2)
. `. t  L* S. J1 t7 ]% i      { z=(y[1]*(t-x0)-y[0]*(t-x0-h))/h;
        return(z);
( N. n6 p% a6 N# ~( ]) A      }
    if (t&gt;x0)
      { p=(t-x0)/h; i=(int)p; q=(float)i;
, R; ~  i3 O- U+ }        if (p&gt;q) i=i+1;
      }
    else i=0;
    k=i-4;
. t6 t" K  T6 \' I" B3 ^& C    if (k&lt;0) k=0;
5 V* }' P# ^8 G8 S+ p    m=i+3;
    if (m&gt;n-1) m=n-1;
' _; f0 r( N# V- f  _2 e& ~, @    for (i=k;i&lt;=m;i++)
      { s=1.0; xi=x0+i*h;
0 M; M) e" r3 I/ l. K1 S        for (j=k; j&lt;=m; j++)
          if (j!=i)
# s( c! u5 I, q3 U* w) @* }            { xj=x0+j*h;
- @9 [6 j! u) B" [3 q$ p              s=s*(t-xj)/(xi-xj);
, P8 y# x  L& s            }
        z=z+s*y;
' w6 d' n+ I- X9 }8 \$ U, N      }
    return(z);
/ f( n8 P# g  H! d+ Y  }
8 r  R* d* \- K- K: |8 h) n向前，向后是一样的思想！！</P><>加权最小二乘法</P><>#include "math.h"
0 q- k) H4 e* h1 V/ M  void hpir1(x,y,n,a,m,dt)
5 i) A7 k/ R& S+ q: A  int n,m;
9 F; O( Z; m  H% F) y; U  double x[],y[],a[],dt[];
; \: U& o2 R: R: `% J  { int i,j,k;
" m' @" i: [. a* q    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;
% O) \& i( k5 S. ^$ Y+ k, E" ~    if (m&gt;20) m=20;
/ Z& m. k8 o) s* e    z=0.0;
    for (i=0; i&lt;=n-1; i++) z=z+x/(1.0*n);
( w6 c+ H9 x4 V5 E  |; Q    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;}
9 A" R1 y( X' ~% T3 l! A    c=c/d1; p=p/d1;
( Q0 r/ E6 M" C    a[0]=c*b[0];
    if (m&gt;1)
      { t[1]=1.0; t[0]=-p;
        d2=0.0; c=0.0; g=0.0;
        for (i=0; i&lt;=n-1; i++)
2 x1 @' J: R) U9 ]          { q=x-z-p; d2=d2+q*q;
            c=c+y*q;
* v( O+ @7 f: D7 d* b            g=g+(x-z)*q*q;
/ G7 e7 B; x. F' _8 @          }
        c=c/d2; p=g/d2; q=d2/d1;
        d1=d2;
        a[1]=c*t[1]; a[0]=c*t[0]+a[0];
0 }! t/ q5 q7 D      }
7 ^: G" F& R4 L5 j* t$ c    for (j=2; j&lt;=m-1; j++)
+ O. b& g! G  z4 q! }' ?1 e      { s[j]=t[j-1];
        s[j-1]=-p*t[j-1]+t[j-2];
3 S) x: _8 U  _0 q3 O        if (j&gt;=3)
          for (k=j-2; k&gt;=1; k--)
            s[k]=-p*t[k]+t[k-1]-q*b[k];
, O$ p( s; u  t1 ]0 w. K        s[0]=-p*t[0]-q*b[0];
        d2=0.0; c=0.0; g=0.0;
. s" z0 ~; |; D1 V        for (i=0; i&lt;=n-1; i++)
          { q=s[j];
5 T! t! }* r+ Q* n, Y7 @* T2 X            for (k=j-1; k&gt;=0; k--)
, S( t) o6 g3 r8 Q$ B              q=q*(x-z)+s[k];
4 c( G1 I1 h5 m% R% q2 E            d2=d2+q*q; c=c+y*q;
* c  Q8 r  D9 k% `+ b4 l% D; A; h            g=g+(x-z)*q*q;
: f9 Z( m" F3 `; ?. p          }
2 U$ H, F5 {! p3 Z0 U        c=c/d2; p=g/d2; q=d2/d1;
1 E* n" S9 P; D( M) o8 k6 _        d1=d2;
        a[j]=c*s[j]; t[j]=s[j];
8 a' E3 q- w* ]. k9 ~6 }        for (k=j-1; k&gt;=0; k--)
1 L) k; n* V" G7 ?/ k( o. z7 ~          { a[k]=c*s[k]+a[k];
            b[k]=t[k]; t[k]=s[k];
1 a# n. n  s/ R% F; Q! T3 K/ ?          }
      }
( q6 t3 h# ~  H" W; ?5 e! W1 w4 S! A    dt[0]=0.0; dt[1]=0.0; dt[2]=0.0;
" h9 ]3 {- n8 X+ z    for (i=0; i&lt;=n-1; i++)
      { q=a[m-1];
        for (k=m-2; k&gt;=0; k--)
          q=a[k]+q*(x-z);
        p=q-y;
        if (fabs(p)&gt;dt[2]) dt[2]=fabs(p);
        dt[0]=dt[0]+p*p;
        dt[1]=dt[1]+fabs(p);
      }
    return;
  N& W2 j0 K5 Z/ O& o7 d  }</P><>龙贝格积分法</P><>#include "math.h"
, ?3 l9 p: u* O/ k3 i% {  double fromb(a,b,eps)
( j8 ]8 w+ j# I! Y  double a,b,eps;
  { extern double frombf();
1 x3 [2 G) F$ B8 z    int m,n,i,k;
    double y[10],h,ep,p,x,s,q;
    h=b-a;
( S6 m2 l- v6 j; Q1 k    y[0]=h*(frombf(a)+frombf(b))/2.0;
    m=1; n=1; ep=eps+1.0;
    while ((ep&gt;=eps)&amp;&amp;(m&lt;=9))
- _2 u" x0 c6 ], ?' ^' A      { p=0.0;
) @- }7 }! A0 {- |  G7 s/ N        for (i=0;i&lt;=n-1;i++)
0 L$ a/ m/ k; ]( p' J4 U          { x=a+(i+0.5)*h;
            p=p+frombf(x);
          }
1 W6 T+ r) `; S4 ?* m        p=(y[0]+h*p)/2.0;
        s=1.0;
        for (k=1;k&lt;=m;k++)
4 m- z& a$ Y7 ]0 o$ O+ v5 h          { s=4.0*s;
% @; @! E4 B' ]& r+ S) @* _            q=(s*p-y[k-1])/(s-1.0);
            y[k-1]=p; p=q;
          }
        ep=fabs(q-y[m-1]);
        m=m+1; y[m-1]=q; n=n+n; h=h/2.0;
% a1 ~3 w4 h& Q- f      }
. C9 G* N) h7 j) ~! m0 H8 Z    return(q);
" z* I* I" ?$ v% L6 L8 I  }</P><>呵呵 希望对你有用！！</P>

西西

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

陈了

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