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

西西

<>本人紧急求助，希望各位大虾们伸出你们友情之手，帮我一下。</P>
<>1。牛顿迭代法</P>
" _  y* R3 L# ?" n$ c3 ?$ U' B<>1。全主消元法</P>
<>3。改进平方根法</P>
3 h' C- A- D* I9 y9 d: G1 ?. g<>4。牛顿向前，向后插值法</P>
<>5。加权最小二乘法</P>
5 o+ {8 A, G, p, Q7 K<>6。龙贝格积分法</P>
7 I6 A( F1 B/ R<>拜托了。。。。</P>

loveyuan

<>1。牛顿迭代法</P><>
  #include "stdio.h"
  #include "math.h"
* f- I; \! c3 I7 U  V+ H  int dnewt(x,eps,js)
  int js;
  double *x,eps;
6 }2 C2 z: A! w6 Z  { extern void dnewtf();
    int k,l;
2 S5 o9 _. m: u: l9 f    double y[2],d,p,x0,x1;
" c" N$ v* D, ?& [- h; w$ S    l=js; k=1; x0=*x;
% I8 `4 p" f3 {0 M' x    dnewtf(x0,y);
' u; O7 m/ h  h    d=eps+1.0;
    while ((d&gt;=eps)&amp;&amp;(l!=0))
# b4 B# x+ B  t9 b3 R      { if (fabs(y[1])+1.0==1.0)
          { printf("err\n"); return(-1);}
0 o$ k4 s+ ?$ C! ^& t4 o' Q        x1=x0-y[0]/y[1];
        dnewtf(x1,y);
# u, K+ ~" U/ e: h2 H        d=fabs(x1-x0); p=fabs(y[0]);
& `/ H) k3 |( A% Y        if (p&gt;d) d=p;
        x0=x1; l=l-1;
& G# n8 o8 G& ]9 n( T0 G      }
    *x=x1;
! {+ j/ d3 s3 s; r% |: O7 f    k=js-l;
    return(k);
2 t" Y7 X7 s/ o  }</P><>全主消元法</P><>#include "stdlib.h"
  #include "stdio.h"
/ s+ u9 k" |& V/ v: D. G  int acgas(ar,ai,n,br,bi)
  int n;
! j, I2 u% d, L. i  double ar[],ai[],br[],bi[];
  { int *js,l,k,i,j,is,u,v;
7 e1 q* A8 L  a4 P9 Z1 r    double p,q,s,d;
+ d9 @' D" P1 ]% s) o    js=malloc(n*sizeof(int));
    for (k=0;k&lt;=n-2;k++)
! o2 O! [3 z* H7 E1 t/ A1 C      { d=0.0;
        for (i=k;i&lt;=n-1;i++)
2 y# L# ~5 `3 c        for (j=k;j&lt;=n-1;j++)
* A& n! q# x8 \4 r; u( ?/ [: ?          { u=i*n+j;
2 P8 x$ M& j- z7 h            p=ar*ar+ai*ai;
            if (p&gt;d) {d=p;js[k]=j;is=i;}
          }
0 f, [* i# H0 K2 X( X        if (d+1.0==1.0)
          { free(js); printf("err**fail\n");
            return(0);
          }
        if (is!=k)
* z* T6 K. H) |- w9 ], P9 L          { for (j=k;j&lt;=n-1;j++)
              { u=k*n+j; v=is*n+j;
                p=ar; ar=ar[v]; ar[v]=p;
# `) ~7 s4 Q9 {2 }                p=ai; ai=ai[v]; ai[v]=p;
* `. Q% N' K! x( g0 |' M              }
8 H0 I# {4 `: D" ?- ]            p=br[k]; br[k]=br[is]; br[is]=p;
            p=bi[k]; bi[k]=bi[is]; bi[is]=p;
          }
        if (js[k]!=k)
& [/ m  K4 q9 H          for (i=0;i&lt;=n-1;i++)
, R* f6 B. `5 n  Z$ u5 I            { u=i*n+k; v=i*n+js[k];
$ h) h: g# E9 |6 s! |* b              p=ar; ar=ar[v]; ar[v]=p;
              p=ai; ai=ai[v]; ai[v]=p;
% P! I  T0 D7 B' J% c+ |/ y/ s1 G            }
        v=k*n+k;
        for (j=k+1;j&lt;=n-1;j++)
          { u=k*n+j;
4 k7 _- R5 o0 {) k% `" g            p=ar*ar[v]; q=-ai*ai[v];
7 [9 B9 j9 k; E8 g2 m' G; j            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];
  q7 |: J- F  T, a0 Q+ X- E8 W# ~        s=(ar[v]-ai[v])*(br[k]+bi[k]);
$ e2 O' G) ?4 G: d4 M7 d2 P+ N        br[k]=(p-q)/d; bi[k]=(s-p-q)/d;
        for (i=k+1;i&lt;=n-1;i++)
          { u=i*n+k;
& h6 F3 _# ~8 ^" A& n' H            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]);
8 [- f8 N8 S6 P- J# ?) Z                ar[l]=ar[l]-p+q;
5 f; R7 g; K! g/ T& D5 Z                ai[l]=ai[l]-s+p+q;
4 U. m4 |- y. T9 Q" {: w. }              }
8 S$ e7 Y3 e) q' v# W8 m( z+ c            p=ar*br[k]; q=ai*bi[k];
            s=(ar+ai)*(br[k]+bi[k]);
            br=br-p+q; bi=bi-s+p+q;
          }
+ L4 u- X4 E: X      }
    u=(n-1)*n+n-1;
0 D0 _- e& J2 A$ v1 B# H: H    d=ar*ar+ai*ai;
3 H5 f3 x" Q: b; ^    if (d+1.0==1.0)
" d& [( S( G2 Z      { free(js); printf("err**fail\n");
2 K" I" X2 Z' W0 j        return(0);
( y8 o0 W( w& X! T      }
    p=ar*br[n-1]; q=-ai*bi[n-1];
    s=(ar-ai)*(br[n-1]+bi[n-1]);
( Z. E3 T. ^5 o2 C6 d: {    br[n-1]=(p-q)/d; bi[n-1]=(s-p-q)/d;
    for (i=n-2;i&gt;=0;i--)
    for (j=i+1;j&lt;=n-1;j++)
9 u9 u7 A4 c6 R8 H% z' Z" G      { u=i*n+j;
' t0 {) C! l. m% d9 i        p=ar*br[j]; q=ai*bi[j];
        s=(ar+ai)*(br[j]+bi[j]);
0 n9 a5 Z4 q" l' g' ]2 @        br=br-p+q;
2 l8 w, K. i/ D. j8 U        bi=bi-s+p+q;
      }
$ f! N" |6 @" H2 W6 s7 _    js[n-1]=n-1;
* M. p( o/ v7 L9 {* {5 B    for (k=n-1;k&gt;=0;k--)
      if (js[k]!=k)
# \+ h  |  e8 F/ w3 s        { p=br[k]; br[k]=br[js[k]]; br[js[k]]=p;
0 ^+ f- |2 D) L/ m' b          p=bi[k]; bi[k]=bi[js[k]]; bi[js[k]]=p;
# _/ K; Y1 s9 f+ k( S        }
2 ~! d$ `: D7 J. e# W" Y. C# c+ X    free(js);
    return(1);
  }</P><>平方根法</P><>#include "math.h"
7 [' ~9 J8 j0 ~; }# r. t! h  #include "stdio.h"
  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))
      { printf("fail\n"); return(-2);}
    a[0]=sqrt(a[0]);
; K- B, U: x- _* ~7 c9 m    for (j=1; j&lt;=n-1; j++) a[j]=a[j]/a[0];
    for (i=1; i&lt;=n-1; i++)
1 h; u6 v# _  j! c. S! x! a9 ?( _  _+ K      { u=i*n+i;
0 L, {- t2 W9 M+ n% g        for (j=1; j&lt;=i; j++)
  U& F( [0 {0 r3 u; i          { v=(j-1)*n+i;
1 ?. E8 U; }3 i1 J9 Q& H            a=a-a[v]*a[v];
2 d# e0 q! C4 V$ _. m$ b          }
4 _5 i7 B% |7 Y: C- W        if ((a+1.0==1.0)||(a&lt;0.0))
4 q& h% @9 p3 a4 C          { printf("fail\n"); return(-2);}
        a=sqrt(a);
) }: T5 H5 N& S4 r% q        if (i!=(n-1))
! N2 d2 e6 ]3 [" A4 N          { for (j=i+1; j&lt;=n-1; j++)
              { v=i*n+j;
                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;
# w7 _0 f" E8 `4 L& U% U9 A              }
          }
! C2 s* m: s/ E/ v      }
    for (j=0; j&lt;=m-1; j++)
      { d[j]=d[j]/a[0];
        for (i=1; i&lt;=n-1; i++)
; T* p  z! Y* R4 U1 d6 J          { u=i*n+i; v=i*m+j;
            for (k=1; k&lt;=i; k++)
              d[v]=d[v]-a[(k-1)*n+i]*d[(k-1)*m+j];
            d[v]=d[v]/a;
          }
      }
    for (j=0; j&lt;=m-1; j++)
      { u=(n-1)*m+j;
$ B9 Z  J* P( ^( w2 h5 w        d=d/a[n*n-1];
% S5 U" O* l& H/ A        for (k=n-1; k&gt;=1; k--)
( B9 o* a' @& Q: U0 j9 Q( |          { u=(k-1)*m+j;
            for (i=k; i&lt;=n-1; i++)
              { v=(k-1)*n+i;
                d=d-a[v]*d[i*m+j];
              }
, o! y, @- _# r- ]: ]            v=(k-1)*n+k-1;
            d=d/a[v];
          }
8 A% ?3 d4 {6 a0 @" |; v7 b9 X1 t      }
    return(2);
  }</P><>牛顿向前，向后插值法</P><>double eelgr(x0,h,n,y,t)
  int n;
* H) S* c/ u  A; A# V" Q  double x0,h,t,y[];
/ R( m4 E5 G- N2 I  { int i,j,k,m;
    double z,s,xi,xj;
% S* F( D. }5 X    float p,q;
$ P& x. K; Q! I0 T5 W; z    z=0.0;
9 W2 i) f4 J4 J. N3 S9 e+ A% i, J1 G    if (n&lt;1) return(z);
  c0 `1 ~) c& x8 m, `    if (n==1) { z=y[0]; return(z);}
' T" Z% l; a6 D    if (n==2)
+ D4 b: Y1 `# J# a0 ?      { z=(y[1]*(t-x0)-y[0]*(t-x0-h))/h;
: \" e* p, |2 h        return(z);
  B+ B5 @* z7 d# B% y2 p      }
% n6 o" x  T; W    if (t&gt;x0)
. e6 f' h! F, j+ e) E. R9 X) H      { p=(t-x0)/h; i=(int)p; q=(float)i;
0 j, g$ z4 _! i: H$ o7 B7 `8 {        if (p&gt;q) i=i+1;
  G7 Y( i: z3 y: Q      }
, l$ \; n7 r" t& x/ O1 R8 I& {$ w7 b    else i=0;
( m/ j# }% h# y# j4 @    k=i-4;
5 K! Q3 K' a' T/ |6 }6 W; n0 e    if (k&lt;0) k=0;
" @" Z0 j, w) l4 X* x    m=i+3;
) K) H9 l* z( N- S5 F2 b& M' R: Y    if (m&gt;n-1) m=n-1;
- w$ Z: |  ]! h5 O9 W9 W    for (i=k;i&lt;=m;i++)
' F4 m9 G/ k0 M) e7 A2 Z1 L      { s=1.0; xi=x0+i*h;
  N, w% p( t% j, C# T        for (j=k; j&lt;=m; j++)
; A5 w' x4 O9 ]4 B$ s% e6 v; l          if (j!=i)
. j& y2 _9 I9 E            { xj=x0+j*h;
              s=s*(t-xj)/(xi-xj);
. t( L% R8 m9 t6 x            }
" v6 L9 O+ v1 m6 l- C" E4 _        z=z+s*y;
$ m9 j9 l8 H) W* [# Z      }
; n% Q% k9 a* ~( X" G    return(z);
  }
+ _0 x# v+ u$ l! U" ~% {( @2 }向前，向后是一样的思想！！</P><>加权最小二乘法</P><>#include "math.h"
  void hpir1(x,y,n,a,m,dt)
  int n,m;
( H' O$ h( O# ^% J  double x[],y[],a[],dt[];
  { int i,j,k;
# P" j6 w5 \- t/ |  w8 t5 i, W1 Z    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;
: Z8 n' y4 }6 q+ L% a, Y    if (m&gt;20) m=20;
    z=0.0;
/ m3 U) E# G3 w: R0 v. n* E    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 t# x* _2 J' i! m1 A* a    c=c/d1; p=p/d1;
5 X4 m4 M. i' g, ^% o    a[0]=c*b[0];
    if (m&gt;1)
      { t[1]=1.0; t[0]=-p;
0 u* W* Z! [( Q& u2 m& D        d2=0.0; c=0.0; g=0.0;
. D/ w  t+ A, t8 {        for (i=0; i&lt;=n-1; i++)
          { q=x-z-p; d2=d2+q*q;
            c=c+y*q;
8 C$ A$ L* b% c/ u& _  c4 }# ^            g=g+(x-z)*q*q;
          }
" m) @2 O( S. C7 X, X' u        c=c/d2; p=g/d2; q=d2/d1;
        d1=d2;
5 G4 g) z- X1 ]& C  T  ?1 p        a[1]=c*t[1]; a[0]=c*t[0]+a[0];
8 r5 m$ y  p" Y' r+ |9 M) I5 p" k      }
2 }; N- w# n" c; ]9 F, \+ j5 q    for (j=2; j&lt;=m-1; j++)
6 R; g0 Y. @3 S9 L      { s[j]=t[j-1];
; o% u5 N& w( i" v& I9 L        s[j-1]=-p*t[j-1]+t[j-2];
% J8 F$ q0 D4 K        if (j&gt;=3)
          for (k=j-2; k&gt;=1; k--)
: L& r* `4 S$ ^            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;
) b1 g# A; A: ]; \% c+ n4 Z& v+ b        for (i=0; i&lt;=n-1; i++)
& p7 a' y  c* w9 T* j4 Z) P          { q=s[j];
            for (k=j-1; k&gt;=0; k--)
2 e) c4 ^  V" o* B- x              q=q*(x-z)+s[k];
            d2=d2+q*q; c=c+y*q;
            g=g+(x-z)*q*q;
          }
% k" d+ M" b  l; c        c=c/d2; p=g/d2; q=d2/d1;
        d1=d2;
4 o4 T# L; R* s9 T) T        a[j]=c*s[j]; t[j]=s[j];
        for (k=j-1; k&gt;=0; k--)
          { a[k]=c*s[k]+a[k];
5 M6 v8 T/ I0 B% u            b[k]=t[k]; t[k]=s[k];
! }$ ~0 Y9 `9 r9 s1 w( a% i9 o          }
      }
    dt[0]=0.0; dt[1]=0.0; dt[2]=0.0;
) m9 |2 O. V" I9 {- [    for (i=0; i&lt;=n-1; i++)
. q9 I. a- D) w- b. B/ f! o      { q=a[m-1];
" |/ f& a( k' \7 C2 a6 h3 u        for (k=m-2; k&gt;=0; k--)
          q=a[k]+q*(x-z);
- C8 Y; K5 V! [% H2 n) ?1 D) b9 r        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);
& e+ T; F0 Z; J      }
    return;
  }</P><>龙贝格积分法</P><>#include "math.h"
  double fromb(a,b,eps)
  double a,b,eps;
  { extern double frombf();
  P: c# d5 ~, Y! u( ?; E  Q5 [    int m,n,i,k;
    double y[10],h,ep,p,x,s,q;
$ M% n9 u6 W1 G    h=b-a;
4 C$ d& G% b" i; k4 O4 `  l- H    y[0]=h*(frombf(a)+frombf(b))/2.0;
    m=1; n=1; ep=eps+1.0;
6 I' m9 o, k7 c6 j0 H    while ((ep&gt;=eps)&amp;&amp;(m&lt;=9))
: u1 F3 f% v* n* k" i7 q      { p=0.0;
( {: K# |+ |$ e& a0 N/ p" @        for (i=0;i&lt;=n-1;i++)
          { x=a+(i+0.5)*h;
1 ^. f$ W/ g) L$ t) c7 A            p=p+frombf(x);
2 k, ^' x) B; @          }
0 H, \' O4 }" D  u        p=(y[0]+h*p)/2.0;
        s=1.0;
        for (k=1;k&lt;=m;k++)
" ~- @9 O& I4 y* Z+ W          { s=4.0*s;
            q=(s*p-y[k-1])/(s-1.0);
3 w5 ^4 J" \2 [' G            y[k-1]=p; p=q;
4 n/ a- A' B: y: x% ~$ ^7 R          }
6 ?( w$ d! J, k        ep=fabs(q-y[m-1]);
: c: s6 Y1 I5 n$ ]5 L" ]' F        m=m+1; y[m-1]=q; n=n+n; h=h/2.0;
      }
    return(q);
  }</P><>呵呵 希望对你有用！！</P>

西西

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

陈了

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