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

西西

<>本人紧急求助，希望各位大虾们伸出你们友情之手，帮我一下。</P>
" K, J3 K! @8 A<>1。牛顿迭代法</P>
<>1。全主消元法</P>
<>3。改进平方根法</P>
<>4。牛顿向前，向后插值法</P>
<>5。加权最小二乘法</P>
<>6。龙贝格积分法</P>
" {# L$ q! _$ c6 M<>拜托了。。。。</P>

loveyuan

<>1。牛顿迭代法</P><>
  #include "stdio.h"
  #include "math.h"
- _0 M' j0 B: z$ N' L+ n  int dnewt(x,eps,js)
+ f0 T+ B' |$ A3 [7 V  a+ ^# L6 _  int js;
  double *x,eps;
  { extern void dnewtf();
    int k,l;
9 e+ D$ x  C! a; {3 c+ Z    double y[2],d,p,x0,x1;
( G3 e: a2 b# O$ {    l=js; k=1; x0=*x;
+ j4 k! e0 H4 M# L8 k1 b( Y    dnewtf(x0,y);
    d=eps+1.0;
    while ((d&gt;=eps)&amp;&amp;(l!=0))
# P% G( ?9 E% p& Z1 m' e- @+ P5 Y3 b6 N+ g      { if (fabs(y[1])+1.0==1.0)
          { printf("err\n"); return(-1);}
        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;
      }
" e7 q, w0 Q* ?5 u) ?    *x=x1;
8 U. q+ K, _/ ?# H    k=js-l;
    return(k);
  }</P><>全主消元法</P><>#include "stdlib.h"
5 ^; v9 o9 ^+ |! B  #include "stdio.h"
. M& j& z  H  ~  int acgas(ar,ai,n,br,bi)
8 N4 S: a5 o$ Y) `1 L1 E  int n;
  double ar[],ai[],br[],bi[];
  { int *js,l,k,i,j,is,u,v;
0 z' q4 v$ \6 j+ {% F; d    double p,q,s,d;
) X! W# p1 g: d) r" `    js=malloc(n*sizeof(int));
( S5 {5 ?1 E5 t  B    for (k=0;k&lt;=n-2;k++)
, F. n2 x8 X6 d; A5 I- T/ N0 ]3 @      { d=0.0;
# B# B: p$ z$ x8 ]. o1 [3 j0 I" E        for (i=k;i&lt;=n-1;i++)
/ }1 n- C9 i0 x3 Y/ f: r        for (j=k;j&lt;=n-1;j++)
          { u=i*n+j;
4 y/ b# D" u+ E, H/ e            p=ar*ar+ai*ai;
            if (p&gt;d) {d=p;js[k]=j;is=i;}
          }
5 {* i1 p$ R2 `! w  n        if (d+1.0==1.0)
          { free(js); printf("err**fail\n");
/ u+ V7 e" i. N# I            return(0);
          }
4 m) T8 Q/ }* i& O7 {9 j/ l        if (is!=k)
          { 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;
, I+ K4 }% b0 b/ A; b$ O              }
% Q& J1 e5 e+ k2 M8 {& W            p=br[k]; br[k]=br[is]; br[is]=p;
2 j) `% R' s+ H! _  j8 }) D            p=bi[k]; bi[k]=bi[is]; bi[is]=p;
          }
        if (js[k]!=k)
          for (i=0;i&lt;=n-1;i++)
            { u=i*n+k; v=i*n+js[k];
              p=ar; ar=ar[v]; ar[v]=p;
4 A- B' l  N4 a) _' o. a' V, ~/ ~7 B              p=ai; ai=ai[v]; ai[v]=p;
            }
        v=k*n+k;
! [) w3 I- D# N* `: k/ {  P, @        for (j=k+1;j&lt;=n-1;j++)
          { u=k*n+j;
$ W, n* \  ^; X% o9 L7 v            p=ar*ar[v]; q=-ai*ai[v];
% X5 r% ]1 z% L            s=(ar[v]-ai[v])*(ar+ai);
- {1 `( n8 R. r& d            ar=(p-q)/d; ai=(s-p-q)/d;
* m7 X7 _9 {- u( e3 A          }
        p=br[k]*ar[v]; q=-bi[k]*ai[v];
        s=(ar[v]-ai[v])*(br[k]+bi[k]);
7 J2 @. M2 t* r; T        br[k]=(p-q)/d; bi[k]=(s-p-q)/d;
* m1 P: u8 D" V) H$ d        for (i=k+1;i&lt;=n-1;i++)
7 K, X# K7 P2 Y1 }          { u=i*n+k;
7 D9 c' v7 y! y; n7 I6 e9 Y            for (j=k+1;j&lt;=n-1;j++)
              { v=k*n+j; l=i*n+j;
                p=ar*ar[v]; q=ai*ai[v];
                s=(ar+ai)*(ar[v]+ai[v]);
6 J% D- Q! D$ D* d) e                ar[l]=ar[l]-p+q;
                ai[l]=ai[l]-s+p+q;
              }
            p=ar*br[k]; q=ai*bi[k];
- x' {" I  G% o. \3 n            s=(ar+ai)*(br[k]+bi[k]);
, A; v$ M9 i1 o% d- Z7 s            br=br-p+q; bi=bi-s+p+q;
( ~) `2 r. g1 ^& O* Z4 S+ r          }
0 z7 I3 L! c! W6 e1 O$ U      }
4 F& }; g9 s: j+ O* j# i- ~; l+ p( F/ J    u=(n-1)*n+n-1;
    d=ar*ar+ai*ai;
    if (d+1.0==1.0)
      { free(js); printf("err**fail\n");
8 A4 X* W; a5 I# h- g# o        return(0);
      }
' q. F" X! n# \    p=ar*br[n-1]; q=-ai*bi[n-1];
% t+ _4 F, t. |( x2 f. z; C    s=(ar-ai)*(br[n-1]+bi[n-1]);
    br[n-1]=(p-q)/d; bi[n-1]=(s-p-q)/d;
    for (i=n-2;i&gt;=0;i--)
& ~. G! x# q; V, F1 z9 E: @+ v: L6 Y5 z    for (j=i+1;j&lt;=n-1;j++)
      { u=i*n+j;
( T  G$ O& K: ~$ V        p=ar*br[j]; q=ai*bi[j];
8 }+ T8 B2 W5 a5 c6 s        s=(ar+ai)*(br[j]+bi[j]);
3 J& Z( _# e# z        br=br-p+q;
1 u9 q4 ]+ V6 o% `& A0 K        bi=bi-s+p+q;
      }
3 Y( }& R* o! a9 R& |3 |    js[n-1]=n-1;
    for (k=n-1;k&gt;=0;k--)
3 R1 D- f" z3 e/ x7 Q6 Z0 r      if (js[k]!=k)
' `, m% B" d; A  E" a        { p=br[k]; br[k]=br[js[k]]; br[js[k]]=p;
          p=bi[k]; bi[k]=bi[js[k]]; bi[js[k]]=p;
  ?5 ^, y. N% P/ ^# N        }
    free(js);
( t: W) R0 W  L$ O3 j2 u    return(1);
& _/ T- h3 C8 D- o% L9 V/ g% H( F. x  }</P><>平方根法</P><>#include "math.h"
  #include "stdio.h"
" B5 A" x; j# p8 }8 X+ Q* H7 C  int achol(a,n,m,d)
) t9 e6 \# s; l) \: w* c8 L  int n,m;
  double a[],d[];
  { 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]);
( D6 E! m. L$ A4 Q    for (j=1; j&lt;=n-1; j++) a[j]=a[j]/a[0];
9 H. L1 s! \3 ?) }; N    for (i=1; i&lt;=n-1; i++)
" @- J+ S3 b# u      { u=i*n+i;
        for (j=1; j&lt;=i; j++)
4 K# r- |. X4 X5 |          { v=(j-1)*n+i;
9 N" `3 C  m8 E, \; d+ {            a=a-a[v]*a[v];
          }
" M6 ~1 x) {' C. f  M# w0 B        if ((a+1.0==1.0)||(a&lt;0.0))
! D! y7 O! t# @          { printf("fail\n"); return(-2);}
7 N' o* X: ^: x! t        a=sqrt(a);
        if (i!=(n-1))
          { for (j=i+1; j&lt;=n-1; j++)
) h. I" A2 X" v6 C3 _3 O# f* {. t              { v=i*n+j;
/ j2 n: V+ e# Z" L& Z1 [4 A( l                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;
  l* ]7 s# p5 @              }
          }
* H. ?2 [: Q1 E6 H      }
+ I% R, h+ ]8 v( W. k* Z    for (j=0; j&lt;=m-1; j++)
; l! K) d7 @; c" S      { d[j]=d[j]/a[0];
/ V! _1 k1 V9 Y" z7 n        for (i=1; i&lt;=n-1; i++)
          { u=i*n+i; v=i*m+j;
7 b( L' t) _1 l! q$ R9 }            for (k=1; k&lt;=i; k++)
0 f. K  Q1 y' P) T" u# x' T              d[v]=d[v]-a[(k-1)*n+i]*d[(k-1)*m+j];
            d[v]=d[v]/a;
          }
      }
    for (j=0; j&lt;=m-1; j++)
      { u=(n-1)*m+j;
. M1 A" A1 J( V$ |% n        d=d/a[n*n-1];
        for (k=n-1; k&gt;=1; k--)
  K6 l6 \1 }2 u          { u=(k-1)*m+j;
6 M8 u9 j5 T6 N! \$ q1 j            for (i=k; i&lt;=n-1; i++)
6 e7 @7 _& j' U6 j              { v=(k-1)*n+i;
) N/ `1 Q" ^9 e0 @6 J                d=d-a[v]*d[i*m+j];
              }
/ \4 B+ X. `, g( \8 e$ G            v=(k-1)*n+k-1;
            d=d/a[v];
          }
3 ^$ h! u' W/ J6 f2 ~      }
0 Q# u' @; B4 j2 b. N6 Q    return(2);
  }</P><>牛顿向前，向后插值法</P><>double eelgr(x0,h,n,y,t)
  int n;
  double x0,h,t,y[];
  { int i,j,k,m;
    double z,s,xi,xj;
    float p,q;
    z=0.0;
) v" ^! O9 `# r% `# H# R' e    if (n&lt;1) return(z);
0 T7 Q. y% G  x    if (n==1) { z=y[0]; return(z);}
    if (n==2)
) V$ L6 C- U& J      { z=(y[1]*(t-x0)-y[0]*(t-x0-h))/h;
3 f1 d. n7 ^1 B/ n        return(z);
      }
+ N  \+ m5 [& P. j3 n& q& @    if (t&gt;x0)
- @: S4 }2 w- e; i      { p=(t-x0)/h; i=(int)p; q=(float)i;
        if (p&gt;q) i=i+1;
1 x( g' Y5 s$ O' g      }
) g4 @  Z) h+ H9 t9 Y6 {    else i=0;
    k=i-4;
, w& u6 L2 T2 v% v  P8 u4 ~    if (k&lt;0) k=0;
    m=i+3;
    if (m&gt;n-1) m=n-1;
- F, F* u: \' v( V. F  e    for (i=k;i&lt;=m;i++)
7 K6 }) {9 x* R# N- b      { s=1.0; xi=x0+i*h;
# f2 K# ~) L8 Y( X) R        for (j=k; j&lt;=m; j++)
$ Y: e% @+ E! W+ [          if (j!=i)
4 @  y( Q+ t1 [            { xj=x0+j*h;
              s=s*(t-xj)/(xi-xj);
            }
8 s- J) ~: H4 v        z=z+s*y;
      }
    return(z);
  }
) [! ], s' M# h2 S向前，向后是一样的思想！！</P><>加权最小二乘法</P><>#include "math.h"
1 i0 c0 b3 Q! k: p& N  i' l# K  void hpir1(x,y,n,a,m,dt)
  int n,m;
  double x[],y[],a[],dt[];
5 s( r7 {1 E% R- `. E  { int i,j,k;
/ v) ]6 d) [1 u0 E8 p6 I- 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;
    if (m&gt;n) m=n;
    if (m&gt;20) m=20;
    z=0.0;
  t# a; R4 i2 E- R- m" q    for (i=0; i&lt;=n-1; i++) z=z+x/(1.0*n);
% i1 k' x9 b8 l# n, ?3 a% j    b[0]=1.0; d1=1.0*n; p=0.0; c=0.0;
    for (i=0; i&lt;=n-1; i++)
      { p=p+(x-z); c=c+y;}
) V( R) S- K) V* X    c=c/d1; p=p/d1;
    a[0]=c*b[0];
9 o. r( E+ @8 C0 x9 P7 o. g+ ?1 K    if (m&gt;1)
      { t[1]=1.0; t[0]=-p;
        d2=0.0; c=0.0; g=0.0;
) T. f0 V' k! e; y        for (i=0; i&lt;=n-1; i++)
) O$ a! z! t- ^) i4 Y, Q9 X1 ^" R7 D          { q=x-z-p; d2=d2+q*q;
0 |6 S5 e& T# ?. q5 @2 S2 Y# l            c=c+y*q;
            g=g+(x-z)*q*q;
* w4 A" W: ^# Z          }
* b! q& H, c* B% j8 {        c=c/d2; p=g/d2; q=d2/d1;
        d1=d2;
0 H% r2 C% m! Q/ L; y, {) b" I        a[1]=c*t[1]; a[0]=c*t[0]+a[0];
      }
# x5 R& P, d2 k% r    for (j=2; j&lt;=m-1; j++)
9 _& ~& L7 g5 j- H) B' W" b& {      { s[j]=t[j-1];
# u2 k$ t( @3 `5 I, [        s[j-1]=-p*t[j-1]+t[j-2];
" z5 K2 S) o: e" _' S        if (j&gt;=3)
8 v- J, P9 x  v7 A+ z          for (k=j-2; k&gt;=1; k--)
/ W" s+ g) L) J) `. ~- ~9 x            s[k]=-p*t[k]+t[k-1]-q*b[k];
        s[0]=-p*t[0]-q*b[0];
. k  E% c9 Q3 D) j. S        d2=0.0; c=0.0; g=0.0;
: s6 w% V0 E; P' P& h/ T; y7 @        for (i=0; i&lt;=n-1; i++)
( m1 G. g7 q; f3 H  m          { q=s[j];
            for (k=j-1; k&gt;=0; k--)
* T0 M0 o& H& S; E: ~; |5 A$ I( c              q=q*(x-z)+s[k];
2 P9 M, P- e0 }0 r, `6 v% ~3 U( C            d2=d2+q*q; c=c+y*q;
: r4 ]; R- T( U1 n- k. P# E            g=g+(x-z)*q*q;
2 Q) P4 h% \( K" k          }
% d# Z' t# K# B! {        c=c/d2; p=g/d2; q=d2/d1;
: _8 w& i9 }- p        d1=d2;
2 U* |5 x! D9 y% u6 z! k$ L1 |        a[j]=c*s[j]; t[j]=s[j];
8 k, @* r* l7 I# g# U# a% {        for (k=j-1; k&gt;=0; k--)
. Q/ m) ?' v1 n) S. p2 B          { a[k]=c*s[k]+a[k];
            b[k]=t[k]; t[k]=s[k];
0 E  Y4 x( z+ i7 A          }
% q& G- L# K/ M' C3 V      }
    dt[0]=0.0; dt[1]=0.0; dt[2]=0.0;
* l. Z* J# E) m: b/ q: Y& c    for (i=0; i&lt;=n-1; i++)
      { q=a[m-1];
        for (k=m-2; k&gt;=0; k--)
          q=a[k]+q*(x-z);
3 Z' H5 x7 h  h1 H4 |* y! M        p=q-y;
5 A7 ?) _/ V# h  ]        if (fabs(p)&gt;dt[2]) dt[2]=fabs(p);
        dt[0]=dt[0]+p*p;
8 @+ u! P% S. M; w* [        dt[1]=dt[1]+fabs(p);
3 y6 ]" g7 Q, v+ D- u" K2 _% m  p      }
    return;
  }</P><>龙贝格积分法</P><>#include "math.h"
, g9 ?: p" ?8 {( c0 |  double fromb(a,b,eps)
  double a,b,eps;
  { extern double frombf();
    int m,n,i,k;
    double y[10],h,ep,p,x,s,q;
$ f7 f$ _% h, n; ?( a3 k    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))
+ j: o2 d+ y* C& A* @      { p=0.0;
        for (i=0;i&lt;=n-1;i++)
          { x=a+(i+0.5)*h;
            p=p+frombf(x);
$ N/ }2 S7 l; R' w# n          }
        p=(y[0]+h*p)/2.0;
        s=1.0;
+ C) s# d  p7 f# w' t* J% Z0 s        for (k=1;k&lt;=m;k++)
          { s=4.0*s;
+ }) o8 b- B2 Y' V: V" d            q=(s*p-y[k-1])/(s-1.0);
            y[k-1]=p; p=q;
          }
        ep=fabs(q-y[m-1]);
: x1 i% b: R- ]) [        m=m+1; y[m-1]=q; n=n+n; h=h/2.0;
      }
    return(q);
  }</P><>呵呵 希望对你有用！！</P>

西西

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

陈了

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