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

西西

<>本人紧急求助，希望各位大虾们伸出你们友情之手，帮我一下。</P>
<>1。牛顿迭代法</P>
- H' ~: }% M8 G1 D+ |2 I<>1。全主消元法</P>
% J/ ]! t# x7 }1 Y<>3。改进平方根法</P>
<>4。牛顿向前，向后插值法</P>
<>5。加权最小二乘法</P>
<>6。龙贝格积分法</P>
<>拜托了。。。。</P>

loveyuan

<>1。牛顿迭代法</P><>
  #include "stdio.h"
# o8 L7 R" }3 x7 q  #include "math.h"
  int dnewt(x,eps,js)
- o1 |, l9 Q/ v4 S  int js;
  double *x,eps;
) H) U1 b* F% |$ F: Y) b  { extern void dnewtf();
6 h; _& b0 I+ w" ?& {    int k,l;
    double y[2],d,p,x0,x1;
    l=js; k=1; x0=*x;
! l8 @1 ]$ a! n" f. E( r7 n    dnewtf(x0,y);
  l- i% W7 W0 W' X    d=eps+1.0;
( `  V8 I, m" L% A4 _8 w    while ((d&gt;=eps)&amp;&amp;(l!=0))
- v! t$ d) n& a7 U% ]$ d      { if (fabs(y[1])+1.0==1.0)
          { printf("err\n"); return(-1);}
        x1=x0-y[0]/y[1];
        dnewtf(x1,y);
+ P$ L8 U! v* g  z2 d" [        d=fabs(x1-x0); p=fabs(y[0]);
        if (p&gt;d) d=p;
        x0=x1; l=l-1;
      }
    *x=x1;
    k=js-l;
, S2 ~  A# \# U) H4 u# B9 m    return(k);
  }</P><>全主消元法</P><>#include "stdlib.h"
6 m6 p2 ^5 l- [4 X# X. a* I! Q8 @  #include "stdio.h"
  int acgas(ar,ai,n,br,bi)
  int n;
  double ar[],ai[],br[],bi[];
  { int *js,l,k,i,j,is,u,v;
' e& F  o" m4 w    double p,q,s,d;
    js=malloc(n*sizeof(int));
    for (k=0;k&lt;=n-2;k++)
6 [/ `) ?2 }9 U, Y      { d=0.0;
        for (i=k;i&lt;=n-1;i++)
- P0 K9 u' f' s3 u3 X) F$ Y* C  [* X        for (j=k;j&lt;=n-1;j++)
1 c) R, Y) j* `8 M$ ]+ w          { u=i*n+j;
. D, Z' U1 n; i& l            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");
+ E. X" |, W: O  S            return(0);
          }
/ C* o6 J1 m$ i) ~# B  x2 R7 u3 U        if (is!=k)
          { for (j=k;j&lt;=n-1;j++)
- y) O, u5 {" J; `" ?              { u=k*n+j; v=is*n+j;
! K8 `8 _( E* y" W                p=ar; ar=ar[v]; ar[v]=p;
2 y6 O2 Q- K9 b4 _$ D                p=ai; ai=ai[v]; ai[v]=p;
: z  r  y: G( b6 t, R              }
            p=br[k]; br[k]=br[is]; br[is]=p;
8 ]0 v4 X1 M0 ?2 C/ k* Y# c            p=bi[k]; bi[k]=bi[is]; bi[is]=p;
5 k8 W& |! ~5 i, x9 u8 ^8 W          }
5 c2 s- N, l/ e! b        if (js[k]!=k)
          for (i=0;i&lt;=n-1;i++)
5 R0 E* d! P( v            { u=i*n+k; v=i*n+js[k];
8 d% K/ E, U) A8 w# ?5 B              p=ar; ar=ar[v]; ar[v]=p;
              p=ai; ai=ai[v]; ai[v]=p;
' y" B* h" [6 b5 O            }
        v=k*n+k;
' e- ^. v+ x. m6 r* B        for (j=k+1;j&lt;=n-1;j++)
' p* p' F+ {2 T7 m% L* g& O          { u=k*n+j;
            p=ar*ar[v]; q=-ai*ai[v];
1 a6 m3 Z. g$ {# ?3 M            s=(ar[v]-ai[v])*(ar+ai);
$ p& k1 S0 n) c& h! s( O- P+ }            ar=(p-q)/d; ai=(s-p-q)/d;
- L1 H4 u6 a  j+ |- L* J8 x          }
9 l+ O5 O" E, D* E! P5 i& l        p=br[k]*ar[v]; q=-bi[k]*ai[v];
/ F  A1 q. h$ P' }        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;
            for (j=k+1;j&lt;=n-1;j++)
! Q% u# C/ V1 p, I& M1 B0 j& u              { v=k*n+j; l=i*n+j;
$ D0 u  [; C9 @$ {                p=ar*ar[v]; q=ai*ai[v];
                s=(ar+ai)*(ar[v]+ai[v]);
                ar[l]=ar[l]-p+q;
; m2 J; [( E' O3 W                ai[l]=ai[l]-s+p+q;
2 D9 D1 r7 B' x% B/ a$ q! c: O              }
3 ^: N% t& ~: m* L1 l            p=ar*br[k]; q=ai*bi[k];
            s=(ar+ai)*(br[k]+bi[k]);
1 m4 _" i+ u% H* Y+ w- _' @            br=br-p+q; bi=bi-s+p+q;
: e8 S  f( c3 c! j) i0 @0 b          }
) o6 K( ~5 g1 w' q. q      }
6 V" f, ]) t6 M" ?( T    u=(n-1)*n+n-1;
5 h. H6 r8 I4 `3 M+ g    d=ar*ar+ai*ai;
* V  S4 b0 T; b- C    if (d+1.0==1.0)
$ T0 |7 [. |* S      { free(js); printf("err**fail\n");
        return(0);
7 B+ C+ t8 m: X      }
; s$ }2 H( d2 w% l    p=ar*br[n-1]; q=-ai*bi[n-1];
' A5 z9 q5 O4 X- R! a7 q    s=(ar-ai)*(br[n-1]+bi[n-1]);
    br[n-1]=(p-q)/d; bi[n-1]=(s-p-q)/d;
8 Y& Q# b# O% j    for (i=n-2;i&gt;=0;i--)
    for (j=i+1;j&lt;=n-1;j++)
      { u=i*n+j;
! `0 Z4 h2 \8 W        p=ar*br[j]; q=ai*bi[j];
        s=(ar+ai)*(br[j]+bi[j]);
        br=br-p+q;
        bi=bi-s+p+q;
      }
3 v: g" f# I3 @+ q2 b4 X3 h    js[n-1]=n-1;
' j+ W6 L) W/ O7 _- p    for (k=n-1;k&gt;=0;k--)
3 T4 K7 S6 t$ }" A      if (js[k]!=k)
" ~6 F5 H$ }' A0 B7 ^6 z        { p=br[k]; br[k]=br[js[k]]; br[js[k]]=p;
+ e, c* E; V0 d5 X          p=bi[k]; bi[k]=bi[js[k]]; bi[js[k]]=p;
        }
    free(js);
    return(1);
# R" t- D: }. g9 C8 t: P5 a- v  }</P><>平方根法</P><>#include "math.h"
  #include "stdio.h"
& {# ?$ W4 v' A  x, N' p' ^  int achol(a,n,m,d)
/ n8 V6 V6 X$ j0 u  int n,m;
+ m+ Z8 K8 f5 D; y  double a[],d[];
* T8 M2 F, e7 ~  { int i,j,k,u,v;
- m" v0 V1 M' M3 ~5 }' b    if ((a[0]+1.0==1.0)||(a[0]&lt;0.0))
      { printf("fail\n"); return(-2);}
0 Q0 Y9 N: x# U* b' S5 k! z    a[0]=sqrt(a[0]);
    for (j=1; j&lt;=n-1; j++) a[j]=a[j]/a[0];
( a' E4 e0 S3 u0 a4 i, f1 w    for (i=1; i&lt;=n-1; i++)
+ Y, Z/ a/ o, R& M$ a      { u=i*n+i;
        for (j=1; j&lt;=i; j++)
          { v=(j-1)*n+i;
            a=a-a[v]*a[v];
: E3 y+ v, G" c% _4 @0 R# {          }
        if ((a+1.0==1.0)||(a&lt;0.0))
; @8 e8 d# J4 |% ~. h3 J          { printf("fail\n"); return(-2);}
6 ]  W2 [/ j8 t- ?' r- P- D! H        a=sqrt(a);
        if (i!=(n-1))
          { for (j=i+1; j&lt;=n-1; j++)
' Y0 p5 y& ^1 l              { v=i*n+j;
                for (k=1; k&lt;=i; k++)
. x$ z7 W( I' g' y6 i: M& p                  a[v]=a[v]-a[(k-1)*n+i]*a[(k-1)*n+j];
" K9 ?/ [) l# |- n; [                a[v]=a[v]/a;
              }
4 n, f9 j! G- n4 r7 G5 F/ s2 s          }
0 |& o% D& E. e2 [      }
    for (j=0; j&lt;=m-1; j++)
      { d[j]=d[j]/a[0];
/ h0 E0 U8 {1 d+ W; B; f/ E( w5 v        for (i=1; i&lt;=n-1; i++)
( [, E" U0 x% `          { u=i*n+i; v=i*m+j;
, h2 ]1 B, M3 L! l            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;
5 ]' P' |# r4 l/ ^/ e          }
7 w5 Z- a5 O# g      }
7 w- L* F4 G( g    for (j=0; j&lt;=m-1; j++)
  m: |' h; C: I      { u=(n-1)*m+j;
        d=d/a[n*n-1];
        for (k=n-1; k&gt;=1; k--)
          { u=(k-1)*m+j;
            for (i=k; i&lt;=n-1; i++)
- P" a% f5 h/ t) Z7 ?              { v=(k-1)*n+i;
                d=d-a[v]*d[i*m+j];
              }
            v=(k-1)*n+k-1;
2 ?/ F! T  n) h2 B            d=d/a[v];
          }
( `- u* Z7 I% m      }
  _* k; c; D+ w3 j    return(2);
4 y" g7 R) S' u( P9 ]  }</P><>牛顿向前，向后插值法</P><>double eelgr(x0,h,n,y,t)
% N3 v1 G9 v& i, f9 y3 T  int n;
  double x0,h,t,y[];
  { int i,j,k,m;
    double z,s,xi,xj;
    float p,q;
. Q$ h4 r, m% C2 l$ _0 x    z=0.0;
4 M& A0 r+ y9 N# u    if (n&lt;1) return(z);
    if (n==1) { z=y[0]; return(z);}
    if (n==2)
1 A0 {+ P/ g, s! c  M      { 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;
7 D: ]& v7 h7 y6 M5 ]        if (p&gt;q) i=i+1;
7 o* O/ c6 B2 F( \. L* V      }
: B7 J  t# S1 E# J    else i=0;
" I* V- d9 A+ V% c# w( `2 c4 Q" a1 G    k=i-4;
    if (k&lt;0) k=0;
    m=i+3;
    if (m&gt;n-1) m=n-1;
! Y. ], H. A2 U4 G) T    for (i=k;i&lt;=m;i++)
6 a+ V& I- z  T) ?- F1 n6 @      { s=1.0; xi=x0+i*h;
/ C2 C3 I& b. l* f  S        for (j=k; j&lt;=m; j++)
$ F7 M5 K  t3 ]* J0 L9 |2 ?  z          if (j!=i)
& k" F! N' o" Y            { xj=x0+j*h;
              s=s*(t-xj)/(xi-xj);
  c. b5 _# ?! e/ J" s/ D            }
        z=z+s*y;
, Q; n) j+ h' e0 P9 a# x      }
    return(z);
  }
向前，向后是一样的思想！！</P><>加权最小二乘法</P><>#include "math.h"
) V7 i: N( J6 ^2 P. r  void hpir1(x,y,n,a,m,dt)
" b! l1 I/ J# G' d  int n,m;
  double x[],y[],a[],dt[];
8 F0 b' b: [! g  { int i,j,k;
/ l0 P4 L: O- h6 M# g' ^    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;
    if (m&gt;20) m=20;
    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;
6 T; R  \8 b  i* I% X( q    for (i=0; i&lt;=n-1; i++)
      { p=p+(x-z); c=c+y;}
    c=c/d1; p=p/d1;
9 ^; U9 [9 N: ~2 \$ V4 ~    a[0]=c*b[0];
    if (m&gt;1)
      { t[1]=1.0; t[0]=-p;
/ f1 d5 V: `* D6 Q6 N        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;
  {8 |& Y. g: ]            c=c+y*q;
# L" C9 {: _  {4 A% u0 H            g=g+(x-z)*q*q;
) g. h6 d5 W9 f: o' P$ V2 D) N          }
: E0 f4 p. f, \+ T1 c9 x        c=c/d2; p=g/d2; q=d2/d1;
: D7 {- l% E: I9 Z5 K        d1=d2;
! Q( `9 I8 U- u+ Z& Y        a[1]=c*t[1]; a[0]=c*t[0]+a[0];
+ Y; X: a" a3 K      }
    for (j=2; j&lt;=m-1; j++)
      { s[j]=t[j-1];
        s[j-1]=-p*t[j-1]+t[j-2];
        if (j&gt;=3)
, M& ^$ X$ P9 [" ~' ?5 o9 F) ?. H9 K          for (k=j-2; k&gt;=1; k--)
1 R# V2 ]; i" f( b            s[k]=-p*t[k]+t[k-1]-q*b[k];
        s[0]=-p*t[0]-q*b[0];
, ^9 I8 M2 k/ t( Z0 g        d2=0.0; c=0.0; g=0.0;
        for (i=0; i&lt;=n-1; i++)
% U5 f' w! i$ Y8 ~& M5 \          { q=s[j];
) I3 A3 S) h( \3 B7 [3 j. W1 m            for (k=j-1; k&gt;=0; k--)
              q=q*(x-z)+s[k];
& Z; K4 _4 m' X8 v; {! R            d2=d2+q*q; c=c+y*q;
6 K. p$ |0 ]! q+ ^: z            g=g+(x-z)*q*q;
! `+ D  k% }1 @& T- O3 ^          }
, C+ B$ K* c) ^  J        c=c/d2; p=g/d2; q=d2/d1;
  P9 t5 [0 f; Z7 g9 W        d1=d2;
        a[j]=c*s[j]; t[j]=s[j];
6 M7 B( F+ D. V/ I6 |' R  \        for (k=j-1; k&gt;=0; k--)
          { a[k]=c*s[k]+a[k];
            b[k]=t[k]; t[k]=s[k];
          }
4 }# V5 B* Q" [      }
+ e8 c. T) c; e2 W    dt[0]=0.0; dt[1]=0.0; dt[2]=0.0;
    for (i=0; i&lt;=n-1; i++)
5 |) P3 @; [( ?' y% q3 {6 F% k      { q=a[m-1];
, N, Y. u( r( J' f' v0 B7 P, V' \        for (k=m-2; k&gt;=0; k--)
          q=a[k]+q*(x-z);
        p=q-y;
" ?. g# y7 ]8 R0 \        if (fabs(p)&gt;dt[2]) dt[2]=fabs(p);
! Q5 _9 S9 n, M' c1 z        dt[0]=dt[0]+p*p;
        dt[1]=dt[1]+fabs(p);
" x; @( t$ B( y- }      }
9 `2 @" D% `; b# J9 Z, k- y0 H    return;
# d# m5 t( b8 S' i6 ?7 l2 z  }</P><>龙贝格积分法</P><>#include "math.h"
' q0 b3 P* c! a" ~; ?/ v  double fromb(a,b,eps)
  double a,b,eps;
: [( f1 F: M) V7 [  { extern double frombf();
    int m,n,i,k;
4 q& k+ Q' L4 _' E, c    double y[10],h,ep,p,x,s,q;
: ?+ P( ~( w# z! x0 J+ U    h=b-a;
    y[0]=h*(frombf(a)+frombf(b))/2.0;
$ M; T7 Q& Y' i+ `8 r' w    m=1; n=1; ep=eps+1.0;
    while ((ep&gt;=eps)&amp;&amp;(m&lt;=9))
0 y& E$ a0 [# c3 H5 n" |% a1 ~      { p=0.0;
0 I9 \3 q: o4 P- i        for (i=0;i&lt;=n-1;i++)
; X. ~: H6 T% U/ y+ y          { x=a+(i+0.5)*h;
/ S5 @4 S$ u8 a0 @5 X            p=p+frombf(x);
- L0 z( }8 r3 t- \* Q, ?. p  i          }
) z, S, n0 Z9 A; N/ {        p=(y[0]+h*p)/2.0;
        s=1.0;
3 @% n0 m5 r$ b4 `) F- t        for (k=1;k&lt;=m;k++)
          { s=4.0*s;
            q=(s*p-y[k-1])/(s-1.0);
            y[k-1]=p; p=q;
& ^) m5 _# Y8 h$ U% J% s% {* d          }
        ep=fabs(q-y[m-1]);
        m=m+1; y[m-1]=q; n=n+n; h=h/2.0;
  J- H$ A/ }6 E2 \9 e      }
7 p4 [. a' o6 K% `    return(q);
  }</P><>呵呵 希望对你有用！！</P>

西西

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

陈了

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