niudun.rar_operation_计算编程

#include <iostream.h> #include <math.h> #include<iomanip.h> #include<stdlib.h> #define MAX 100 #define eps 0.0000001 void niud (int n,double pa[],double x0); void main() { //cout<<setiosflags(ios::fixed); //cout<<setiosflags(ios::scientific); double x0; double *pa; int n; cout<<"输入方程式的阶数："; cin>>n; cout<<endl; cout<<"输入方程式的各项系数："; pa=new double[n+1]; for(int i=0;i<=n;i++) cin>>pa[i]; cout<<"输入x的初值: "; cin>>x0; niud(n,pa,x0); cout<<endl; delete []pa; } void niud (int n,double pa[],double x0) { double neps; double *px,*pb,*pc; pb=new double[n+1]; px=new double[n+1]; pc=new double[n]; px[0]=x0; cout<<"k= "; cout<<"0"<<endl; cout<<"xk="; cout<<px[0]<<endl; for(int k=0;k<MAX;k++)//开始计算 { pb[0]=pa[0]; pc[0]=pa[0]; cout<<"bi= "; cout<<pb[0]<<" "; for(int i=1;i<=n;i++) //计算 bi { pb[i]=pa[i]+px[k]*pb[i-1]; cout<<pb[i]<<" "; } cout<<endl; cout<<"ci= "; cout<<pc[0]<<" "; for(i=1;i<n;i++) //计算 ci { pc[i]=pb[i]+px[k]*pc[i-1]; cout<<pc[i]<<" "; } cout<<endl; neps=pb[n]/pc[n-1];//计算精度 bn/cn-1 cout<<"bn/cn-1= "; cout<<neps<<endl; cout<<endl; cout<<"k="; cout<<k+1<<endl; cout<<"xk= "; px[k+1]=px[k]-neps;//计算下一个x cout<<px[k+1]<<endl; if(fabs(neps)<eps) //判断精度 { free(px); /*delete []pb; delete []pc;*/ break ; } } cout<<endl; }