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

西西

<>本人紧急求助，希望各位大虾们伸出你们友情之手，帮我一下。</P>
! \: E1 I) X3 p8 w! W<>1。牛顿迭代法</P>
* p; Y. z+ ^* \" `<>1。全主消元法</P>
/ c* y) U8 S/ t' B9 s0 g<>3。改进平方根法</P>
% r, k$ ?0 z* e( _- [! U" R% a<>4。牛顿向前，向后插值法</P>
  m' L/ y; h2 b1 E8 v  i<>5。加权最小二乘法</P>
<>6。龙贝格积分法</P>
<>拜托了。。。。</P>

loveyuan

<>1。牛顿迭代法</P><>
! v0 V8 x( I7 f8 }6 u' B! B  #include "stdio.h"
  #include "math.h"
1 _2 r4 o+ b' {& J# J0 K3 _  int dnewt(x,eps,js)
9 }; w; Z4 k$ x9 n* e7 ^  int js;
8 |9 F$ ~! y8 q% x7 l; a0 W" K9 i  double *x,eps;
- Q9 q" R3 U: e7 K  { extern void dnewtf();
" ], a" i) h* f    int k,l;
    double y[2],d,p,x0,x1;
    l=js; k=1; x0=*x;
' N$ Z) c5 v, S1 b    dnewtf(x0,y);
    d=eps+1.0;
    while ((d&gt;=eps)&amp;&amp;(l!=0))
      { if (fabs(y[1])+1.0==1.0)
$ E, h$ y  O/ d" ~          { printf("err\n"); return(-1);}
        x1=x0-y[0]/y[1];
) ^, g# G3 Y$ A1 V1 S% \' ?% F1 M        dnewtf(x1,y);
0 Y1 S9 k  K" B! F: [9 t% [        d=fabs(x1-x0); p=fabs(y[0]);
9 K7 e1 Q4 M  n        if (p&gt;d) d=p;
        x0=x1; l=l-1;
/ K  L. C, X; F1 _3 R      }
* I& T; O- o1 d* o9 a    *x=x1;
    k=js-l;
& t7 N1 t. i# |$ l0 Y    return(k);
  }</P><>全主消元法</P><>#include "stdlib.h"
# q9 b4 G! ^# @0 S( @% c  #include "stdio.h"
2 y% z0 P* K+ M  w0 }  int acgas(ar,ai,n,br,bi)
9 X; ^6 k5 w& ]. \- R  int n;
" D9 G& v  h6 I" s* S  double ar[],ai[],br[],bi[];
  { int *js,l,k,i,j,is,u,v;
8 @1 L' |3 M7 _1 ?3 v    double p,q,s,d;
, _" j! F9 |/ E0 r$ m9 C    js=malloc(n*sizeof(int));
    for (k=0;k&lt;=n-2;k++)
      { d=0.0;
        for (i=k;i&lt;=n-1;i++)
" `8 @; H# O9 R3 @1 a) {        for (j=k;j&lt;=n-1;j++)
6 c6 U+ [* b6 F3 i+ L& r          { u=i*n+j;
            p=ar*ar+ai*ai;
6 Y' D5 k, F. x            if (p&gt;d) {d=p;js[k]=j;is=i;}
8 _4 p7 r  [9 A" @+ o" ^7 P+ |/ R          }
' E4 S. y' x& Z1 P5 R2 ^$ E        if (d+1.0==1.0)
          { free(js); printf("err**fail\n");
: d1 f- C& l: f& p) A% N            return(0);
          }
, Q- `, r, J  E        if (is!=k)
          { for (j=k;j&lt;=n-1;j++)
: R- V' Y" a* F2 k              { u=k*n+j; v=is*n+j;
                p=ar; ar=ar[v]; ar[v]=p;
0 Z. _6 L% U8 P3 N                p=ai; ai=ai[v]; ai[v]=p;
              }
            p=br[k]; br[k]=br[is]; br[is]=p;
; O$ G0 x  I& o, {% h4 O            p=bi[k]; bi[k]=bi[is]; bi[is]=p;
          }
0 ~" X# O& j, T        if (js[k]!=k)
; |. D0 j7 o3 ]& X8 ?          for (i=0;i&lt;=n-1;i++)
            { u=i*n+k; v=i*n+js[k];
( E7 C( A. N/ e5 k: a& A* a              p=ar; ar=ar[v]; ar[v]=p;
0 u, b* U3 j3 i1 w0 s, `0 k* j              p=ai; ai=ai[v]; ai[v]=p;
! L% s- s0 j( ]5 A0 p            }
' J  v4 G# `6 x2 A. ~) A" W        v=k*n+k;
0 y) L' l  X) ~/ U" k3 g/ E        for (j=k+1;j&lt;=n-1;j++)
9 y$ N9 U5 C( m/ B, |9 k          { u=k*n+j;
" m' X. V0 ]: X" J) j            p=ar*ar[v]; q=-ai*ai[v];
% W* c; k7 B9 N8 x9 A            s=(ar[v]-ai[v])*(ar+ai);
! L/ t* d- m  }; m            ar=(p-q)/d; ai=(s-p-q)/d;
          }
, i0 R7 [( X% Z: J6 A+ w        p=br[k]*ar[v]; q=-bi[k]*ai[v];
        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++)
& o! v5 D" w! _1 P. R& L! f              { v=k*n+j; l=i*n+j;
& I! z, C# M( j% S0 ~% {% o3 }# q                p=ar*ar[v]; q=ai*ai[v];
                s=(ar+ai)*(ar[v]+ai[v]);
                ar[l]=ar[l]-p+q;
- P9 i$ X4 F# \8 p  h, m; a                ai[l]=ai[l]-s+p+q;
0 H& ?' V0 s% `. {7 s2 F9 y              }
4 p7 |4 S6 W8 J# N            p=ar*br[k]; q=ai*bi[k];
            s=(ar+ai)*(br[k]+bi[k]);
            br=br-p+q; bi=bi-s+p+q;
          }
$ ?3 q/ V7 r7 r& b3 C4 \+ ~% Y' M      }
3 J# @7 T# N/ i; P  K    u=(n-1)*n+n-1;
( a1 D9 Y$ n: @2 p, V+ K    d=ar*ar+ai*ai;
9 p. ]0 J* ~/ q/ a* e- X  @( n    if (d+1.0==1.0)
      { free(js); printf("err**fail\n");
" W1 j% X8 Q. }        return(0);
! a! V. I8 T' X4 z      }
    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;
    for (i=n-2;i&gt;=0;i--)
    for (j=i+1;j&lt;=n-1;j++)
3 S! Q% A+ {' m  U" p% o3 E# Z, }      { u=i*n+j;
7 `7 B( }# N+ b$ B* n% Y) D" K) l6 i* Y        p=ar*br[j]; q=ai*bi[j];
: u% N2 K: z7 \* u        s=(ar+ai)*(br[j]+bi[j]);
        br=br-p+q;
        bi=bi-s+p+q;
& v. i: Q$ [+ K0 ]% E7 f% {      }
    js[n-1]=n-1;
% t# [0 C8 d) k; y7 L4 @3 f" a: a3 F    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;
4 i, _. t) k+ Y3 K9 H/ J        }
    free(js);
    return(1);
  }</P><>平方根法</P><>#include "math.h"
  #include "stdio.h"
  int achol(a,n,m,d)
3 a3 w# |/ Q0 w9 A" w# l6 U; y/ \  int n,m;
* t) Y( ]! A8 Z  double a[],d[];
  { int i,j,k,u,v;
8 R- P7 n% [( q1 r' a/ @    if ((a[0]+1.0==1.0)||(a[0]&lt;0.0))
      { printf("fail\n"); return(-2);}
6 v. G  K' f4 @2 P/ G    a[0]=sqrt(a[0]);
    for (j=1; j&lt;=n-1; j++) a[j]=a[j]/a[0];
+ k7 c% R9 I& \$ T+ W4 r- }    for (i=1; i&lt;=n-1; i++)
" s9 b9 s8 ]8 K6 k      { u=i*n+i;
( U6 S0 V: y9 V* `8 [/ u1 S- A0 Q: z        for (j=1; j&lt;=i; j++)
          { v=(j-1)*n+i;
! r% T' ~! K) N# E# X+ O            a=a-a[v]*a[v];
          }
        if ((a+1.0==1.0)||(a&lt;0.0))
3 W0 e& Q( O- v, f1 ?5 _7 p          { printf("fail\n"); return(-2);}
        a=sqrt(a);
        if (i!=(n-1))
: h) s$ f4 F' g! f, t, y" e% J1 K5 E          { for (j=i+1; j&lt;=n-1; j++)
( R9 u% v% D3 v8 K; t5 }. d* T5 N              { v=i*n+j;
                for (k=1; k&lt;=i; k++)
                  a[v]=a[v]-a[(k-1)*n+i]*a[(k-1)*n+j];
, B; P9 I* \# Y7 h. s" g+ Q+ H2 h                a[v]=a[v]/a;
( e/ Z. W; j- ]8 O4 i# i% k/ S              }
          }
3 J5 H  V  X; P6 P      }
    for (j=0; j&lt;=m-1; j++)
      { d[j]=d[j]/a[0];
& v$ i; W5 N8 ?+ D) M        for (i=1; i&lt;=n-1; i++)
          { u=i*n+i; v=i*m+j;
  q) {$ k2 n$ `( D! p1 c$ g            for (k=1; k&lt;=i; k++)
' l# l! p/ k7 X, Z              d[v]=d[v]-a[(k-1)*n+i]*d[(k-1)*m+j];
! E8 V/ E; M- ^6 \& U8 [  o            d[v]=d[v]/a;
          }
      }
7 _- L" m" Y% F* ^( {9 Y    for (j=0; j&lt;=m-1; j++)
      { u=(n-1)*m+j;
        d=d/a[n*n-1];
8 ^  B" x; U8 W+ s8 b7 \( I        for (k=n-1; k&gt;=1; k--)
          { u=(k-1)*m+j;
3 ]2 \: Y+ N; O' M# @4 O# k            for (i=k; i&lt;=n-1; i++)
: ?7 |: X# D$ `1 H              { v=(k-1)*n+i;
                d=d-a[v]*d[i*m+j];
* Q  Z& v9 l9 p7 L1 ?1 x; I              }
            v=(k-1)*n+k-1;
7 x4 U% b8 V4 I9 U            d=d/a[v];
          }
- |; h" k% c$ ~9 o6 o2 I. r      }
4 d+ }% v, ~: |  n/ p3 q: ]    return(2);
" k4 _/ w5 {0 z+ P% u2 a  }</P><>牛顿向前，向后插值法</P><>double eelgr(x0,h,n,y,t)
# ]) z) c& |" ~4 [& e6 Q: a  int n;
  double x0,h,t,y[];
1 B  u  ]# b) W' O5 t  { int i,j,k,m;
    double z,s,xi,xj;
% O" m7 A7 I. g  `; d( z) W! Q5 q    float p,q;
# ~1 ^1 x: t! K  e% a7 s    z=0.0;
    if (n&lt;1) return(z);
    if (n==1) { z=y[0]; return(z);}
    if (n==2)
2 g" N3 ^7 ]+ X% o4 l      { z=(y[1]*(t-x0)-y[0]*(t-x0-h))/h;
; `0 M: S1 R$ C8 p) i" |        return(z);
      }
) T, i& L0 W) @9 F/ l    if (t&gt;x0)
; j- X! L; h+ b. m. y$ t# d      { p=(t-x0)/h; i=(int)p; q=(float)i;
9 V; Y* d8 U3 g! f8 g8 T) T        if (p&gt;q) i=i+1;
& _$ x) q  L+ X7 E7 P+ X1 o      }
    else i=0;
    k=i-4;
: Y6 A7 q5 |8 L# r    if (k&lt;0) k=0;
; v- n- p% s) g8 s7 M    m=i+3;
4 d3 K0 |1 k5 F' O% Q    if (m&gt;n-1) m=n-1;
( p. {$ `5 V; q0 J! J    for (i=k;i&lt;=m;i++)
      { s=1.0; xi=x0+i*h;
- S4 c6 Y8 u6 n        for (j=k; j&lt;=m; j++)
- d6 l: e: K2 u& o, A" z) X          if (j!=i)
; [2 g9 T8 C6 v3 S% @            { xj=x0+j*h;
              s=s*(t-xj)/(xi-xj);
            }
        z=z+s*y;
      }
    return(z);
  }
" B' `9 Q  }( ^向前，向后是一样的思想！！</P><>加权最小二乘法</P><>#include "math.h"
( i8 ^8 x+ o* E9 _" u$ S  void hpir1(x,y,n,a,m,dt)
7 Q2 w0 R& @. l: h! U) r  int n,m;
  double x[],y[],a[],dt[];
  { int i,j,k;
6 @# O: o* I1 }, b0 X# }    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;
7 q9 V' v( y8 R# C3 \" \" x    if (m&gt;n) m=n;
" |  D# m6 N; B/ _' _1 e& b    if (m&gt;20) m=20;
8 K0 j1 Z4 i5 |$ J6 T& |6 w; a    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;
    for (i=0; i&lt;=n-1; i++)
      { p=p+(x-z); c=c+y;}
    c=c/d1; p=p/d1;
    a[0]=c*b[0];
    if (m&gt;1)
      { t[1]=1.0; t[0]=-p;
        d2=0.0; c=0.0; g=0.0;
# H5 M9 m* I- {5 J' B! G        for (i=0; i&lt;=n-1; i++)
: q8 W  ?% L0 n1 g9 o" q) M          { q=x-z-p; d2=d2+q*q;
3 r3 n% c/ ?6 V+ ^& E# b: U            c=c+y*q;
            g=g+(x-z)*q*q;
7 a6 s9 K* C& S0 r- i          }
# b" o5 k& \# z) L        c=c/d2; p=g/d2; q=d2/d1;
        d1=d2;
2 _% x) K; Y+ M- Q8 v( \        a[1]=c*t[1]; a[0]=c*t[0]+a[0];
. B: s; {3 X# X- ?      }
* A3 m" z, r. N$ I9 t    for (j=2; j&lt;=m-1; j++)
# ^2 ~2 B( d! q' L3 A5 X      { s[j]=t[j-1];
0 n0 Q* L* Q7 s5 R+ S  p        s[j-1]=-p*t[j-1]+t[j-2];
) _& j% O; Q& T+ H        if (j&gt;=3)
6 ^9 D' t2 ^, w7 [9 T0 ?          for (k=j-2; k&gt;=1; k--)
            s[k]=-p*t[k]+t[k-1]-q*b[k];
, ]) \" X0 h. g        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++)
          { q=s[j];
            for (k=j-1; k&gt;=0; k--)
( G$ g" ~4 I) _              q=q*(x-z)+s[k];
            d2=d2+q*q; c=c+y*q;
' j# K6 w  {& S* t, k3 }            g=g+(x-z)*q*q;
          }
, W$ K  R6 q( v6 S$ Q2 [/ T; w        c=c/d2; p=g/d2; q=d2/d1;
        d1=d2;
        a[j]=c*s[j]; t[j]=s[j];
- Y5 n- ]: p& y        for (k=j-1; k&gt;=0; k--)
& A4 @5 P: W7 \          { a[k]=c*s[k]+a[k];
            b[k]=t[k]; t[k]=s[k];
6 G: F+ H/ G# T4 C: D+ N          }
% ~  ~7 }/ \. R! f" A      }
5 d# C' _9 a+ Q% e9 u    dt[0]=0.0; dt[1]=0.0; dt[2]=0.0;
+ l( B( J( r8 ]. f9 t$ ]2 I5 {    for (i=0; i&lt;=n-1; i++)
2 S* N$ I4 a2 F- n      { q=a[m-1];
        for (k=m-2; k&gt;=0; k--)
          q=a[k]+q*(x-z);
        p=q-y;
        if (fabs(p)&gt;dt[2]) dt[2]=fabs(p);
        dt[0]=dt[0]+p*p;
        dt[1]=dt[1]+fabs(p);
8 ^% r$ t8 F& Y$ o      }
. l* A& {- U& x. p7 u3 Q1 u% m    return;
/ A" o' W3 k- M( r; {5 F: U  }</P><>龙贝格积分法</P><>#include "math.h"
  double fromb(a,b,eps)
! Y# T+ V( j2 J' c5 s0 m3 H  double a,b,eps;
  { extern double frombf();
    int m,n,i,k;
  d7 x* z8 u* H) H    double y[10],h,ep,p,x,s,q;
2 s$ h' r7 E  o0 z    h=b-a;
' B" }. \. j6 i! R7 r" Y5 q    y[0]=h*(frombf(a)+frombf(b))/2.0;
" u, \& ^- s/ Z! p! d    m=1; n=1; ep=eps+1.0;
    while ((ep&gt;=eps)&amp;&amp;(m&lt;=9))
( b- c  f  c$ M2 H; ]8 V3 W      { p=0.0;
; S3 e8 p2 r6 N! y3 g  `        for (i=0;i&lt;=n-1;i++)
          { x=a+(i+0.5)*h;
* W7 T: G% @9 |; u7 S5 e" [" A            p=p+frombf(x);
          }
        p=(y[0]+h*p)/2.0;
        s=1.0;
+ Q& D' t3 l( l9 l        for (k=1;k&lt;=m;k++)
          { s=4.0*s;
            q=(s*p-y[k-1])/(s-1.0);
2 \5 Z, j2 y) t/ a) w& r            y[k-1]=p; p=q;
          }
# C! h% D/ z/ X' i' l8 c        ep=fabs(q-y[m-1]);
        m=m+1; y[m-1]=q; n=n+n; h=h/2.0;
      }
    return(q);
8 d4 c* @! _8 P3 b! X* h  }</P><>呵呵 希望对你有用！！</P>

西西

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

陈了

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