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

西西

<>本人紧急求助，希望各位大虾们伸出你们友情之手，帮我一下。</P>
<>1。牛顿迭代法</P>
<>1。全主消元法</P>
<>3。改进平方根法</P>
% Q- ~8 R- I5 M' @2 ~<>4。牛顿向前，向后插值法</P>
, F. C# m+ Q& y" f: r6 H<>5。加权最小二乘法</P>
<>6。龙贝格积分法</P>
<>拜托了。。。。</P>

loveyuan

<>1。牛顿迭代法</P><>
! P, z2 l; n- ]- k! m- `7 P9 }  #include "stdio.h"
  #include "math.h"
- @: e. E! ?3 ~  int dnewt(x,eps,js)
  int js;
  double *x,eps;
  N  Z2 _0 T: I5 n, c7 {9 t: U, m7 d  { extern void dnewtf();
, D4 |: B/ |: Y; `' B    int k,l;
$ G& S& i0 E0 `2 o$ k0 \1 r    double y[2],d,p,x0,x1;
    l=js; k=1; x0=*x;
/ F' v, l1 a' K! E/ f* n    dnewtf(x0,y);
    d=eps+1.0;
+ o$ J; U! z2 t6 x5 W5 t4 S    while ((d&gt;=eps)&amp;&amp;(l!=0))
      { if (fabs(y[1])+1.0==1.0)
          { printf("err\n"); return(-1);}
. A; A* l9 }" [! M* u3 H* w        x1=x0-y[0]/y[1];
6 J! I5 S% {# W; e        dnewtf(x1,y);
5 r9 \) D0 z) {% W7 Y$ [        d=fabs(x1-x0); p=fabs(y[0]);
  Z8 n! n: [, @        if (p&gt;d) d=p;
        x0=x1; l=l-1;
4 G$ {, s; ?/ c1 l      }
  M  q  i9 N) k0 K    *x=x1;
& q, q  l! r- c    k=js-l;
6 K+ o& D$ C7 o8 u; K4 A    return(k);
3 w' C4 r9 N8 z- b7 z  }</P><>全主消元法</P><>#include "stdlib.h"
' b, z6 |3 |0 U' C4 P; w1 K: i! c# c  #include "stdio.h"
  int acgas(ar,ai,n,br,bi)
  int n;
  double ar[],ai[],br[],bi[];
& [0 C- y) y! m" m+ T  { int *js,l,k,i,j,is,u,v;
# \) F. b: ~- Z% e/ g    double p,q,s,d;
    js=malloc(n*sizeof(int));
, X3 ?  T3 x( m% P$ l/ y    for (k=0;k&lt;=n-2;k++)
      { d=0.0;
5 G& C$ G/ a2 _1 ?" s        for (i=k;i&lt;=n-1;i++)
        for (j=k;j&lt;=n-1;j++)
          { u=i*n+j;
            p=ar*ar+ai*ai;
            if (p&gt;d) {d=p;js[k]=j;is=i;}
          }
        if (d+1.0==1.0)
          { free(js); printf("err**fail\n");
            return(0);
  F, E8 C; K% w7 c. k- a. p# O6 s          }
8 b* F- X& T0 @/ l3 P; R; C        if (is!=k)
& L/ Q; x9 b+ [4 o/ K( I1 v          { for (j=k;j&lt;=n-1;j++)
              { u=k*n+j; v=is*n+j;
                p=ar; ar=ar[v]; ar[v]=p;
                p=ai; ai=ai[v]; ai[v]=p;
% t' k, [7 n/ h% t2 V              }
            p=br[k]; br[k]=br[is]; br[is]=p;
            p=bi[k]; bi[k]=bi[is]; bi[is]=p;
) l' I* ]8 y; d& N. C/ F          }
        if (js[k]!=k)
6 B& h! b2 T( ?: }9 i          for (i=0;i&lt;=n-1;i++)
/ i& _* l& C' K- i' D8 L# p            { 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;
2 N* S' w5 u! A# H4 E. E/ G' _0 Z            }
        v=k*n+k;
7 T' b+ j& C5 I' s        for (j=k+1;j&lt;=n-1;j++)
          { u=k*n+j;
            p=ar*ar[v]; q=-ai*ai[v];
            s=(ar[v]-ai[v])*(ar+ai);
            ar=(p-q)/d; ai=(s-p-q)/d;
  P" {; c! z; j' U9 ^  E0 {          }
        p=br[k]*ar[v]; q=-bi[k]*ai[v];
        s=(ar[v]-ai[v])*(br[k]+bi[k]);
2 |7 g, v5 y6 _. m! R, _4 \3 I0 z        br[k]=(p-q)/d; bi[k]=(s-p-q)/d;
        for (i=k+1;i&lt;=n-1;i++)
          { u=i*n+k;
! g- Y6 a9 X  G! p. v/ k2 J            for (j=k+1;j&lt;=n-1;j++)
2 c2 }% j; a8 l. P              { v=k*n+j; l=i*n+j;
                p=ar*ar[v]; q=ai*ai[v];
                s=(ar+ai)*(ar[v]+ai[v]);
# b! R( q* A  t. k, U                ar[l]=ar[l]-p+q;
% W% {" b. ?5 m2 k# A4 J% B1 ^                ai[l]=ai[l]-s+p+q;
              }
& ?0 r% j- a# c, k% O. `$ {' b            p=ar*br[k]; q=ai*bi[k];
! q- N$ P! y) @' W            s=(ar+ai)*(br[k]+bi[k]);
* E/ t% ~4 B: [% x& o% M            br=br-p+q; bi=bi-s+p+q;
) v: w* ?5 E: A4 V7 M- F: ]          }
* _9 G# U; q- p9 l  ^      }
) S5 a  k' \6 W' {    u=(n-1)*n+n-1;
    d=ar*ar+ai*ai;
    if (d+1.0==1.0)
      { free(js); printf("err**fail\n");
6 u1 b: o+ h, {; `/ N, i4 W        return(0);
      }
    p=ar*br[n-1]; q=-ai*bi[n-1];
    s=(ar-ai)*(br[n-1]+bi[n-1]);
. i- b; J' y$ V" x( k( s    br[n-1]=(p-q)/d; bi[n-1]=(s-p-q)/d;
( W( K, Y( N3 t; `9 O4 f. Z    for (i=n-2;i&gt;=0;i--)
/ D, C4 e  X- F$ v% ]3 x$ N    for (j=i+1;j&lt;=n-1;j++)
      { u=i*n+j;
& d, U4 T9 ^" l( a        p=ar*br[j]; q=ai*bi[j];
& b7 |5 y7 B: W        s=(ar+ai)*(br[j]+bi[j]);
        br=br-p+q;
# W5 t6 V1 a: {$ c6 V$ @        bi=bi-s+p+q;
8 A  \0 f0 q4 _' z, j9 X      }
    js[n-1]=n-1;
    for (k=n-1;k&gt;=0;k--)
      if (js[k]!=k)
. c/ Z1 Q, S( a6 ]8 X% |        { p=br[k]; br[k]=br[js[k]]; br[js[k]]=p;
& ?2 G! Q- @; {% Q1 t          p=bi[k]; bi[k]=bi[js[k]]; bi[js[k]]=p;
        }
0 |* H" g$ P( ^    free(js);
    return(1);
7 I0 S) v& U) k5 `1 x  }</P><>平方根法</P><>#include "math.h"
  #include "stdio.h"
  int achol(a,n,m,d)
( E& u+ O) H3 s8 r; i  A  int n,m;
- z+ f/ @2 v/ k+ ~# n  F. m  double a[],d[];
# Y% x) v8 T9 k  { int i,j,k,u,v;
    if ((a[0]+1.0==1.0)||(a[0]&lt;0.0))
& h3 A( P  }# `0 Y      { printf("fail\n"); return(-2);}
$ I/ \+ z$ P6 P6 N3 t6 g    a[0]=sqrt(a[0]);
. F" d4 E5 R$ M4 T/ h' I- g    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++)
          { v=(j-1)*n+i;
            a=a-a[v]*a[v];
' I- t2 {4 ?, U# L! ~! z. r          }
        if ((a+1.0==1.0)||(a&lt;0.0))
! S3 a5 M2 [' ?* k, }          { printf("fail\n"); return(-2);}
. X) `; M$ M2 g6 x        a=sqrt(a);
0 z% E# ^9 J0 c" o6 G        if (i!=(n-1))
          { for (j=i+1; j&lt;=n-1; j++)
0 ?' K* L6 `/ _( }) i: {! J2 Q              { v=i*n+j;
/ C, }# \2 k, f5 [" 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;
              }
) u% B/ \# ~% v6 Y8 I- [          }
" |" @) W# g$ B! p+ n      }
* E- g& X4 e- A& S* A  B3 n    for (j=0; j&lt;=m-1; j++)
) P: c4 M5 i4 s* i0 x      { d[j]=d[j]/a[0];
        for (i=1; i&lt;=n-1; i++)
          { 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];
- w3 ?* j8 o2 s            d[v]=d[v]/a;
          }
      }
    for (j=0; j&lt;=m-1; j++)
      { u=(n-1)*m+j;
5 S/ f+ f" |4 O        d=d/a[n*n-1];
/ C+ s9 x' a9 t. A' j+ S        for (k=n-1; k&gt;=1; k--)
          { u=(k-1)*m+j;
' o& w. H/ T8 ]: y            for (i=k; i&lt;=n-1; i++)
2 s2 ^2 C4 _$ V              { v=(k-1)*n+i;
                d=d-a[v]*d[i*m+j];
  A7 }% `0 U% a3 O4 J3 ~              }
            v=(k-1)*n+k-1;
            d=d/a[v];
* {& s3 ?6 O0 v3 i: {0 t          }
3 I6 H& Z" ^2 r- w) L( U      }
    return(2);
$ V7 R6 k3 A  p1 M6 a4 y* t  }</P><>牛顿向前，向后插值法</P><>double eelgr(x0,h,n,y,t)
* @5 Z' x) M6 i/ d/ S5 x  int n;
  double x0,h,t,y[];
  { int i,j,k,m;
    double z,s,xi,xj;
' K- r- p! H8 D    float p,q;
1 P( s/ I8 A- _$ j    z=0.0;
) c8 N9 }: r2 f; t4 b* |    if (n&lt;1) return(z);
    if (n==1) { z=y[0]; return(z);}
8 ]4 J. C' R) l) ^: ~/ @! z    if (n==2)
$ P9 _7 ~8 O1 D* l% ]      { z=(y[1]*(t-x0)-y[0]*(t-x0-h))/h;
        return(z);
' i# M5 T! p! F( n* H8 k4 H      }
    if (t&gt;x0)
, h3 Z6 |5 l; A  y      { p=(t-x0)/h; i=(int)p; q=(float)i;
) T  l1 ]% c$ ]+ C        if (p&gt;q) i=i+1;
7 a' f. {: _2 y+ J) Q      }
( v. k* l& x* K- x8 E    else i=0;
* M3 z4 L2 `* o    k=i-4;
    if (k&lt;0) k=0;
    m=i+3;
. O3 x) Y# L0 Y, y* Q; Q    if (m&gt;n-1) m=n-1;
, e# R& g% g8 q7 J. p4 l0 z    for (i=k;i&lt;=m;i++)
      { s=1.0; xi=x0+i*h;
        for (j=k; j&lt;=m; j++)
  G1 {' Q! f4 u# X& \          if (j!=i)
            { xj=x0+j*h;
              s=s*(t-xj)/(xi-xj);
- z1 G0 \& c( P3 T9 f% q. X) {7 a/ U            }
/ U: \. {! s; v; l        z=z+s*y;
( _: v  O3 u. {7 ~5 ^      }
# @! P4 x8 s# Y7 [  C9 X    return(z);
  }
' I5 ?: _. |+ P+ w( ]6 \向前，向后是一样的思想！！</P><>加权最小二乘法</P><>#include "math.h"
+ x. b1 p7 m4 w) |" k* G  void hpir1(x,y,n,a,m,dt)
) _8 B+ g) q+ J5 M" f  int n,m;
3 j) ~2 K  U. w  double x[],y[],a[],dt[];
  { int i,j,k;
3 f9 T5 g  V* y, }5 G7 M    double z,p,c,g,q,d1,d2,s[20],t[20],b[20];
, M* d1 I) E7 L9 m$ o    for (i=0; i&lt;=m-1; i++) a=0.0;
    if (m&gt;n) m=n;
! e0 m# V; u! v( w0 L; ]    if (m&gt;20) m=20;
5 L) o3 ~. @0 x" I5 v    z=0.0;
+ f) I9 Z. L7 X" d5 l$ }    for (i=0; i&lt;=n-1; i++) z=z+x/(1.0*n);
. ]; R  M8 q2 t0 K& p5 Y    b[0]=1.0; d1=1.0*n; p=0.0; c=0.0;
' N8 }* o% p6 |9 x) d$ ?$ R    for (i=0; i&lt;=n-1; i++)
. L  N. ?* E9 V  ?  }+ P: w- q      { p=p+(x-z); c=c+y;}
, v' D6 a& O/ k    c=c/d1; p=p/d1;
2 y9 K2 ]/ h7 O: r; W    a[0]=c*b[0];
! b2 }: D6 t- @: i  [0 P1 a    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++)
( `! n3 G% u# D6 Q' w$ |, ^          { q=x-z-p; d2=d2+q*q;
            c=c+y*q;
            g=g+(x-z)*q*q;
* D& {' }, T( d6 a          }
. q3 {* r3 Q" d- Z9 i( _        c=c/d2; p=g/d2; q=d2/d1;
        d1=d2;
0 W- e! V' v# a+ Q- d* s5 P$ e        a[1]=c*t[1]; a[0]=c*t[0]+a[0];
      }
2 }; F# J* O( r' m" x+ y    for (j=2; j&lt;=m-1; j++)
      { s[j]=t[j-1];
, I1 [6 E- N7 K/ s4 o: R% q        s[j-1]=-p*t[j-1]+t[j-2];
* O  `2 E0 G% v0 g4 E  `) h        if (j&gt;=3)
8 G$ q. S" ]' m7 t( L- ?# h          for (k=j-2; k&gt;=1; k--)
            s[k]=-p*t[k]+t[k-1]-q*b[k];
) c* z; m/ g. z- Q8 e        s[0]=-p*t[0]-q*b[0];
        d2=0.0; c=0.0; g=0.0;
        for (i=0; i&lt;=n-1; i++)
4 O. [) E5 \# {+ W$ o4 M          { q=s[j];
            for (k=j-1; k&gt;=0; k--)
              q=q*(x-z)+s[k];
            d2=d2+q*q; c=c+y*q;
            g=g+(x-z)*q*q;
          }
        c=c/d2; p=g/d2; q=d2/d1;
        d1=d2;
        a[j]=c*s[j]; t[j]=s[j];
        for (k=j-1; k&gt;=0; k--)
          { a[k]=c*s[k]+a[k];
            b[k]=t[k]; t[k]=s[k];
/ B7 l6 C, y! v7 c0 _) ^; t          }
# r2 r0 y+ ~5 S- I- J      }
" C. E4 r2 W- F% @+ \. T- s. Y    dt[0]=0.0; dt[1]=0.0; dt[2]=0.0;
3 y. W9 x6 `4 ?! |    for (i=0; i&lt;=n-1; i++)
      { q=a[m-1];
: Q3 ^. w6 L6 d7 v5 R4 w! g' e! w        for (k=m-2; k&gt;=0; k--)
          q=a[k]+q*(x-z);
: e/ Y5 [" j! o& ^2 R- M        p=q-y;
7 y: j; Q/ m; n8 |9 P        if (fabs(p)&gt;dt[2]) dt[2]=fabs(p);
        dt[0]=dt[0]+p*p;
        dt[1]=dt[1]+fabs(p);
) e+ O" N- Q% K8 W      }
    return;
  }</P><>龙贝格积分法</P><>#include "math.h"
  double fromb(a,b,eps)
. F; ]. u- O8 A' j# ]1 t: ^  double a,b,eps;
  { extern double frombf();
* p5 _$ g/ @/ Q    int m,n,i,k;
    double y[10],h,ep,p,x,s,q;
( g. G0 t$ q& S, X6 \+ A! ?' v    h=b-a;
. I  h# U$ `) N! C# t* ^# j    y[0]=h*(frombf(a)+frombf(b))/2.0;
$ Z: D0 s9 d. v) w; G. Z5 b; ]    m=1; n=1; ep=eps+1.0;
% c) L) G  Z# ~4 L$ b' M( V/ `    while ((ep&gt;=eps)&amp;&amp;(m&lt;=9))
      { p=0.0;
9 z8 @# y% ?$ E5 ?5 L        for (i=0;i&lt;=n-1;i++)
$ u/ Q+ ]7 T. l3 _+ _! a          { x=a+(i+0.5)*h;
$ R4 S3 q' S' Q, Q            p=p+frombf(x);
* S+ t. U# _1 r: M. q. u          }
        p=(y[0]+h*p)/2.0;
        s=1.0;
$ I: }7 d6 f* }" ?( c" l: L        for (k=1;k&lt;=m;k++)
; e: H; b1 o/ D4 |1 Q* I          { s=4.0*s;
1 |  Q3 E6 l) }4 w  ?2 q% Z7 @2 ?            q=(s*p-y[k-1])/(s-1.0);
            y[k-1]=p; p=q;
          }
; j7 o( H) A/ ~- }+ y+ w& t        ep=fabs(q-y[m-1]);
        m=m+1; y[m-1]=q; n=n+n; h=h/2.0;
      }
    return(q);
) K9 ?$ v! v) B% c4 @  }</P><>呵呵 希望对你有用！！</P>

西西

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

陈了

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