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

西西

<>本人紧急求助，希望各位大虾们伸出你们友情之手，帮我一下。</P>
5 [% `5 o: I4 L  C, b1 ^& t* V<>1。牛顿迭代法</P>
3 c( g, J* u% n4 }: v* U5 n; l<>1。全主消元法</P>
1 b! U6 g% K( X<>3。改进平方根法</P>
<>4。牛顿向前，向后插值法</P>
<>5。加权最小二乘法</P>
% j/ d* ~' T1 s" _3 Z<>6。龙贝格积分法</P>
) D: [/ Y' N/ N: s' s4 q<>拜托了。。。。</P>

loveyuan

<>1。牛顿迭代法</P><>
  #include "stdio.h"
  #include "math.h"
2 d6 I3 ?9 E* C, j6 @# ?* Y  int dnewt(x,eps,js)
3 @5 N% Z  T8 [  int js;
  double *x,eps;
( ]2 j5 v, a- L& r1 U; {& v  { extern void dnewtf();
    int k,l;
    double y[2],d,p,x0,x1;
    l=js; k=1; x0=*x;
# O+ ~: D3 u) Y7 d. T( V5 R    dnewtf(x0,y);
, G2 J" K$ `  O1 z6 J    d=eps+1.0;
    while ((d&gt;=eps)&amp;&amp;(l!=0))
      { if (fabs(y[1])+1.0==1.0)
4 T7 ~# q0 J) _* P% k          { printf("err\n"); return(-1);}
5 x% P- a3 Q* q        x1=x0-y[0]/y[1];
        dnewtf(x1,y);
        d=fabs(x1-x0); p=fabs(y[0]);
        if (p&gt;d) d=p;
        x0=x1; l=l-1;
      }
* t* c) U# Y6 R" ~8 X    *x=x1;
    k=js-l;
    return(k);
2 V9 q& @, N; ^  }</P><>全主消元法</P><>#include "stdlib.h"
4 i) w% e+ r9 b; d! }0 o4 \  #include "stdio.h"
  int acgas(ar,ai,n,br,bi)
. r9 _* V0 Q& g, `+ ?( R: N  int n;
2 {" p) f6 z  x  double ar[],ai[],br[],bi[];
0 {- g& B, H) f( S  { int *js,l,k,i,j,is,u,v;
# M! j; m8 ~) m    double p,q,s,d;
    js=malloc(n*sizeof(int));
    for (k=0;k&lt;=n-2;k++)
. G" r* [; V: Q/ O( R) e! {' [      { d=0.0;
- O& j: K8 W9 o        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;
$ U. F$ [+ a2 f" z8 {+ |8 `            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);
6 T! c! X7 B. U+ N6 L          }
        if (is!=k)
          { for (j=k;j&lt;=n-1;j++)
              { u=k*n+j; v=is*n+j;
. T7 v1 }2 I. M3 A: X' G/ L' `0 u                p=ar; ar=ar[v]; ar[v]=p;
6 l" q" C: ^! @3 R) E) _, G! B                p=ai; ai=ai[v]; ai[v]=p;
              }
  C* {" p$ a: v6 e! j. ]0 }* I            p=br[k]; br[k]=br[is]; br[is]=p;
; w  @; x7 f/ @+ A" G            p=bi[k]; bi[k]=bi[is]; bi[is]=p;
9 C, K2 ~# `+ t( y! f          }
. ?  D! w  m8 e7 B9 T  g7 ~% N5 D        if (js[k]!=k)
$ G# X, ?, E2 F$ b9 ^          for (i=0;i&lt;=n-1;i++)
* [2 T4 A$ x9 F& k$ I2 L7 D            { u=i*n+k; v=i*n+js[k];
, p8 h8 M( n$ C2 s+ v              p=ar; ar=ar[v]; ar[v]=p;
0 A" ?8 R2 ^9 Y4 E$ _' h              p=ai; ai=ai[v]; ai[v]=p;
            }
        v=k*n+k;
        for (j=k+1;j&lt;=n-1;j++)
          { u=k*n+j;
            p=ar*ar[v]; q=-ai*ai[v];
% [2 r; ]1 G7 P3 V$ |  @            s=(ar[v]-ai[v])*(ar+ai);
" Z' ]7 D2 g/ E/ Q2 z; m% a$ f$ B            ar=(p-q)/d; ai=(s-p-q)/d;
( W, {2 s2 k% f4 j) P4 w          }
& [2 o- j4 f& {2 ^        p=br[k]*ar[v]; q=-bi[k]*ai[v];
7 g% |2 @) q: ~+ R( z2 s5 k/ A        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;
# E7 m# P/ c# b: E2 a            for (j=k+1;j&lt;=n-1;j++)
- |2 @# p  A- E# j+ j              { v=k*n+j; l=i*n+j;
# Y* W# V( d5 p9 B" C$ ^                p=ar*ar[v]; q=ai*ai[v];
2 c6 ?: I4 `- c( m                s=(ar+ai)*(ar[v]+ai[v]);
9 F% @% o- i6 b/ g  }                ar[l]=ar[l]-p+q;
+ |1 W. M1 y& P                ai[l]=ai[l]-s+p+q;
$ W8 J) \4 S3 U  t8 V              }
            p=ar*br[k]; q=ai*bi[k];
, q, Z7 p! N1 U            s=(ar+ai)*(br[k]+bi[k]);
6 r: J$ N; Q& ~- @            br=br-p+q; bi=bi-s+p+q;
          }
3 z: W- C$ G( T$ D) r! y3 E+ a      }
. q8 J1 a1 ^: n9 J. ?, I) v    u=(n-1)*n+n-1;
# G3 n: z# Y: o2 W8 u  D- W    d=ar*ar+ai*ai;
    if (d+1.0==1.0)
      { free(js); printf("err**fail\n");
/ e; S, j) L% S        return(0);
      }
    p=ar*br[n-1]; q=-ai*bi[n-1];
  P* h6 z4 u4 i    s=(ar-ai)*(br[n-1]+bi[n-1]);
; _1 b5 s3 o+ k* ], J" 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++)
4 }3 `. ?) I5 A      { u=i*n+j;
2 k" _' X! _* p& Q8 `/ V        p=ar*br[j]; q=ai*bi[j];
/ @) I7 W$ f/ h7 h% B8 y        s=(ar+ai)*(br[j]+bi[j]);
        br=br-p+q;
: ~3 T& {7 O' }' ^+ a  F! \9 @' y, N        bi=bi-s+p+q;
      }
1 W8 ^1 u4 i2 O0 `( Z3 t6 e5 b3 R    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;
5 N5 Z8 Q7 D2 C          p=bi[k]; bi[k]=bi[js[k]]; bi[js[k]]=p;
        }
    free(js);
    return(1);
1 O" A' l2 `+ F( X4 r( C! d0 f+ v1 r  }</P><>平方根法</P><>#include "math.h"
4 \/ A$ X5 E# A. T  #include "stdio.h"
  int achol(a,n,m,d)
  int n,m;
8 ?! n  l' u% U) U2 H; ?  double a[],d[];
5 ^/ }, Y  Y6 x( f. ?  { int i,j,k,u,v;
    if ((a[0]+1.0==1.0)||(a[0]&lt;0.0))
  r% K7 K" b. V+ z8 o5 p, s' G& `      { printf("fail\n"); return(-2);}
. w* j1 ]2 |7 \" c# ~: e- v$ C    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++)
+ u7 a/ h- y7 p( a      { u=i*n+i;
. w, L; X" e) R+ s        for (j=1; j&lt;=i; j++)
          { v=(j-1)*n+i;
% S6 B" [! w* ~            a=a-a[v]*a[v];
( f1 [" S0 s: W5 A( [  A          }
3 D; M1 N9 K! e, O9 \& \$ N        if ((a+1.0==1.0)||(a&lt;0.0))
- w/ M! F) s) I; W( }          { printf("fail\n"); return(-2);}
        a=sqrt(a);
        if (i!=(n-1))
. \) M0 k3 i$ w) j- n          { for (j=i+1; j&lt;=n-1; j++)
              { v=i*n+j;
7 q5 o( o+ a" R! m) q$ k, E. l9 i                for (k=1; k&lt;=i; k++)
                  a[v]=a[v]-a[(k-1)*n+i]*a[(k-1)*n+j];
: {0 q7 G! M' W% A4 o3 I. b9 u                a[v]=a[v]/a;
$ B/ a. \" o, F6 ?. Y+ {              }
          }
      }
! G6 s' R% n6 H2 P" i9 A    for (j=0; j&lt;=m-1; j++)
      { d[j]=d[j]/a[0];
        for (i=1; i&lt;=n-1; i++)
3 D$ f& v" V( C) ?3 e' |, }- y          { u=i*n+i; v=i*m+j;
# Y& f8 b8 p+ c/ E            for (k=1; k&lt;=i; k++)
5 p9 o- ~5 s4 ]! Q: T3 A* d6 Z              d[v]=d[v]-a[(k-1)*n+i]*d[(k-1)*m+j];
; L# n6 [$ m9 N+ |8 C# j- _( m$ ]            d[v]=d[v]/a;
          }
      }
    for (j=0; j&lt;=m-1; j++)
      { u=(n-1)*m+j;
& J2 t" W7 x. w+ n6 `        d=d/a[n*n-1];
) \# A& B( k! g/ h        for (k=n-1; k&gt;=1; k--)
          { u=(k-1)*m+j;
8 }- h- A# Y1 M$ D- M            for (i=k; i&lt;=n-1; i++)
: F: D9 I  z- T' L  q4 j              { v=(k-1)*n+i;
                d=d-a[v]*d[i*m+j];
              }
! l; B  D7 Q! D1 `- k2 ^            v=(k-1)*n+k-1;
            d=d/a[v];
          }
4 N. n# h5 g# J% I3 n      }
    return(2);
4 _1 \! h4 A2 R; ~8 _4 ~  }</P><>牛顿向前，向后插值法</P><>double eelgr(x0,h,n,y,t)
  int n;
! y6 g) |. \6 ~  }6 F1 Z7 }3 M  double x0,h,t,y[];
8 I7 g! `) t( I  { int i,j,k,m;
    double z,s,xi,xj;
    float p,q;
) ^. |* |$ |3 E' T8 C4 v' i/ t3 Y  X  ~, Q    z=0.0;
5 U9 O/ j0 i$ `9 c( f& f0 q3 Y    if (n&lt;1) return(z);
/ e' l* R8 Y, z! U    if (n==1) { z=y[0]; return(z);}
$ q. T8 T+ Y. s1 D2 c  Y+ g, Y    if (n==2)
  f( ?0 w: B6 K9 ]      { 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;
4 x3 q" h1 z/ A. A        if (p&gt;q) i=i+1;
      }
, T8 L1 [. A4 @! `5 d! G/ r    else i=0;
    k=i-4;
    if (k&lt;0) k=0;
    m=i+3;
    if (m&gt;n-1) m=n-1;
7 M6 F- R" }; z  Z4 d! D/ o+ u: H1 A    for (i=k;i&lt;=m;i++)
: v0 }& k0 p3 M' J7 _      { s=1.0; xi=x0+i*h;
        for (j=k; j&lt;=m; j++)
4 w6 \1 f& z: [* ]) R, S, ^9 i& }6 n          if (j!=i)
, {# t9 e. c1 z9 @( v" r$ Z            { xj=x0+j*h;
$ @9 ?9 m9 I8 c) _% p) F4 p6 G              s=s*(t-xj)/(xi-xj);
9 H# E# g2 E2 B" `            }
4 {3 i% x5 t/ X# q        z=z+s*y;
      }
    return(z);
/ M. g' C& _# e8 S7 c) |  }
向前，向后是一样的思想！！</P><>加权最小二乘法</P><>#include "math.h"
9 T9 }1 J5 o5 j6 L) L0 M2 }  void hpir1(x,y,n,a,m,dt)
% t4 U( U4 O% f( r  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;
9 [  T1 C/ \( K2 M  s! I& }    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);
( g5 ~* V5 m+ J) v2 B    b[0]=1.0; d1=1.0*n; p=0.0; c=0.0;
    for (i=0; i&lt;=n-1; i++)
$ b9 E" k0 p" w' V+ E/ x2 S' G, i      { p=p+(x-z); c=c+y;}
    c=c/d1; p=p/d1;
    a[0]=c*b[0];
$ A# A7 h! }  G9 N4 Q  c    if (m&gt;1)
$ {- @7 a) C) F      { t[1]=1.0; t[0]=-p;
7 Z  B1 y1 S# e: E1 ~3 s        d2=0.0; c=0.0; g=0.0;
        for (i=0; i&lt;=n-1; i++)
$ A4 _7 _, C& n/ {1 e, @          { q=x-z-p; d2=d2+q*q;
/ M1 V) w; e, A9 b            c=c+y*q;
            g=g+(x-z)*q*q;
          }
        c=c/d2; p=g/d2; q=d2/d1;
        d1=d2;
        a[1]=c*t[1]; a[0]=c*t[0]+a[0];
      }
+ _' y9 G+ {  V( p    for (j=2; j&lt;=m-1; j++)
      { s[j]=t[j-1];
4 I; d: v6 `! t/ w        s[j-1]=-p*t[j-1]+t[j-2];
" n" X! e7 R- T' f( [, `1 n+ S        if (j&gt;=3)
* V/ g0 H9 g- c1 J7 v8 ?          for (k=j-2; k&gt;=1; k--)
            s[k]=-p*t[k]+t[k-1]-q*b[k];
" R2 E5 E, I, o7 d1 c( a9 k        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++)
          { q=s[j];
, s8 q$ w0 t; [            for (k=j-1; k&gt;=0; k--)
; Z# f# j( y/ ^              q=q*(x-z)+s[k];
            d2=d2+q*q; c=c+y*q;
            g=g+(x-z)*q*q;
          }
; A9 C7 u0 Y0 P" B7 c1 c        c=c/d2; p=g/d2; q=d2/d1;
$ o& X/ G: F6 r( W4 ]        d1=d2;
) k% Q6 L2 ^2 ]% A        a[j]=c*s[j]; t[j]=s[j];
8 Y: D1 ~9 K) N3 [2 M! y' l; F        for (k=j-1; k&gt;=0; k--)
) H" Y& T2 i8 t9 o0 s$ m          { a[k]=c*s[k]+a[k];
% q1 E1 h% z% p2 p3 Q( m            b[k]=t[k]; t[k]=s[k];
          }
. S# q) J9 t- J7 Y0 B3 B& ?, U      }
    dt[0]=0.0; dt[1]=0.0; dt[2]=0.0;
    for (i=0; i&lt;=n-1; i++)
  H7 w6 L% h8 d* \) \6 x* Q1 @      { q=a[m-1];
        for (k=m-2; k&gt;=0; k--)
          q=a[k]+q*(x-z);
( M8 j- ^4 i# o- @) V( T9 X/ c        p=q-y;
        if (fabs(p)&gt;dt[2]) dt[2]=fabs(p);
3 o* \4 w1 H9 S        dt[0]=dt[0]+p*p;
        dt[1]=dt[1]+fabs(p);
+ S5 _8 u6 K" V1 ~6 @/ a      }
    return;
# C4 q! I8 D/ l: ?  }</P><>龙贝格积分法</P><>#include "math.h"
* s( [. X- }* j  double fromb(a,b,eps)
  double a,b,eps;
; l) Q; Z+ X' T) ?& f8 C. @  { extern double frombf();
& ~, m7 p6 s' S    int m,n,i,k;
3 p2 |% a" U" c+ ^1 l2 n  ?4 \8 q    double y[10],h,ep,p,x,s,q;
; I/ e- u! _0 \4 h' {- N% n* ?1 G* A    h=b-a;
    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))
      { p=0.0;
* _  C* y2 N. k1 o        for (i=0;i&lt;=n-1;i++)
          { x=a+(i+0.5)*h;
            p=p+frombf(x);
          }
. m! H9 B" g( M1 t2 r        p=(y[0]+h*p)/2.0;
& g& u3 G; `1 c% W        s=1.0;
6 d1 `1 R" X! U        for (k=1;k&lt;=m;k++)
3 x/ z) l; M. ^$ V$ _  p" B          { s=4.0*s;
  T  p% w8 `2 Y/ c            q=(s*p-y[k-1])/(s-1.0);
+ h1 `9 D/ ^8 V& S% p8 K' Z            y[k-1]=p; p=q;
          }
        ep=fabs(q-y[m-1]);
3 F# l9 _. j& F. q5 _        m=m+1; y[m-1]=q; n=n+n; h=h/2.0;
8 i/ B& h- H- l: z; c% l1 X      }
1 u9 G4 q& |% q, c    return(q);
  }</P><>呵呵 希望对你有用！！</P>

西西

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

陈了

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