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

西西

<>本人紧急求助，希望各位大虾们伸出你们友情之手，帮我一下。</P>
<>1。牛顿迭代法</P>
! l" Q( t1 Z5 g7 h7 |0 S<>1。全主消元法</P>
<>3。改进平方根法</P>
<>4。牛顿向前，向后插值法</P>
% @$ @" B9 N7 f' ?) A# K<>5。加权最小二乘法</P>
  ^6 E2 C! z' d% K$ \! u1 S<>6。龙贝格积分法</P>
3 ~9 r- R3 S2 r# X0 J" y<>拜托了。。。。</P>

loveyuan

<>1。牛顿迭代法</P><>
* u* K- e. _# P3 z+ o9 y  #include "stdio.h"
& ?; f0 P$ e  q+ l% s  #include "math.h"
  int dnewt(x,eps,js)
  int js;
  double *x,eps;
) h) y; m! l/ J/ Y) A( \  { extern void dnewtf();
2 }7 f# [$ ]/ e- }/ }) j, V    int k,l;
6 G, Q# s0 O+ m8 {  W3 h6 n- ^    double y[2],d,p,x0,x1;
    l=js; k=1; x0=*x;
+ [9 X; ?6 c  x! m/ d    dnewtf(x0,y);
    d=eps+1.0;
    while ((d&gt;=eps)&amp;&amp;(l!=0))
' G5 B- k3 p( ~- x      { if (fabs(y[1])+1.0==1.0)
          { printf("err\n"); return(-1);}
0 n$ Z* ~0 L' a* W5 A) j/ c6 D        x1=x0-y[0]/y[1];
( V/ n+ c; ^  M2 \0 T        dnewtf(x1,y);
        d=fabs(x1-x0); p=fabs(y[0]);
5 [- c+ ~, t( L/ z: _+ o        if (p&gt;d) d=p;
        x0=x1; l=l-1;
      }
    *x=x1;
( ~  _: B6 L6 E    k=js-l;
    return(k);
6 S1 D3 w- y4 [9 G: u' G; s- O8 k  }</P><>全主消元法</P><>#include "stdlib.h"
9 _) h% W" U. v3 ^  #include "stdio.h"
4 s' T" V  |: Y' T$ O: x) U- M  int acgas(ar,ai,n,br,bi)
  int n;
3 p7 S. E4 ^, \, k5 \6 B; ]  double ar[],ai[],br[],bi[];
  { int *js,l,k,i,j,is,u,v;
    double p,q,s,d;
    js=malloc(n*sizeof(int));
    for (k=0;k&lt;=n-2;k++)
      { d=0.0;
        for (i=k;i&lt;=n-1;i++)
        for (j=k;j&lt;=n-1;j++)
, r8 C9 o& J1 x, @! w- v          { u=i*n+j;
; w3 H9 s+ L% Z8 D/ h; y* 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");
: ]; a. `- C- e$ D. X            return(0);
/ C" l- R3 E+ ?- S- B7 u# j          }
        if (is!=k)
8 b9 C0 j4 H4 [2 |  \, Z8 G          { 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;
+ U: G1 g, h% g( Y: F! c* h              }
            p=br[k]; br[k]=br[is]; br[is]=p;
            p=bi[k]; bi[k]=bi[is]; bi[is]=p;
          }
        if (js[k]!=k)
7 t) i9 x1 j# |5 r6 h          for (i=0;i&lt;=n-1;i++)
9 @( m* I( }" w            { u=i*n+k; v=i*n+js[k];
6 A% C9 c7 V6 ?              p=ar; ar=ar[v]; ar[v]=p;
  A4 F5 t: R8 P) d* q              p=ai; ai=ai[v]; ai[v]=p;
            }
+ X3 ]* U0 w# p' ~5 y2 L        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];
            s=(ar[v]-ai[v])*(ar+ai);
- M% P5 |; C( z, ]7 j  y            ar=(p-q)/d; ai=(s-p-q)/d;
0 h# w2 ^+ i% Q$ F* {6 s% O          }
        p=br[k]*ar[v]; q=-bi[k]*ai[v];
        s=(ar[v]-ai[v])*(br[k]+bi[k]);
: a% e/ T' L# `- j. L9 y        br[k]=(p-q)/d; bi[k]=(s-p-q)/d;
        for (i=k+1;i&lt;=n-1;i++)
          { u=i*n+k;
            for (j=k+1;j&lt;=n-1;j++)
+ }+ `1 _  \) l5 G  U3 r  g$ O              { v=k*n+j; l=i*n+j;
3 {4 [9 T1 c, c+ g: i                p=ar*ar[v]; q=ai*ai[v];
                s=(ar+ai)*(ar[v]+ai[v]);
" c1 _" C: O4 _                ar[l]=ar[l]-p+q;
( y3 {$ H6 B! I  \3 a, h2 m                ai[l]=ai[l]-s+p+q;
              }
, q! [" F% ]! u6 r. @! t. a            p=ar*br[k]; q=ai*bi[k];
$ @5 t" }3 v# ?' A2 V% o            s=(ar+ai)*(br[k]+bi[k]);
; j; T+ ^# i0 ?1 s0 M3 a2 H' f$ M9 Z3 T4 g            br=br-p+q; bi=bi-s+p+q;
          }
# L. M8 |, m) e( ^6 q% T; H      }
    u=(n-1)*n+n-1;
+ S  E& [) h+ N' G6 e4 K0 i" `    d=ar*ar+ai*ai;
+ t% u( N- H& D) S! C    if (d+1.0==1.0)
( \8 g  b6 L' W; C' l) C      { free(js); printf("err**fail\n");
/ d" @+ ?; q" {( N8 j" k        return(0);
      }
6 Z0 T9 B2 {8 |+ x7 Q0 j    p=ar*br[n-1]; q=-ai*bi[n-1];
    s=(ar-ai)*(br[n-1]+bi[n-1]);
; `- G9 ]$ g; {. r. c* f) A    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++)
      { u=i*n+j;
3 t  e% d4 y6 R' `! K5 u# ?        p=ar*br[j]; q=ai*bi[j];
        s=(ar+ai)*(br[j]+bi[j]);
        br=br-p+q;
        bi=bi-s+p+q;
1 S  y, i( @6 }5 ~, b8 m; t      }
2 w% V1 r0 |8 C6 x- f    js[n-1]=n-1;
    for (k=n-1;k&gt;=0;k--)
      if (js[k]!=k)
6 |, u# n& S( e$ F8 w        { p=br[k]; br[k]=br[js[k]]; br[js[k]]=p;
9 C! L8 C: p9 {5 i. X          p=bi[k]; bi[k]=bi[js[k]]; bi[js[k]]=p;
" o3 a# W" Q# _9 F9 L        }
    free(js);
    return(1);
  }</P><>平方根法</P><>#include "math.h"
  #include "stdio.h"
  int achol(a,n,m,d)
0 s; T2 T4 h+ ?7 U0 N1 _  E  int n,m;
  I" a$ x4 x8 I+ Y+ R0 ~7 G9 q  double a[],d[];
, b& U6 _5 _1 W7 B. b7 N! T; t( a  { 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]);
    for (j=1; j&lt;=n-1; j++) a[j]=a[j]/a[0];
- M7 e# D, |" {    for (i=1; i&lt;=n-1; i++)
# k$ U8 a) z  w) o) u* Q      { u=i*n+i;
        for (j=1; j&lt;=i; j++)
          { v=(j-1)*n+i;
( O: ^4 i: h2 s4 O  o5 O# z1 W3 R            a=a-a[v]*a[v];
8 Z; a; s" g) {( Y+ F9 @& Y          }
        if ((a+1.0==1.0)||(a&lt;0.0))
          { printf("fail\n"); return(-2);}
        a=sqrt(a);
0 ]0 m: \. X0 V# y        if (i!=(n-1))
          { for (j=i+1; j&lt;=n-1; j++)
8 _' a! D8 E% D/ u              { 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;
              }
0 G6 _& D) x6 Z% N( P          }
& B9 s3 j0 N+ ^3 s. o2 B      }
$ k/ L" u9 N3 {) f/ w    for (j=0; j&lt;=m-1; j++)
      { d[j]=d[j]/a[0];
! e* M+ {3 X) v+ e" u        for (i=1; i&lt;=n-1; i++)
9 }6 R8 r" m1 b' R          { 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];
  c8 a# }3 w. E! i. n1 }  M            d[v]=d[v]/a;
          }
0 N% }3 c5 {0 u$ q4 \      }
    for (j=0; j&lt;=m-1; j++)
      { u=(n-1)*m+j;
- y% ^0 r+ o" p8 m        d=d/a[n*n-1];
6 a( B6 @/ v+ [% S& ^% o- H2 A+ m8 T        for (k=n-1; k&gt;=1; k--)
. B; T% @, [8 y& ~* x          { u=(k-1)*m+j;
# i: N: Y- @# Q8 j" R            for (i=k; i&lt;=n-1; i++)
              { v=(k-1)*n+i;
                d=d-a[v]*d[i*m+j];
4 F/ O4 v5 Y( u: W- v* v. O              }
            v=(k-1)*n+k-1;
, t* R* H+ K' h. `4 P            d=d/a[v];
/ ~7 f3 n- x- E) C8 T          }
7 {; }( ]3 G6 ]# W      }
    return(2);
" g3 u1 y8 d, i0 ~" g* g/ Y  }</P><>牛顿向前，向后插值法</P><>double eelgr(x0,h,n,y,t)
4 T/ ]- g" m6 P* [# y* p& y$ u, Y  int n;
  double x0,h,t,y[];
  { int i,j,k,m;
" x+ ]& @0 B/ F/ j    double z,s,xi,xj;
1 {8 g. O9 R3 J4 F8 l2 F    float p,q;
    z=0.0;
' f. m2 f. c% }, P# G    if (n&lt;1) return(z);
, j- y( f" f! g9 N* y8 p5 u! y  a    if (n==1) { z=y[0]; return(z);}
    if (n==2)
8 ?8 u5 _' O% g3 c7 a; N      { z=(y[1]*(t-x0)-y[0]*(t-x0-h))/h;
6 C/ k( h) y* Q+ G( J5 ~3 t* |        return(z);
& {+ x8 P* V. ^) ]      }
, E( _- e" l& M2 a4 |$ t    if (t&gt;x0)
) g7 N, d; y+ w" R: X3 y      { p=(t-x0)/h; i=(int)p; q=(float)i;
        if (p&gt;q) i=i+1;
      }
& L8 J& L. V) V( K. P3 s    else i=0;
    k=i-4;
  G6 B7 ~  L: W$ E; W  X    if (k&lt;0) k=0;
    m=i+3;
    if (m&gt;n-1) m=n-1;
2 }& V5 _+ r, O    for (i=k;i&lt;=m;i++)
6 l# |3 \: L& o3 V      { s=1.0; xi=x0+i*h;
        for (j=k; j&lt;=m; j++)
7 J3 D0 N" a1 I+ C+ p( j  E" R          if (j!=i)
- d( {- E$ R5 F8 l5 Q( v            { xj=x0+j*h;
              s=s*(t-xj)/(xi-xj);
            }
        z=z+s*y;
      }
% @( S/ m- @9 l* A    return(z);
  M% Q' Q8 K& N# S: n: D3 n$ l  }
向前，向后是一样的思想！！</P><>加权最小二乘法</P><>#include "math.h"
3 J5 W" g  d- R4 U7 d  G  void hpir1(x,y,n,a,m,dt)
" L0 ]3 L# U9 N  int n,m;
4 F- y/ X  Z6 M  double x[],y[],a[],dt[];
  { int i,j,k;
4 V+ y, h+ A4 |, y& T6 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;
9 M3 Z7 z3 \/ u3 N+ f1 n    if (m&gt;20) m=20;
    z=0.0;
* O$ \! B$ R. m( [, o; m) [    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++)
" k* y- C! T  h' ~" ^      { p=p+(x-z); c=c+y;}
    c=c/d1; p=p/d1;
    a[0]=c*b[0];
    if (m&gt;1)
0 X$ H6 b7 H, i9 }( G& y  W      { t[1]=1.0; t[0]=-p;
; v* \0 n4 i3 {0 J        d2=0.0; c=0.0; g=0.0;
        for (i=0; i&lt;=n-1; i++)
! F8 Q7 D) R5 }7 `+ ~- e          { q=x-z-p; d2=d2+q*q;
            c=c+y*q;
            g=g+(x-z)*q*q;
$ K; O. d& a( L7 K          }
        c=c/d2; p=g/d2; q=d2/d1;
        d1=d2;
        a[1]=c*t[1]; a[0]=c*t[0]+a[0];
      }
) K. @5 O# Z2 T4 y    for (j=2; j&lt;=m-1; j++)
      { s[j]=t[j-1];
# |+ w- j/ d, g        s[j-1]=-p*t[j-1]+t[j-2];
        if (j&gt;=3)
% z8 \  P6 d. {3 P* w6 Z( r          for (k=j-2; k&gt;=1; k--)
) Z9 _& v9 f" L: K1 O, Z            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;
        for (i=0; i&lt;=n-1; i++)
0 R: X% \/ Q4 x9 M( M1 {          { q=s[j];
  S6 w3 G# Z3 D' q- `            for (k=j-1; k&gt;=0; k--)
  z4 O% b$ O0 I* ?3 g              q=q*(x-z)+s[k];
& L3 o; H9 L  b3 ^2 Q8 S            d2=d2+q*q; c=c+y*q;
            g=g+(x-z)*q*q;
: \+ r1 ?0 n; b+ I5 {1 ]          }
4 J$ F3 a$ s- N! T# @        c=c/d2; p=g/d2; q=d2/d1;
        d1=d2;
: Z' `4 A% I" F5 {* D0 q" m        a[j]=c*s[j]; t[j]=s[j];
        for (k=j-1; k&gt;=0; k--)
          { a[k]=c*s[k]+a[k];
3 Y" K8 p0 n7 z2 p, \3 c1 W            b[k]=t[k]; t[k]=s[k];
2 R- K% y# y# y8 l' M" Z          }
      }
    dt[0]=0.0; dt[1]=0.0; dt[2]=0.0;
5 R+ d& a+ O: g6 D1 j    for (i=0; i&lt;=n-1; i++)
) w* G& G, q& w7 Y      { q=a[m-1];
        for (k=m-2; k&gt;=0; k--)
          q=a[k]+q*(x-z);
: l# z6 y% r: h2 h        p=q-y;
        if (fabs(p)&gt;dt[2]) dt[2]=fabs(p);
        dt[0]=dt[0]+p*p;
& G4 V9 {7 e1 K+ l; A; X/ d        dt[1]=dt[1]+fabs(p);
' v1 h. @( P5 }/ L1 `& Q8 _5 B      }
, E# F& u- d% U$ c2 L# S: `+ ]    return;
) o. U6 A' f/ J: V  }</P><>龙贝格积分法</P><>#include "math.h"
  double fromb(a,b,eps)
  double a,b,eps;
% H' s2 j9 {# J3 a0 a  { extern double frombf();
7 J: ]2 M" R1 H* ^& U    int m,n,i,k;
' c; S6 |* R6 z    double y[10],h,ep,p,x,s,q;
    h=b-a;
    y[0]=h*(frombf(a)+frombf(b))/2.0;
! z7 I5 P, J6 W: J4 R  h9 Z' z! M+ b    m=1; n=1; ep=eps+1.0;
    while ((ep&gt;=eps)&amp;&amp;(m&lt;=9))
" S% G" D4 P" I5 x4 R+ [+ @0 `$ B      { p=0.0;
, ]) i& M, J6 \- ]7 ]& G# z        for (i=0;i&lt;=n-1;i++)
          { x=a+(i+0.5)*h;
2 z% c* r1 b. a! Y& L$ }            p=p+frombf(x);
# ^- a0 h4 {; u. n' \( J% U* c6 ]          }
        p=(y[0]+h*p)/2.0;
        s=1.0;
# u- r. i6 M( l4 ^4 A2 i        for (k=1;k&lt;=m;k++)
& ?* g( q4 a; D0 c# q; ~) R          { s=4.0*s;
9 ^5 s( I3 d# j) F. g. {. \            q=(s*p-y[k-1])/(s-1.0);
( r# B4 D7 e+ Y/ x            y[k-1]=p; p=q;
* i- \. k9 F, g1 M          }
        ep=fabs(q-y[m-1]);
3 U2 p0 \0 t, ?. _! _5 N) L        m=m+1; y[m-1]=q; n=n+n; h=h/2.0;
; c7 w0 Y3 D$ a      }
6 a, T- [3 J% x4 u    return(q);
) {( L& r& P0 G: b  \  }</P><>呵呵 希望对你有用！！</P>

西西

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

陈了

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