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

西西

<>本人紧急求助，希望各位大虾们伸出你们友情之手，帮我一下。</P>
8 B8 N- n! P9 \' y- t<>1。牛顿迭代法</P>
8 F3 u6 M- m6 S: F4 ?) n$ \<>1。全主消元法</P>
<>3。改进平方根法</P>
- f) J4 Q: e$ u! @! h<>4。牛顿向前，向后插值法</P>
<>5。加权最小二乘法</P>
. r- I' ?% g9 W& e<>6。龙贝格积分法</P>
- F# U1 b0 Z+ }* o<>拜托了。。。。</P>

loveyuan

<>1。牛顿迭代法</P><>
  #include "stdio.h"
, ^: i! s+ w- o$ z1 Y  #include "math.h"
9 f. o: N7 y$ m* Q  int dnewt(x,eps,js)
: u/ Y3 `" h' {. ?9 Q  int js;
' j4 P- {  C' x/ j2 z, ~  double *x,eps;
  { extern void dnewtf();
    int k,l;
+ |3 g  c! T) A/ b: P- A    double y[2],d,p,x0,x1;
    l=js; k=1; x0=*x;
( k+ @+ v/ ^) J2 b  T) C    dnewtf(x0,y);
2 s* [; ?) c) E* p- _    d=eps+1.0;
/ A$ {. f- j. @- P4 ^% W8 L2 Z    while ((d&gt;=eps)&amp;&amp;(l!=0))
      { if (fabs(y[1])+1.0==1.0)
          { printf("err\n"); return(-1);}
        x1=x0-y[0]/y[1];
        dnewtf(x1,y);
1 Q+ N7 |3 k- l: n* z# e2 ?3 Y% y        d=fabs(x1-x0); p=fabs(y[0]);
1 n/ y& u7 x$ Y- l) Q5 Y6 P        if (p&gt;d) d=p;
; K0 a3 x! z( O2 D- q2 |( P$ n        x0=x1; l=l-1;
3 J4 z2 K0 l7 Q' j/ Q      }
    *x=x1;
    k=js-l;
" {. N  D+ G. I+ y3 W  u    return(k);
1 r6 |* j" g* c. _. Y* h( A& a  }</P><>全主消元法</P><>#include "stdlib.h"
; K% ]$ C3 `( L9 Q5 h/ e  #include "stdio.h"
! t* r/ }' q# c/ v( N  int acgas(ar,ai,n,br,bi)
  int n;
( i9 Y: @7 `0 s: c  double ar[],ai[],br[],bi[];
* |: b- d7 |# v1 B& L% H( {  { int *js,l,k,i,j,is,u,v;
- W9 }$ J- [% O# W# ]9 b* w/ K    double p,q,s,d;
    js=malloc(n*sizeof(int));
    for (k=0;k&lt;=n-2;k++)
) F/ W2 U, s  L, J, x9 j      { d=0.0;
        for (i=k;i&lt;=n-1;i++)
$ C- S+ ?; z" r5 q        for (j=k;j&lt;=n-1;j++)
          { u=i*n+j;
            p=ar*ar+ai*ai;
5 Q) U- J$ k% l4 ?7 ~            if (p&gt;d) {d=p;js[k]=j;is=i;}
          }
+ U9 x# r# P( F* i' P        if (d+1.0==1.0)
  z+ |- T8 @8 h9 R          { free(js); printf("err**fail\n");
, Y5 M' a6 q4 E. T  z; ^            return(0);
          }
        if (is!=k)
          { for (j=k;j&lt;=n-1;j++)
1 L2 n, e( u, X: i; _% p( @$ T* z$ Z              { u=k*n+j; v=is*n+j;
                p=ar; ar=ar[v]; ar[v]=p;
                p=ai; ai=ai[v]; ai[v]=p;
! G5 M1 c' G) F4 |# q% A2 u3 y              }
: u, b; m+ P( i5 c6 C0 Q            p=br[k]; br[k]=br[is]; br[is]=p;
/ n; e! i" E: n, p. B9 W! M            p=bi[k]; bi[k]=bi[is]; bi[is]=p;
          }
        if (js[k]!=k)
( t9 ~1 |( Y% Q# p' {          for (i=0;i&lt;=n-1;i++)
9 `9 u2 u0 V' h/ P  ]" h            { 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;
            }
        v=k*n+k;
$ c% k! I- B0 y        for (j=k+1;j&lt;=n-1;j++)
1 I- n" A8 p3 b          { u=k*n+j;
" R( X. U# I0 o6 E0 q) ]6 q. {; u            p=ar*ar[v]; q=-ai*ai[v];
            s=(ar[v]-ai[v])*(ar+ai);
            ar=(p-q)/d; ai=(s-p-q)/d;
3 [: ^; s: |* p" S% x) Z  j          }
        p=br[k]*ar[v]; q=-bi[k]*ai[v];
9 b3 O1 S& A/ O( u1 L! i        s=(ar[v]-ai[v])*(br[k]+bi[k]);
" G3 [, r7 W1 {- j1 c        br[k]=(p-q)/d; bi[k]=(s-p-q)/d;
3 c+ q& e. f, y9 |  B0 j        for (i=k+1;i&lt;=n-1;i++)
% \" @  f( L  `' o0 ?0 f  M          { u=i*n+k;
* U  B: ?( y# n            for (j=k+1;j&lt;=n-1;j++)
              { v=k*n+j; l=i*n+j;
4 z5 f" {: ]5 Q' g' q2 P: a                p=ar*ar[v]; q=ai*ai[v];
) d5 b) o5 g# ~9 B7 C                s=(ar+ai)*(ar[v]+ai[v]);
8 |' O% t( J$ e- E* V                ar[l]=ar[l]-p+q;
* m. t0 [- L9 R& J% e8 u, F, H                ai[l]=ai[l]-s+p+q;
$ L3 l) \6 B5 ], u/ f' l              }
& c# b& {& Z8 C0 C- X            p=ar*br[k]; q=ai*bi[k];
            s=(ar+ai)*(br[k]+bi[k]);
            br=br-p+q; bi=bi-s+p+q;
          }
      }
    u=(n-1)*n+n-1;
0 ~' f* L) }8 e  K8 r( k; R    d=ar*ar+ai*ai;
    if (d+1.0==1.0)
      { free(js); printf("err**fail\n");
        return(0);
      }
, l' q& n! ]% n9 {! y    p=ar*br[n-1]; q=-ai*bi[n-1];
    s=(ar-ai)*(br[n-1]+bi[n-1]);
    br[n-1]=(p-q)/d; bi[n-1]=(s-p-q)/d;
* W7 o" w" w' V6 g4 A' f( @/ S' l    for (i=n-2;i&gt;=0;i--)
( S0 U7 [$ u: X2 @1 Q    for (j=i+1;j&lt;=n-1;j++)
! Y1 ?7 M0 a: ]( i6 \      { u=i*n+j;
8 M/ T8 M4 Z0 z, V; |        p=ar*br[j]; q=ai*bi[j];
# s" T8 [# P: y# k# d        s=(ar+ai)*(br[j]+bi[j]);
        br=br-p+q;
; I- D) p; a5 u+ ?! J) _        bi=bi-s+p+q;
      }
% d$ @+ r8 L1 q    js[n-1]=n-1;
& H& f" n* f" E2 K    for (k=n-1;k&gt;=0;k--)
# J0 y  r5 z- \' A8 n: |8 r6 V      if (js[k]!=k)
7 k# `1 R( F  q$ H. _" e: B8 M$ T! I: d        { p=br[k]; br[k]=br[js[k]]; br[js[k]]=p;
          p=bi[k]; bi[k]=bi[js[k]]; bi[js[k]]=p;
        }
; ?* g8 Q3 r% d# @    free(js);
# D# j1 F7 y7 [/ O/ z% F0 T    return(1);
% g: [7 U2 Z2 h  G5 q  Q  }</P><>平方根法</P><>#include "math.h"
" |4 N8 e4 K" M! w) H, a' @  #include "stdio.h"
$ X6 T+ F' W1 x7 r* w  int achol(a,n,m,d)
0 t7 ^- T# x7 P4 o$ T* j2 h  int n,m;
  double a[],d[];
9 W: L% X  i. P( a  { int i,j,k,u,v;
    if ((a[0]+1.0==1.0)||(a[0]&lt;0.0))
      { printf("fail\n"); return(-2);}
/ k3 ?4 ~' o4 Z- b+ ^8 ?) N, }: K    a[0]=sqrt(a[0]);
7 t5 h) ]2 w: X# }    for (j=1; j&lt;=n-1; j++) a[j]=a[j]/a[0];
1 h# s5 P# i9 D0 F! u    for (i=1; i&lt;=n-1; i++)
' @% [8 ^' U2 b; C" P+ H, t      { u=i*n+i;
: e* t0 W5 o" O! R        for (j=1; j&lt;=i; j++)
5 ^! U, A& f5 n- V0 R          { v=(j-1)*n+i;
/ r+ P, W# {5 n            a=a-a[v]*a[v];
          }
        if ((a+1.0==1.0)||(a&lt;0.0))
          { printf("fail\n"); return(-2);}
        a=sqrt(a);
        if (i!=(n-1))
3 L* V) l# S7 o$ X4 w! Y          { for (j=i+1; j&lt;=n-1; j++)
              { v=i*n+j;
                for (k=1; k&lt;=i; k++)
  t" e2 v7 b- x1 x, |/ a                  a[v]=a[v]-a[(k-1)*n+i]*a[(k-1)*n+j];
* Z' ?! n+ G- d; c, e* T                a[v]=a[v]/a;
2 W, k( [4 C9 }5 ^/ M' W4 {              }
8 H* W. b+ W3 l( j6 U/ x          }
      }
    for (j=0; j&lt;=m-1; j++)
      { 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++)
- C& @" o3 z2 v& E0 b, f9 J9 V              d[v]=d[v]-a[(k-1)*n+i]*d[(k-1)*m+j];
' l: ?: h, y9 l2 P, u0 h            d[v]=d[v]/a;
          }
      }
8 }, m( B9 Z. r. _    for (j=0; j&lt;=m-1; j++)
      { u=(n-1)*m+j;
' w- m+ _  @) d1 {9 b& W        d=d/a[n*n-1];
        for (k=n-1; k&gt;=1; k--)
          { u=(k-1)*m+j;
) r# n0 U$ v& G8 g4 s            for (i=k; i&lt;=n-1; i++)
  C: s) f5 D  F& U" u3 t7 W0 M1 D9 N              { v=(k-1)*n+i;
9 r( n, F: z8 i$ a- F+ p9 t7 y                d=d-a[v]*d[i*m+j];
              }
+ q. n- X* x) |% Y" q2 Y/ |            v=(k-1)*n+k-1;
) y2 g1 @5 E; _- B  W& L  V            d=d/a[v];
          }
, n8 N1 M. p9 w! x. t+ c' y+ T      }
% w7 J/ ]) ]% T! A+ A    return(2);
2 Y( S1 ~# B5 O6 D  }</P><>牛顿向前，向后插值法</P><>double eelgr(x0,h,n,y,t)
  int n;
' F; H4 Z) x( l  double x0,h,t,y[];
  s; Q+ q/ m4 a1 W  { int i,j,k,m;
    double z,s,xi,xj;
    float p,q;
7 E. W4 Y2 |* G4 {$ U* N    z=0.0;
& x9 E+ ?3 l9 X8 E$ E; u. w    if (n&lt;1) return(z);
& U3 ^, M; J6 U( w8 v3 _( k, ]    if (n==1) { z=y[0]; return(z);}
; @( M$ Z2 _! ^    if (n==2)
/ Q& z; H2 |+ u      { z=(y[1]*(t-x0)-y[0]*(t-x0-h))/h;
4 b& y6 x% I3 |2 `        return(z);
      }
    if (t&gt;x0)
      { p=(t-x0)/h; i=(int)p; q=(float)i;
        if (p&gt;q) i=i+1;
      }
    else i=0;
    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++)
      { s=1.0; xi=x0+i*h;
* x" x$ N# r" t# {7 b5 Q        for (j=k; j&lt;=m; j++)
          if (j!=i)
$ [( f4 v# d- U% P' A            { xj=x0+j*h;
              s=s*(t-xj)/(xi-xj);
            }
( G' `+ G! t1 ]% E- N        z=z+s*y;
: ?1 I" X* u1 A. m8 Q      }
; ?2 W$ T% ^9 [" j    return(z);
  }
) s3 u! p7 q: D6 n0 c+ b8 T6 r向前，向后是一样的思想！！</P><>加权最小二乘法</P><>#include "math.h"
  void hpir1(x,y,n,a,m,dt)
% G% ~0 ~; p; @  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];
8 @: q, V! ?  B8 |0 u  D    for (i=0; i&lt;=m-1; i++) a=0.0;
    if (m&gt;n) m=n;
    if (m&gt;20) m=20;
! _3 O( w3 I  \7 a- p    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;
$ d8 s* C, F& Z7 V$ N) }( V    for (i=0; i&lt;=n-1; i++)
      { p=p+(x-z); c=c+y;}
1 r) |2 Z9 y, I2 t# G    c=c/d1; p=p/d1;
    a[0]=c*b[0];
! _; @* N/ t+ K9 T! K    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++)
          { q=x-z-p; d2=d2+q*q;
            c=c+y*q;
! o7 ~0 `8 V8 Q% w/ Y9 I            g=g+(x-z)*q*q;
          }
        c=c/d2; p=g/d2; q=d2/d1;
+ [* M8 I0 b9 Z4 A' q" u        d1=d2;
        a[1]=c*t[1]; a[0]=c*t[0]+a[0];
      }
* r1 z4 H3 r1 n7 U5 e    for (j=2; j&lt;=m-1; j++)
      { s[j]=t[j-1];
/ t3 m# R' ?( M        s[j-1]=-p*t[j-1]+t[j-2];
2 Z+ W; c( t# N. G; b. S        if (j&gt;=3)
          for (k=j-2; k&gt;=1; k--)
            s[k]=-p*t[k]+t[k-1]-q*b[k];
1 O: J% B( _3 [& j- [        s[0]=-p*t[0]-q*b[0];
0 e  m. }2 Z. W5 s$ b% T        d2=0.0; c=0.0; g=0.0;
        for (i=0; i&lt;=n-1; i++)
          { q=s[j];
( k& F* [& ^& g  c            for (k=j-1; k&gt;=0; k--)
3 O, c9 \/ M8 N; |2 g              q=q*(x-z)+s[k];
            d2=d2+q*q; c=c+y*q;
2 S* F2 m1 e  @            g=g+(x-z)*q*q;
; c1 `2 ]3 {6 Z  u0 n          }
        c=c/d2; p=g/d2; q=d2/d1;
8 n+ o0 E0 f/ Y        d1=d2;
6 H0 z  j! r9 x        a[j]=c*s[j]; t[j]=s[j];
! v6 _7 |$ P. ^9 M2 F        for (k=j-1; k&gt;=0; k--)
          { a[k]=c*s[k]+a[k];
; S$ {" ~. K7 {. u0 r; R            b[k]=t[k]; t[k]=s[k];
          }
      }
    dt[0]=0.0; dt[1]=0.0; dt[2]=0.0;
    for (i=0; i&lt;=n-1; i++)
      { q=a[m-1];
. _+ d% T, q- L7 V        for (k=m-2; k&gt;=0; k--)
; G* j5 y2 j) U' @          q=a[k]+q*(x-z);
7 d! C1 F$ ^* u/ a        p=q-y;
' p1 ^  l9 e7 \5 E+ ]        if (fabs(p)&gt;dt[2]) dt[2]=fabs(p);
. f# \4 ^& t& }; w        dt[0]=dt[0]+p*p;
        dt[1]=dt[1]+fabs(p);
      }
    return;
  }</P><>龙贝格积分法</P><>#include "math.h"
# {. v1 ], Q5 P7 R) R  double fromb(a,b,eps)
+ Y* t( {& X) O' I" X8 V4 C  double a,b,eps;
8 _9 A; B; z# ]  b# I" o4 \' m  { extern double frombf();
# S' v/ T: F6 o! `    int m,n,i,k;
    double y[10],h,ep,p,x,s,q;
8 X" k# v5 z  |9 l    h=b-a;
) d8 n: v; h8 y9 m+ t    y[0]=h*(frombf(a)+frombf(b))/2.0;
    m=1; n=1; ep=eps+1.0;
, @% G5 j. ]: w; @    while ((ep&gt;=eps)&amp;&amp;(m&lt;=9))
      { p=0.0;
  M! W! C( K( T5 c4 N. k# w' h. d        for (i=0;i&lt;=n-1;i++)
          { x=a+(i+0.5)*h;
            p=p+frombf(x);
          }
        p=(y[0]+h*p)/2.0;
        s=1.0;
0 X3 G  d4 r; x. S) j' v        for (k=1;k&lt;=m;k++)
          { s=4.0*s;
            q=(s*p-y[k-1])/(s-1.0);
# b" E- N8 p% Q7 m, ]# d            y[k-1]=p; p=q;
+ m! m' {; Y. N0 c1 G( a# i          }
        ep=fabs(q-y[m-1]);
        m=m+1; y[m-1]=q; n=n+n; h=h/2.0;
      }
    return(q);
6 ~# J0 e) a0 g( h" b  }</P><>呵呵 希望对你有用！！</P>

西西

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

陈了

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