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

西西

<>本人紧急求助，希望各位大虾们伸出你们友情之手，帮我一下。</P>
<>1。牛顿迭代法</P>
/ Z+ j( E; X3 k' v4 N& J<>1。全主消元法</P>
$ ?: }5 X# K, E. M5 K; o<>3。改进平方根法</P>
/ ~% k  H3 e8 r<>4。牛顿向前，向后插值法</P>
7 ~# @5 `; T* Q, b<>5。加权最小二乘法</P>
/ b# _- {" O! U5 `/ v1 Q<>6。龙贝格积分法</P>
' o; N& @% @; b<>拜托了。。。。</P>

loveyuan

<>1。牛顿迭代法</P><>
  #include "stdio.h"
  #include "math.h"
* C. l" A' d, \( @" r- |  int dnewt(x,eps,js)
  int js;
  double *x,eps;
  { extern void dnewtf();
    int k,l;
1 O0 d- o: N: x    double y[2],d,p,x0,x1;
    l=js; k=1; x0=*x;
7 ^9 g- _- _) M( M1 o6 t    dnewtf(x0,y);
    d=eps+1.0;
9 a( y0 x- F3 d2 ^4 }    while ((d&gt;=eps)&amp;&amp;(l!=0))
      { if (fabs(y[1])+1.0==1.0)
, q# }, m6 P) B& a          { printf("err\n"); return(-1);}
2 F/ I7 A4 z6 A- s        x1=x0-y[0]/y[1];
$ d' t8 s( ]- A; q4 x- G. D  z; K        dnewtf(x1,y);
        d=fabs(x1-x0); p=fabs(y[0]);
        if (p&gt;d) d=p;
) a/ w' P$ d' r* ]) f        x0=x1; l=l-1;
      }
    *x=x1;
1 `! n4 O6 N7 ?/ H2 d5 m( D9 W    k=js-l;
& V5 x' a" Q4 T0 _    return(k);
  }</P><>全主消元法</P><>#include "stdlib.h"
  O: X. ?: b# z  #include "stdio.h"
0 R) c+ E! e* q/ F  int acgas(ar,ai,n,br,bi)
  int n;
  double ar[],ai[],br[],bi[];
  { int *js,l,k,i,j,is,u,v;
    double p,q,s,d;
    js=malloc(n*sizeof(int));
, G$ c4 \6 z0 r( ^# @    for (k=0;k&lt;=n-2;k++)
0 {& F& g  @' z; N# ]      { d=0.0;
3 G4 a- ~' ^) H. G8 X# Q        for (i=k;i&lt;=n-1;i++)
        for (j=k;j&lt;=n-1;j++)
. z0 y7 t3 A8 G" h          { u=i*n+j;
            p=ar*ar+ai*ai;
6 }% ?9 V, N! h3 U; {- h% P9 x            if (p&gt;d) {d=p;js[k]=j;is=i;}
* j3 w2 }+ k: A7 b, y5 e6 F          }
* f" e( ~! D6 p* `  E8 I        if (d+1.0==1.0)
+ ?4 e3 @1 V# [          { free(js); printf("err**fail\n");
            return(0);
          }
0 A2 D; `/ w- s* G; C$ X        if (is!=k)
          { for (j=k;j&lt;=n-1;j++)
              { u=k*n+j; v=is*n+j;
( i7 E, U( Q6 p6 \- y* T* _- j: W                p=ar; ar=ar[v]; ar[v]=p;
                p=ai; ai=ai[v]; ai[v]=p;
4 w! L( [" U- j+ I! q$ p              }
            p=br[k]; br[k]=br[is]; br[is]=p;
$ `5 `* A3 I" A            p=bi[k]; bi[k]=bi[is]; bi[is]=p;
. r$ f: u- C8 P( F# B! R          }
/ }' ^: q( c4 @2 j1 G* `        if (js[k]!=k)
) \. O+ J% P5 u$ [          for (i=0;i&lt;=n-1;i++)
! j' ^  b& G, d4 u6 k$ F% I! @            { u=i*n+k; v=i*n+js[k];
( _0 f( Q' h* F, ]8 A% j! R4 {              p=ar; ar=ar[v]; ar[v]=p;
              p=ai; ai=ai[v]; ai[v]=p;
7 A1 A2 S. A! W9 ?8 N+ K3 C) l            }
        v=k*n+k;
        for (j=k+1;j&lt;=n-1;j++)
# g! E1 @! P  }; R          { u=k*n+j;
1 ?3 ~- Q4 P$ O            p=ar*ar[v]; q=-ai*ai[v];
8 s1 t6 z4 V$ V2 N  B& r7 e            s=(ar[v]-ai[v])*(ar+ai);
            ar=(p-q)/d; ai=(s-p-q)/d;
          }
        p=br[k]*ar[v]; q=-bi[k]*ai[v];
! y2 K$ h1 H. g7 n" G, ?- g% |        s=(ar[v]-ai[v])*(br[k]+bi[k]);
        br[k]=(p-q)/d; bi[k]=(s-p-q)/d;
        for (i=k+1;i&lt;=n-1;i++)
          { u=i*n+k;
: b3 K; m4 c: O' o# j            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];
0 X: i5 t' N. _) y* j                s=(ar+ai)*(ar[v]+ai[v]);
8 e1 k3 D% O  a/ K  n                ar[l]=ar[l]-p+q;
                ai[l]=ai[l]-s+p+q;
, a2 w2 R  `& H4 S% S              }
            p=ar*br[k]; q=ai*bi[k];
            s=(ar+ai)*(br[k]+bi[k]);
: ?& y; G" `$ D* x9 Z            br=br-p+q; bi=bi-s+p+q;
          }
5 N# n% n# g: Y  |      }
; o# v) E! X5 @* G+ P: k    u=(n-1)*n+n-1;
    d=ar*ar+ai*ai;
& ^6 W3 J$ _7 T% D8 @    if (d+1.0==1.0)
- e/ j% {0 r! |) L  J" c      { free(js); printf("err**fail\n");
        return(0);
4 J1 F# g# B& z6 I, k. a6 r6 Y      }
1 D' T0 k9 b, e" [) U    p=ar*br[n-1]; q=-ai*bi[n-1];
0 |, t$ t; E% D7 A! r" o    s=(ar-ai)*(br[n-1]+bi[n-1]);
3 G4 c/ x# M  h1 V    br[n-1]=(p-q)/d; bi[n-1]=(s-p-q)/d;
& L1 z. y% M9 X/ m3 O8 r    for (i=n-2;i&gt;=0;i--)
+ E  m# `: z9 Z    for (j=i+1;j&lt;=n-1;j++)
- X# [0 f. L: S  o" e      { u=i*n+j;
9 `' q- \6 _1 T  O9 h& }        p=ar*br[j]; q=ai*bi[j];
        s=(ar+ai)*(br[j]+bi[j]);
( Q! n* p0 S% e& Y+ D        br=br-p+q;
        bi=bi-s+p+q;
      }
    js[n-1]=n-1;
4 _% G. C3 Q# N7 I/ c    for (k=n-1;k&gt;=0;k--)
* d- D: {& c2 g, t: y. j      if (js[k]!=k)
* s6 ^: r9 N0 M- ~' v/ \! z% g4 e        { p=br[k]; br[k]=br[js[k]]; br[js[k]]=p;
          p=bi[k]; bi[k]=bi[js[k]]; bi[js[k]]=p;
# ~* ^* I$ R5 A9 s$ P        }
    free(js);
    return(1);
3 o" g& h9 H$ E4 D* B$ L  }</P><>平方根法</P><>#include "math.h"
  #include "stdio.h"
  int achol(a,n,m,d)
$ h' F7 @! }/ O8 O6 a* g  int n,m;
4 m# f/ ]2 l% p8 t) H  double a[],d[];
  { int i,j,k,u,v;
6 d# d' n1 r# V; B; j' U& I    if ((a[0]+1.0==1.0)||(a[0]&lt;0.0))
      { printf("fail\n"); return(-2);}
  _8 B+ y& B# ^0 E' D: B    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++)
      { u=i*n+i;
        for (j=1; j&lt;=i; j++)
4 Z, v/ {5 z$ r) N5 x          { v=(j-1)*n+i;
            a=a-a[v]*a[v];
4 y; C( S6 y" ?/ Q4 w          }
        if ((a+1.0==1.0)||(a&lt;0.0))
          { printf("fail\n"); return(-2);}
5 g, m1 T+ u3 P# x! n/ R        a=sqrt(a);
/ I3 C7 s, J6 w# J  h& @        if (i!=(n-1))
          { for (j=i+1; j&lt;=n-1; j++)
* G- ]- O( F1 E+ N3 p              { v=i*n+j;
                for (k=1; k&lt;=i; k++)
1 P2 `6 \! x6 h. n3 @                  a[v]=a[v]-a[(k-1)*n+i]*a[(k-1)*n+j];
  G. ]6 i+ Z6 u                a[v]=a[v]/a;
              }
          }
      }
    for (j=0; j&lt;=m-1; j++)
      { d[j]=d[j]/a[0];
$ B' }' F. u  n% U0 L        for (i=1; i&lt;=n-1; i++)
          { u=i*n+i; v=i*m+j;
            for (k=1; k&lt;=i; k++)
3 ^4 }6 i" D+ Z5 F  e( h              d[v]=d[v]-a[(k-1)*n+i]*d[(k-1)*m+j];
            d[v]=d[v]/a;
0 C( B: {, ?5 r! x- I9 `          }
      }
    for (j=0; j&lt;=m-1; j++)
      { u=(n-1)*m+j;
        d=d/a[n*n-1];
        for (k=n-1; k&gt;=1; k--)
          { u=(k-1)*m+j;
) ^, J% h2 V5 z! P% O8 l* V' q            for (i=k; i&lt;=n-1; i++)
              { v=(k-1)*n+i;
                d=d-a[v]*d[i*m+j];
+ _4 l' n8 }* z9 v$ A& `& P& m; C              }
  F- E! o: Z- F* c1 ?! H* `! I3 W            v=(k-1)*n+k-1;
4 O) [3 U) L$ [" s            d=d/a[v];
2 h% r! `) P* |2 Q          }
* S/ l8 O. R6 s' j. G0 G      }
    return(2);
% I, |( r  T5 k& J  }</P><>牛顿向前，向后插值法</P><>double eelgr(x0,h,n,y,t)
, |( e5 }7 u% N9 R* R$ a4 L1 ~$ s  int n;
7 _7 B: d; D- }- C  double x0,h,t,y[];
  { int i,j,k,m;
) e- H0 H/ A( G1 S3 F% c# \    double z,s,xi,xj;
    float p,q;
    z=0.0;
" p9 V, l, h' y; |3 e  r5 C* L- M    if (n&lt;1) return(z);
' K1 n; M: y' V1 m/ ?5 p, T' ?3 d    if (n==1) { z=y[0]; return(z);}
: l) O4 V* ^$ t, Z; x8 Y    if (n==2)
7 X* c! n- }- E0 q9 z      { z=(y[1]*(t-x0)-y[0]*(t-x0-h))/h;
5 x" M0 n1 m& p        return(z);
  E. L9 w: g( x      }
    if (t&gt;x0)
. J7 e( @% G. t      { p=(t-x0)/h; i=(int)p; q=(float)i;
        if (p&gt;q) i=i+1;
8 T( G5 w' E% h8 |6 A3 e/ l      }
: S/ `! E. Y( t  o% _    else i=0;
% L! V' t3 P1 I- ]. R6 P2 g    k=i-4;
    if (k&lt;0) k=0;
    m=i+3;
    if (m&gt;n-1) m=n-1;
    for (i=k;i&lt;=m;i++)
# v) a! H6 C. |( h# G2 |  |      { s=1.0; xi=x0+i*h;
4 e' H1 y3 ]  Z: z( `        for (j=k; j&lt;=m; j++)
6 G5 O+ G) \0 d- w( a          if (j!=i)
            { xj=x0+j*h;
              s=s*(t-xj)/(xi-xj);
            }
        z=z+s*y;
      }
: a5 _& T. b9 {' [0 ?" O/ A    return(z);
  }
向前，向后是一样的思想！！</P><>加权最小二乘法</P><>#include "math.h"
  J* x6 K9 S  u4 F5 t; v2 U  void hpir1(x,y,n,a,m,dt)
3 t  G+ ~( n# r5 f  int n,m;
  double x[],y[],a[],dt[];
: {7 s. k9 B3 U9 d9 o  { int i,j,k;
: G2 v: \; N, B6 r    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;
, t8 c, D& }9 F$ X5 V8 r9 D    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;
" u, X4 {1 C* g% ^+ G1 [1 B) U    for (i=0; i&lt;=n-1; i++)
) H/ h! k' s% i# Z4 w      { p=p+(x-z); c=c+y;}
4 n- J5 q. K8 q& q: Y    c=c/d1; p=p/d1;
6 G+ `/ t& T) @0 l; x% r    a[0]=c*b[0];
+ ?" c0 J& V* R4 `; s5 C    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++)
/ `- {% Y! K; |: z' z$ T7 k" N          { q=x-z-p; d2=d2+q*q;
            c=c+y*q;
' ]9 Y4 {7 L( A2 n            g=g+(x-z)*q*q;
+ I* t7 m* M8 t8 P0 K# C          }
4 b; {0 Q+ f5 Z7 \, R& W- S7 _        c=c/d2; p=g/d2; q=d2/d1;
        d1=d2;
        a[1]=c*t[1]; a[0]=c*t[0]+a[0];
      }
    for (j=2; j&lt;=m-1; j++)
      { s[j]=t[j-1];
1 u: f" q. T6 Q$ V- B4 u        s[j-1]=-p*t[j-1]+t[j-2];
        if (j&gt;=3)
          for (k=j-2; k&gt;=1; k--)
, `8 ?+ D2 O' A            s[k]=-p*t[k]+t[k-1]-q*b[k];
% r5 O% O, [( H' k        s[0]=-p*t[0]-q*b[0];
- \# ]) z0 x2 O3 C& }( i        d2=0.0; c=0.0; g=0.0;
) h5 X$ T4 a0 o) a5 k* R: j$ v4 I4 a        for (i=0; i&lt;=n-1; i++)
, O9 z/ H( R! A1 l( x* k' n          { q=s[j];
            for (k=j-1; k&gt;=0; k--)
+ E" A3 m& j" R              q=q*(x-z)+s[k];
* j5 y' \+ S" c4 R* }+ v  M! e7 [' F            d2=d2+q*q; c=c+y*q;
. {  Z* b& \; l+ X; ]            g=g+(x-z)*q*q;
          }
        c=c/d2; p=g/d2; q=d2/d1;
; U( `5 @( h4 [. [7 u' u        d1=d2;
" s( R) G- _; B, ^; j        a[j]=c*s[j]; t[j]=s[j];
7 X  O. o# D9 p" |- _) z" D& Y. I        for (k=j-1; k&gt;=0; k--)
          { a[k]=c*s[k]+a[k];
            b[k]=t[k]; t[k]=s[k];
; n+ |* o; r" x: C7 m9 P          }
: N+ O2 a: U% w3 Z! s+ t      }
, M; t& ^  m; }+ _1 ^6 b+ m    dt[0]=0.0; dt[1]=0.0; dt[2]=0.0;
" ]+ J; `8 b5 @% g$ C, G  A) q; N    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);
5 ^9 T& e3 t; ~; c! Y        dt[0]=dt[0]+p*p;
        dt[1]=dt[1]+fabs(p);
7 w/ |" Y/ N$ [" M5 ^( ?      }
    return;
6 P' V; l$ j2 D  }</P><>龙贝格积分法</P><>#include "math.h"
  double fromb(a,b,eps)
$ g& F. ~, `, h8 |  z" c  double a,b,eps;
3 K6 P$ E( k% u  { extern double frombf();
: q/ S. [  O$ ?" X! E; C    int m,n,i,k;
    double y[10],h,ep,p,x,s,q;
    h=b-a;
    y[0]=h*(frombf(a)+frombf(b))/2.0;
6 z( f! x- j' M" I9 L& b    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++)
5 m$ n) c4 P7 _, l; S  P5 r  s/ [          { x=a+(i+0.5)*h;
  h9 Y3 ~5 L2 w' N9 R( n            p=p+frombf(x);
' A  {4 [" w5 @$ t. F          }
        p=(y[0]+h*p)/2.0;
        s=1.0;
! O5 c$ p0 h. x6 W1 q: E4 Z6 |* X( C        for (k=1;k&lt;=m;k++)
8 D/ y9 y  |# ~9 y+ |          { s=4.0*s;
            q=(s*p-y[k-1])/(s-1.0);
1 O  C7 k' T* h) @4 r            y[k-1]=p; p=q;
, C; l5 H( O" V2 u% S* P          }
        ep=fabs(q-y[m-1]);
  J$ P% T# c$ ?        m=m+1; y[m-1]=q; n=n+n; h=h/2.0;
: F4 }* k3 C/ s8 I  W      }
    return(q);
  }</P><>呵呵 希望对你有用！！</P>

西西

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

陈了

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