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

西西

<>本人紧急求助，希望各位大虾们伸出你们友情之手，帮我一下。</P>
<>1。牛顿迭代法</P>
4 E! d4 j9 L9 w- W( ^) z% r<>1。全主消元法</P>
<>3。改进平方根法</P>
<>4。牛顿向前，向后插值法</P>
<>5。加权最小二乘法</P>
<>6。龙贝格积分法</P>
<>拜托了。。。。</P>

loveyuan

<>1。牛顿迭代法</P><>
  #include "stdio.h"
6 L  i$ K* ]4 Y  #include "math.h"
# V3 \( I$ W1 F( J- A  int dnewt(x,eps,js)
! X1 @! u- n  u! a3 u! T; ~  int js;
  double *x,eps;
  { extern void dnewtf();
    int k,l;
1 ?& i9 y6 J4 k, M6 L: j    double y[2],d,p,x0,x1;
    l=js; k=1; x0=*x;
    dnewtf(x0,y);
    d=eps+1.0;
    while ((d&gt;=eps)&amp;&amp;(l!=0))
. s8 |+ b, t- r  U      { if (fabs(y[1])+1.0==1.0)
" u- k+ G5 C! f( J7 y- |3 p9 u; }          { printf("err\n"); return(-1);}
2 ]* R; t; z7 Y7 [, L        x1=x0-y[0]/y[1];
1 j) C+ l4 [& l! M% X        dnewtf(x1,y);
        d=fabs(x1-x0); p=fabs(y[0]);
        if (p&gt;d) d=p;
0 [3 S! ^; k  d        x0=x1; l=l-1;
0 f2 Y) _# r3 c) @% t      }
    *x=x1;
    k=js-l;
    return(k);
: G& Q8 L: {& P1 B( f6 u0 k  }</P><>全主消元法</P><>#include "stdlib.h"
  #include "stdio.h"
  int acgas(ar,ai,n,br,bi)
2 N$ b2 o: {9 A' X! j+ F  int n;
  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;
* m# }9 U5 h8 I& q2 K' W/ Z( K        for (i=k;i&lt;=n-1;i++)
2 ^7 Z$ f* u9 X" _& V        for (j=k;j&lt;=n-1;j++)
          { u=i*n+j;
            p=ar*ar+ai*ai;
4 }. k0 \0 b; I. Y$ p4 ^% {            if (p&gt;d) {d=p;js[k]=j;is=i;}
, V& S& s3 [! Q# }: T/ V# H          }
        if (d+1.0==1.0)
# p8 M; r2 r. }0 k& f7 V7 S" C          { free(js); printf("err**fail\n");
0 p6 d4 o' c" g            return(0);
          }
2 @& k! N6 a3 x4 E& Z  k; y        if (is!=k)
  g+ y3 R% F: o2 J. o4 W* ]/ R          { for (j=k;j&lt;=n-1;j++)
7 b8 |* o  f3 |! a) p* }9 y              { u=k*n+j; v=is*n+j;
                p=ar; ar=ar[v]; ar[v]=p;
                p=ai; ai=ai[v]; ai[v]=p;
              }
            p=br[k]; br[k]=br[is]; br[is]=p;
            p=bi[k]; bi[k]=bi[is]; bi[is]=p;
1 [5 s0 Q- b4 V4 P' d  y) H          }
0 \& [5 E( P5 C1 H' f4 G, _  F) `' N        if (js[k]!=k)
          for (i=0;i&lt;=n-1;i++)
  X. z  V) d1 s3 a9 [3 j            { u=i*n+k; v=i*n+js[k];
: O: P! r- U3 n; S" a" W" Z              p=ar; ar=ar[v]; ar[v]=p;
$ A+ }# ^& x' G7 M$ _              p=ai; ai=ai[v]; ai[v]=p;
; k: \/ q' B: S( {+ Z            }
        v=k*n+k;
        for (j=k+1;j&lt;=n-1;j++)
6 i" I9 Y  l$ ^2 g          { u=k*n+j;
* x+ K) E; \- B3 c  x+ [% W/ Q            p=ar*ar[v]; q=-ai*ai[v];
            s=(ar[v]-ai[v])*(ar+ai);
            ar=(p-q)/d; ai=(s-p-q)/d;
          }
' S, ^4 ~4 H, m        p=br[k]*ar[v]; q=-bi[k]*ai[v];
        s=(ar[v]-ai[v])*(br[k]+bi[k]);
, [% F! L6 K, Y6 n        br[k]=(p-q)/d; bi[k]=(s-p-q)/d;
1 m! @/ E; a$ R2 X        for (i=k+1;i&lt;=n-1;i++)
          { u=i*n+k;
! g& U. v" C& ?* R3 L0 e- R            for (j=k+1;j&lt;=n-1;j++)
) E3 v- l  p! d4 _: V              { v=k*n+j; l=i*n+j;
# h9 |1 _  d; ]6 ]4 ]. n                p=ar*ar[v]; q=ai*ai[v];
                s=(ar+ai)*(ar[v]+ai[v]);
4 f; f* B. n( b/ p+ Z/ B: p                ar[l]=ar[l]-p+q;
, l5 G" w8 Q+ t) i9 }                ai[l]=ai[l]-s+p+q;
              }
+ R' |2 I2 V/ n# h, z. |            p=ar*br[k]; q=ai*bi[k];
            s=(ar+ai)*(br[k]+bi[k]);
1 A" ^% B" x. G+ E% I            br=br-p+q; bi=bi-s+p+q;
          }
      }
/ `6 a2 O) q6 {" F7 @# P. d    u=(n-1)*n+n-1;
( M1 _& M# h# i9 ~- \6 _    d=ar*ar+ai*ai;
# m- S2 f0 S4 D0 c    if (d+1.0==1.0)
5 G$ b: u, U: n5 p! z" H( ]  o      { free(js); printf("err**fail\n");
        return(0);
      }
2 X( s5 J2 r, C$ \    p=ar*br[n-1]; q=-ai*bi[n-1];
+ b, u% f1 e7 Y    s=(ar-ai)*(br[n-1]+bi[n-1]);
3 j" [  x* a! h$ a9 I) P7 E  Z7 {    br[n-1]=(p-q)/d; bi[n-1]=(s-p-q)/d;
    for (i=n-2;i&gt;=0;i--)
# T. d: x( X5 x; V    for (j=i+1;j&lt;=n-1;j++)
      { u=i*n+j;
        p=ar*br[j]; q=ai*bi[j];
' X! F: t- n& T6 L        s=(ar+ai)*(br[j]+bi[j]);
        br=br-p+q;
0 v" I' l/ I" R+ k0 C        bi=bi-s+p+q;
4 D: d" N8 v5 g) K      }
    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;
          p=bi[k]; bi[k]=bi[js[k]]; bi[js[k]]=p;
6 X0 T6 _2 ~, ?7 ?; G5 q7 [        }
    free(js);
    return(1);
  }</P><>平方根法</P><>#include "math.h"
. r/ }0 I+ @6 ?  #include "stdio.h"
  int achol(a,n,m,d)
  int n,m;
  double a[],d[];
, n1 r6 N) i, q1 {, @5 O  { int i,j,k,u,v;
. `* S1 e/ K7 z: I    if ((a[0]+1.0==1.0)||(a[0]&lt;0.0))
* i6 k- {6 q* u  A      { printf("fail\n"); return(-2);}
+ x3 M1 D0 T; U+ h( Q! Z    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++)
% R1 E$ }2 ?1 N# S      { u=i*n+i;
+ q0 F4 E/ c( O* j2 R$ M: C        for (j=1; j&lt;=i; j++)
          { v=(j-1)*n+i;
) P( J6 R7 `: S7 D- Q! w            a=a-a[v]*a[v];
          }
        if ((a+1.0==1.0)||(a&lt;0.0))
" D9 C2 {2 f- a/ J: O          { printf("fail\n"); return(-2);}
6 ^* i; ^6 u( M+ |6 [- ?3 W, k        a=sqrt(a);
" u7 y/ V/ B' l4 |; o        if (i!=(n-1))
% v9 H1 g: X$ B  m7 _          { for (j=i+1; j&lt;=n-1; j++)
/ o  v( R5 V0 T0 f              { v=i*n+j;
                for (k=1; k&lt;=i; k++)
5 V/ Q$ y; b! H4 d$ Q1 `! ~                  a[v]=a[v]-a[(k-1)*n+i]*a[(k-1)*n+j];
                a[v]=a[v]/a;
              }
          }
      }
) G; a' Q, S# S6 A( o    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;
7 e9 Y2 b0 T  \            for (k=1; k&lt;=i; k++)
4 E5 S8 B& l: `8 u) W              d[v]=d[v]-a[(k-1)*n+i]*d[(k-1)*m+j];
, H& ?, j( z4 s( m. g* H! T# o            d[v]=d[v]/a;
          }
      }
    for (j=0; j&lt;=m-1; j++)
  L/ t( e+ d! E2 ]8 I8 g      { u=(n-1)*m+j;
% U4 y* o! ?/ U) `, P4 ^/ S% V* W: {        d=d/a[n*n-1];
        for (k=n-1; k&gt;=1; k--)
          { u=(k-1)*m+j;
: W1 ~4 n) y4 ]# [* j! w8 R            for (i=k; i&lt;=n-1; i++)
& U0 i' j' a3 G              { v=(k-1)*n+i;
# r! q" o! `8 P                d=d-a[v]*d[i*m+j];
              }
9 ]  c' S; U2 T- S& ?* n% ]            v=(k-1)*n+k-1;
* T0 }( h" q1 E8 m( T- A( n# W            d=d/a[v];
          }
6 V3 W, a- q( m' J' V- ]      }
* R% k# y6 e' ?( ?% k2 Q    return(2);
  }</P><>牛顿向前，向后插值法</P><>double eelgr(x0,h,n,y,t)
9 b+ p( N. Q+ ], G) }4 V* e5 @, s  int n;
9 T. F; P$ U$ h8 z7 v( }6 M" ^  double x0,h,t,y[];
  { int i,j,k,m;
    double z,s,xi,xj;
    float p,q;
    z=0.0;
    if (n&lt;1) return(z);
    if (n==1) { z=y[0]; return(z);}
    if (n==2)
      { z=(y[1]*(t-x0)-y[0]*(t-x0-h))/h;
/ }3 F8 T9 h& s8 N        return(z);
- [0 V* Z8 m' H4 f8 Q1 t! s1 l. ~      }
    if (t&gt;x0)
6 l9 u0 u$ u" R- Q% u2 i      { p=(t-x0)/h; i=(int)p; q=(float)i;
" T' q9 y3 C* U  O2 M. R        if (p&gt;q) i=i+1;
      }
+ \- z1 C! h; T- c4 E, N    else i=0;
0 ~4 _& O4 A3 q( _! K+ q7 ]    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;
% @- ]9 Z4 }# |& ^        for (j=k; j&lt;=m; j++)
$ Y1 [6 n7 J3 S1 R( @$ M$ w          if (j!=i)
7 H  Q$ Q& Q, M" ~0 d- I3 ]            { xj=x0+j*h;
# e! m; w. e' n& i/ D              s=s*(t-xj)/(xi-xj);
            }
2 ^2 E/ r8 F0 J  R4 q        z=z+s*y;
      }
7 i: ^7 u4 ~2 M- }9 h% _- ^/ [' p    return(z);
' [% w* W9 j  k# i- y  }
( Q1 M/ `3 S# Y( J0 {向前，向后是一样的思想！！</P><>加权最小二乘法</P><>#include "math.h"
  void hpir1(x,y,n,a,m,dt)
  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;
* c' a6 X3 h8 I3 _    if (m&gt;n) m=n;
6 C  p0 r5 O! b/ D$ H6 W7 C    if (m&gt;20) m=20;
& D  {$ D- t$ M& y    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;
* Y8 j7 B0 I# e    for (i=0; i&lt;=n-1; i++)
% s7 I- S* l1 S5 k      { p=p+(x-z); c=c+y;}
    c=c/d1; p=p/d1;
3 J% y* H% Z. [! w0 d& C    a[0]=c*b[0];
    if (m&gt;1)
4 c( _8 @. r  W0 i$ `, W      { t[1]=1.0; t[0]=-p;
        d2=0.0; c=0.0; g=0.0;
        for (i=0; i&lt;=n-1; i++)
1 H- b2 p  j& s, u$ i          { q=x-z-p; d2=d2+q*q;
            c=c+y*q;
, V$ N5 n5 ~6 x( S% {! {0 B            g=g+(x-z)*q*q;
          }
        c=c/d2; p=g/d2; q=d2/d1;
        d1=d2;
: x  k, B9 v6 W! x7 F        a[1]=c*t[1]; a[0]=c*t[0]+a[0];
      }
. D' d' R: Z8 R8 P+ Y, h0 Q7 y9 d, [    for (j=2; j&lt;=m-1; j++)
- J) B& b4 T+ S, t. M      { s[j]=t[j-1];
        s[j-1]=-p*t[j-1]+t[j-2];
2 k# i8 S$ W0 W( Y        if (j&gt;=3)
          for (k=j-2; k&gt;=1; k--)
            s[k]=-p*t[k]+t[k-1]-q*b[k];
$ {" p* O) m" ~) ]) }+ x& y        s[0]=-p*t[0]-q*b[0];
, r: x" P5 _5 A        d2=0.0; c=0.0; g=0.0;
        for (i=0; i&lt;=n-1; i++)
          { q=s[j];
            for (k=j-1; k&gt;=0; k--)
              q=q*(x-z)+s[k];
$ R  t. Y. W, V2 ?+ l            d2=d2+q*q; c=c+y*q;
            g=g+(x-z)*q*q;
          }
5 J/ t5 h- n" s+ _        c=c/d2; p=g/d2; q=d2/d1;
        d1=d2;
/ ~5 {. m- r  {' \& S6 S/ p9 n        a[j]=c*s[j]; t[j]=s[j];
        for (k=j-1; k&gt;=0; k--)
          { a[k]=c*s[k]+a[k];
: \) q, t+ K9 e" U, O# S1 P( [            b[k]=t[k]; t[k]=s[k];
; q* O- {! i- _1 Y, a5 O- r          }
      }
    dt[0]=0.0; dt[1]=0.0; dt[2]=0.0;
    for (i=0; i&lt;=n-1; i++)
      { q=a[m-1];
- ^: U( F8 w* O5 ]3 D        for (k=m-2; k&gt;=0; k--)
          q=a[k]+q*(x-z);
& b# v# E, J5 w; R/ @0 n* c9 L, W8 `        p=q-y;
; x3 {# b5 F% n) _        if (fabs(p)&gt;dt[2]) dt[2]=fabs(p);
        dt[0]=dt[0]+p*p;
3 @: D* n0 \3 S8 a8 z, ~        dt[1]=dt[1]+fabs(p);
      }
# r/ n1 G: ]2 ]  u    return;
  }</P><>龙贝格积分法</P><>#include "math.h"
8 K; J/ f# _# k& `  double fromb(a,b,eps)
$ `+ Q' t4 X5 l6 V$ e- u  double a,b,eps;
$ i9 S* Z4 {( j$ e  n1 W  { extern double frombf();
    int m,n,i,k;
5 z, d: H# Z4 [" ~- a! D6 w% k    double y[10],h,ep,p,x,s,q;
8 h% L% e& V7 U    h=b-a;
* ]% M, O$ M$ F9 B    y[0]=h*(frombf(a)+frombf(b))/2.0;
    m=1; n=1; ep=eps+1.0;
' S4 ]1 J) g1 l& z& \" c% x    while ((ep&gt;=eps)&amp;&amp;(m&lt;=9))
      { p=0.0;
0 B8 T& \# S) t+ C( j# b        for (i=0;i&lt;=n-1;i++)
          { x=a+(i+0.5)*h;
; Q. `% `* Z7 k) D: }5 w9 O            p=p+frombf(x);
4 ^+ Q' P4 p# i# U! ?; B9 R2 S0 |5 F          }
# A! k5 A% n9 @4 `        p=(y[0]+h*p)/2.0;
5 a" \" U  z; X+ b" O        s=1.0;
        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;
& O; p2 e' y7 T# y6 v5 m          }
0 f8 M. r; G' J$ b8 V" v" e1 W7 u        ep=fabs(q-y[m-1]);
" F& x  k# @/ d7 `0 P! T- @        m=m+1; y[m-1]=q; n=n+n; h=h/2.0;
      }
+ m( r9 c/ a6 H1 N2 {    return(q);
! i3 G1 b! j+ r  }</P><>呵呵 希望对你有用！！</P>

西西

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

陈了

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