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

西西

<>本人紧急求助，希望各位大虾们伸出你们友情之手，帮我一下。</P>
) U) y' a, {( a  {<>1。牛顿迭代法</P>
& S  Q! w! Q( \<>1。全主消元法</P>
" U& A% B! B" y/ `; U! s& o. [<>3。改进平方根法</P>
<>4。牛顿向前，向后插值法</P>
<>5。加权最小二乘法</P>
<>6。龙贝格积分法</P>
<>拜托了。。。。</P>

loveyuan

<>1。牛顿迭代法</P><>
, v* z. f( R6 X! J2 Q0 Q  #include "stdio.h"
  #include "math.h"
( U; q8 O7 A( M. w1 N0 }& Y0 Q  int dnewt(x,eps,js)
- J, [8 S1 B# t, u% i# k+ P/ ^  int js;
% c9 @( Z$ l2 B' K% D( Q( C7 n/ a# R  double *x,eps;
  { extern void dnewtf();
    int k,l;
2 a) p1 A8 ]; L) v    double y[2],d,p,x0,x1;
. _+ z) V3 ~! X* G    l=js; k=1; x0=*x;
    dnewtf(x0,y);
    d=eps+1.0;
    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);
        d=fabs(x1-x0); p=fabs(y[0]);
        if (p&gt;d) d=p;
6 k2 U( Q  E$ \6 Z0 p; ~4 ]        x0=x1; l=l-1;
6 v1 X& ~. P- ?% }5 k. s      }
    *x=x1;
: z: j& J& n5 _. T    k=js-l;
4 j7 B+ x: ~# T* k8 e    return(k);
  }</P><>全主消元法</P><>#include "stdlib.h"
  #include "stdio.h"
, G3 \9 z* T; U  int acgas(ar,ai,n,br,bi)
6 W5 t. @+ n9 I- v/ U+ |+ w* q  int n;
8 Q: d% B+ r  i) V+ ]3 z  double ar[],ai[],br[],bi[];
  { int *js,l,k,i,j,is,u,v;
) S4 h* U" f& D4 z/ S8 ^    double p,q,s,d;
% o  x( D3 e9 v5 J    js=malloc(n*sizeof(int));
: l) C' S6 z" z$ A' _( v    for (k=0;k&lt;=n-2;k++)
      { d=0.0;
5 [2 J7 @: k; N- Z* H8 N' L        for (i=k;i&lt;=n-1;i++)
9 O) l5 n$ g( S1 _) W! a2 Q6 K        for (j=k;j&lt;=n-1;j++)
          { u=i*n+j;
# P' t5 ?" s' Z" b) d            p=ar*ar+ai*ai;
/ J, B5 _% G* H1 I* H            if (p&gt;d) {d=p;js[k]=j;is=i;}
* E1 d$ v) N! [: S          }
2 Y8 C5 W% V5 {- ~# o( \        if (d+1.0==1.0)
          { free(js); printf("err**fail\n");
            return(0);
# P  J9 z2 N, B, j/ @          }
  I! R' l& F  @; P6 G        if (is!=k)
          { for (j=k;j&lt;=n-1;j++)
2 I" {' b. a- b; p              { 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;
' F" @2 |4 }& p% O( m+ Z- j: z. s          }
        if (js[k]!=k)
  h" I" M( \1 F. t          for (i=0;i&lt;=n-1;i++)
3 W! K; ]1 Z. h( N% T" ^  C. K            { u=i*n+k; v=i*n+js[k];
( S) \2 k3 V' H              p=ar; ar=ar[v]; ar[v]=p;
" E. e' b, S: C6 M/ K7 u4 w              p=ai; ai=ai[v]; ai[v]=p;
            }
        v=k*n+k;
: h. j1 N4 M' p$ V! Z! T        for (j=k+1;j&lt;=n-1;j++)
          { u=k*n+j;
            p=ar*ar[v]; q=-ai*ai[v];
0 p# _3 i' }) Y9 X            s=(ar[v]-ai[v])*(ar+ai);
& q" W- G0 V# d1 Y            ar=(p-q)/d; ai=(s-p-q)/d;
          }
        p=br[k]*ar[v]; q=-bi[k]*ai[v];
% A9 P7 _1 a2 _$ U: i        s=(ar[v]-ai[v])*(br[k]+bi[k]);
  r1 Y2 w9 \" ?0 C        br[k]=(p-q)/d; bi[k]=(s-p-q)/d;
        for (i=k+1;i&lt;=n-1;i++)
1 a2 [; l" g# c. Y: Q          { u=i*n+k;
- i: m% U, ?. p! W1 b: X: p            for (j=k+1;j&lt;=n-1;j++)
. {; H% i+ v" c# P              { v=k*n+j; l=i*n+j;
                p=ar*ar[v]; q=ai*ai[v];
6 D+ P  Z+ t1 {1 F                s=(ar+ai)*(ar[v]+ai[v]);
& F$ z2 t. {) k% m8 i5 o                ar[l]=ar[l]-p+q;
                ai[l]=ai[l]-s+p+q;
              }
            p=ar*br[k]; q=ai*bi[k];
            s=(ar+ai)*(br[k]+bi[k]);
2 V, @4 d7 x& ?8 o, h9 H0 L0 o            br=br-p+q; bi=bi-s+p+q;
          }
      }
    u=(n-1)*n+n-1;
    d=ar*ar+ai*ai;
- u. x3 @# n8 N! U    if (d+1.0==1.0)
! e% [! q! W- L4 p( n      { free(js); printf("err**fail\n");
( ^! s% _" A, t% A. W" J$ Z  C        return(0);
: |+ I& [* W$ l      }
    p=ar*br[n-1]; q=-ai*bi[n-1];
% z% j" \9 ?# X- Q4 H0 U    s=(ar-ai)*(br[n-1]+bi[n-1]);
    br[n-1]=(p-q)/d; bi[n-1]=(s-p-q)/d;
/ L  X- o+ q; U: g    for (i=n-2;i&gt;=0;i--)
    for (j=i+1;j&lt;=n-1;j++)
  j+ [1 v( F  `6 W- ?      { u=i*n+j;
        p=ar*br[j]; q=ai*bi[j];
        s=(ar+ai)*(br[j]+bi[j]);
3 s3 u/ h1 X- N        br=br-p+q;
        bi=bi-s+p+q;
! J; e/ s* A. L8 p9 p      }
    js[n-1]=n-1;
  J, i. t$ Z9 M. }/ g8 w$ W    for (k=n-1;k&gt;=0;k--)
      if (js[k]!=k)
        { p=br[k]; br[k]=br[js[k]]; br[js[k]]=p;
9 d' I0 G2 Q( W; M          p=bi[k]; bi[k]=bi[js[k]]; bi[js[k]]=p;
        }
: V- c$ V6 y+ J8 Z- d0 p/ r4 C    free(js);
    return(1);
/ D' `; v9 D/ E) z; c9 g! `  }</P><>平方根法</P><>#include "math.h"
  #include "stdio.h"
' l4 Q+ ?) p! q  int achol(a,n,m,d)
  int n,m;
  double a[],d[];
# g& A% p. I( w( j" k  { int i,j,k,u,v;
    if ((a[0]+1.0==1.0)||(a[0]&lt;0.0))
) [; P+ r, j# ^" Z3 S; j7 {      { printf("fail\n"); return(-2);}
    a[0]=sqrt(a[0]);
; p9 e  g$ I% w+ t/ g4 ?    for (j=1; j&lt;=n-1; j++) a[j]=a[j]/a[0];
: `% H  x/ Z/ Q: P8 A    for (i=1; i&lt;=n-1; i++)
' k: k- d2 C" T+ x' }! }. h3 E# q      { u=i*n+i;
$ q% U/ ~8 \2 l$ n# H        for (j=1; j&lt;=i; j++)
; `7 R- X- a  z( q          { v=(j-1)*n+i;
            a=a-a[v]*a[v];
          }
        if ((a+1.0==1.0)||(a&lt;0.0))
          { printf("fail\n"); return(-2);}
; z8 n. h+ u0 g: q        a=sqrt(a);
2 q# {) ~) @( V0 |( ^3 X* T        if (i!=(n-1))
          { for (j=i+1; j&lt;=n-1; j++)
              { v=i*n+j;
3 Z1 q/ q. g7 V* o                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;
              }
( @, T4 _  g7 f4 w: {$ J          }
+ y9 x: y* D+ M$ [$ P, o      }
7 L0 N' J2 r# [1 b  y  V    for (j=0; j&lt;=m-1; j++)
      { d[j]=d[j]/a[0];
7 T- I) b9 _1 T0 i( v        for (i=1; i&lt;=n-1; i++)
: q# K! m; ?7 {' X0 U' o          { u=i*n+i; v=i*m+j;
            for (k=1; k&lt;=i; k++)
% B" w; z' j. L. t( D' p              d[v]=d[v]-a[(k-1)*n+i]*d[(k-1)*m+j];
            d[v]=d[v]/a;
/ s5 s; S7 @- [# h. K          }
      }
8 B5 n! ]5 I8 @" `0 K% b8 L+ U    for (j=0; j&lt;=m-1; j++)
. |' b  X. N3 z0 V3 M, R" E+ S      { u=(n-1)*m+j;
& l8 s0 M0 o" A# |  Y* x        d=d/a[n*n-1];
6 T$ Q2 Z, \$ |3 f2 }6 k: }' Q        for (k=n-1; k&gt;=1; k--)
6 ]- ^$ A$ [) p  t6 I" f          { u=(k-1)*m+j;
/ C3 U) G9 W# O* d: [# s( }            for (i=k; i&lt;=n-1; i++)
              { v=(k-1)*n+i;
; U5 g  W8 S( X9 ?" Y7 A                d=d-a[v]*d[i*m+j];
              }
* C$ y6 y8 I" r            v=(k-1)*n+k-1;
6 O$ M4 h; `9 b/ Z( c            d=d/a[v];
          }
' W. S2 G6 q9 B" T4 W      }
2 Z& n1 }2 \/ L* e8 F& i    return(2);
  }</P><>牛顿向前，向后插值法</P><>double eelgr(x0,h,n,y,t)
  int n;
& Z  ?, M% V3 o6 n% m9 n  double x0,h,t,y[];
  { int i,j,k,m;
    double z,s,xi,xj;
# Q2 q1 \2 z0 B$ \    float p,q;
, Q0 G# }8 D5 i1 o* G' Z) m    z=0.0;
+ v& _3 n' C, X0 V) d    if (n&lt;1) return(z);
% P% |# Z* L/ e3 G; s# A    if (n==1) { z=y[0]; return(z);}
6 g( I- T7 h& i    if (n==2)
$ _/ c) n5 k0 @' f3 x) e+ ~/ O      { z=(y[1]*(t-x0)-y[0]*(t-x0-h))/h;
4 K! o) e! [( i; ?0 @        return(z);
2 t6 u+ t6 w6 g/ l6 P      }
    if (t&gt;x0)
      { p=(t-x0)/h; i=(int)p; q=(float)i;
" m2 ~8 b& W* x        if (p&gt;q) i=i+1;
      }
. h+ S9 q$ W1 \. m. D/ W7 n& X    else i=0;
    k=i-4;
    if (k&lt;0) k=0;
' ?8 h0 o' e( ^$ _+ Y- X    m=i+3;
# \  s3 I6 u6 W' o  Q- \9 t: f& d8 e    if (m&gt;n-1) m=n-1;
    for (i=k;i&lt;=m;i++)
      { s=1.0; xi=x0+i*h;
        for (j=k; j&lt;=m; j++)
          if (j!=i)
. _" Y  w* j! f8 T7 B1 p            { xj=x0+j*h;
              s=s*(t-xj)/(xi-xj);
+ j2 C3 z6 d. `2 N            }
        z=z+s*y;
      }
3 G( Q6 S3 N. c! i2 V- S    return(z);
  }
. j4 A. \" f5 i+ @) Y8 l: A向前，向后是一样的思想！！</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];
6 {0 d; y, K1 f    for (i=0; i&lt;=m-1; i++) a=0.0;
    if (m&gt;n) m=n;
! v& w# ~( g. P8 \0 b; l' M8 o, }4 x    if (m&gt;20) m=20;
    z=0.0;
- W6 G# T7 L  p; p- i" }    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++)
      { p=p+(x-z); c=c+y;}
* j! }2 h+ j% C7 d% Z    c=c/d1; p=p/d1;
1 A7 R1 M* E- N2 E1 }$ E6 m" I    a[0]=c*b[0];
  V# ]" y3 Q$ U/ Z5 S1 b; k  O    if (m&gt;1)
      { t[1]=1.0; t[0]=-p;
( q6 }. n: Y/ K# u        d2=0.0; c=0.0; g=0.0;
, ?" l2 {: @5 N        for (i=0; i&lt;=n-1; i++)
$ ?4 k7 B* O. Z( V          { q=x-z-p; d2=d2+q*q;
            c=c+y*q;
            g=g+(x-z)*q*q;
) j. p9 [) g( M+ A% V7 g          }
        c=c/d2; p=g/d2; q=d2/d1;
        d1=d2;
        a[1]=c*t[1]; a[0]=c*t[0]+a[0];
      }
. p4 k! S8 I) }& r    for (j=2; j&lt;=m-1; j++)
      { s[j]=t[j-1];
: w& ~7 u6 q1 g4 b9 D- K  J  r        s[j-1]=-p*t[j-1]+t[j-2];
        if (j&gt;=3)
: t- a+ y' d2 K) y          for (k=j-2; k&gt;=1; k--)
            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++)
6 p6 r* I3 n1 d, D9 G1 o7 B          { q=s[j];
8 N0 C2 A) @; d+ K6 V8 f$ |" S            for (k=j-1; k&gt;=0; k--)
- ^: m8 f# Y( {! H: U              q=q*(x-z)+s[k];
9 z  F7 K2 }3 e/ {            d2=d2+q*q; c=c+y*q;
) ~2 x- V, X8 v$ d- p4 \" f  I4 @            g=g+(x-z)*q*q;
$ h" T. o0 S6 A+ x" B6 ^0 M- h          }
        c=c/d2; p=g/d2; q=d2/d1;
1 i" j+ e  j  F" c1 p$ ~+ g2 }        d1=d2;
. F1 ?3 h1 d6 T% a; a6 j) @        a[j]=c*s[j]; t[j]=s[j];
        for (k=j-1; k&gt;=0; k--)
          { a[k]=c*s[k]+a[k];
            b[k]=t[k]; t[k]=s[k];
; W# ]  y3 ~9 G2 H+ [! n3 q2 h) ]          }
; _# T/ u# a( \& p5 Z( \6 `      }
6 K; m# Z$ V! m. q    dt[0]=0.0; dt[1]=0.0; dt[2]=0.0;
) K0 X  L5 H" e    for (i=0; i&lt;=n-1; i++)
8 G7 V6 R( e$ |5 x# l; o& D: E      { q=a[m-1];
/ L) k2 j! G/ y, h: I8 a        for (k=m-2; k&gt;=0; k--)
          q=a[k]+q*(x-z);
) \9 F" N. ^* }, d5 X' f        p=q-y;
        if (fabs(p)&gt;dt[2]) dt[2]=fabs(p);
% v) j9 ~: ?5 L0 R- C' L        dt[0]=dt[0]+p*p;
        dt[1]=dt[1]+fabs(p);
0 v3 @2 w& f" `      }
* `9 d- o6 D9 F$ j# Y, M; ~    return;
3 z( l' a; V: Y8 t  }</P><>龙贝格积分法</P><>#include "math.h"
  double fromb(a,b,eps)
  double a,b,eps;
  { extern double frombf();
    int m,n,i,k;
1 v+ V; W' i$ u) u9 m* H    double y[10],h,ep,p,x,s,q;
    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))
9 T5 V  @$ h) _! m2 U( l7 k      { p=0.0;
        for (i=0;i&lt;=n-1;i++)
/ U+ @+ g" D. i3 `! Y5 U          { x=a+(i+0.5)*h;
            p=p+frombf(x);
" f# h, s+ n3 M9 N1 {# n0 n: t          }
9 s+ t8 d9 @7 W7 d  \. ?" u        p=(y[0]+h*p)/2.0;
. a# q2 `4 l2 ]        s=1.0;
        for (k=1;k&lt;=m;k++)
          { s=4.0*s;
            q=(s*p-y[k-1])/(s-1.0);
' @0 b( }7 D; ~/ v            y[k-1]=p; p=q;
5 m' U: m% K: ^" l          }
        ep=fabs(q-y[m-1]);
  n6 e7 ?: K$ k) k7 r& C0 s7 c        m=m+1; y[m-1]=q; n=n+n; h=h/2.0;
2 H: i, ]1 T+ W/ @7 y      }
    return(q);
/ L" b" s4 i8 W. K" e% W  }</P><>呵呵 希望对你有用！！</P>

西西

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

陈了

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