开发工具:
文件大小: 2kb
下载次数: 0
上传时间: 2012-06-12
详细说明： #include #include #include using namespace std; class scyt {private: int i,n; double s;//插入点的计算值 double X;//输入的插入点 double z1,z2;//起点和终点的一阶导数值 double *xx,*yy,*a,*b,*c,*f,*h;//b用来存主对角线上的值， //a,c用来存下，上次对角线上的值 double *x,*y,*l,*u;//求解三对角方程所需元素 public: void input();//数据输入 void qiudao();//求二阶导数 void sdj();//解三对角方程 void jscrd();//计算插入点的值 ~scyt()//释放内存 {delete []xx,yy,a,b,c,h,x,y,l,u,f;} }; void main() {scyt myscyt; myscyt.input();//读入x,y数据 myscyt.qiudao();//计算二阶导数 myscyt.jscrd();//计算 插入点的值 } void scyt::input()//读入x,y数据 {ifstream fin("ty.dat"); fin>>n; xx=new double[n]; yy=new double[n]; for(i=0;i>xx[i]>>yy[i];} fin>>z1>>z2; fin.close(); } void scyt::qiudao()//求二阶导数 {a=new double[n-1]; b=new double[n]; c=new double[n-1]; f=new double[n]; h=new double[n-1]; for(i=0;i=0;i--) {x[i]=y[i]-u[i]*x[i+1];} cout<>X; if(X!=1000) {for(i=0;i=X) break; } s=pow((xx[i+1]-X),3)*x[i]/(6*h[i])+pow\ ((X-xx[i]),3)*x[i+1]/(6*h[i])+\ (xx[i+1]-X)*(yy[i]-h[i]*h[i]*x[i]/6)/h[i]\ +(X-xx[i])*(yy[i+1]-h[i]*h[i]*x[i+1]/6)/h[i]; cout<<"s["<
(系统自动生成,下载前可以参看下载内容)
下载文件列表
相关说明
- 本站资源为会员上传分享交流与学习,如有侵犯您的权益,请联系我们删除.
- 本站是交换下载平台，提供交流渠道，下载内容来自于网络，除下载问题外，其它问题请自行百度。
- 本站已设置防盗链，请勿用迅雷、QQ旋风等多线程下载软件下载资源，下载后用WinRAR最新版进行解压.
- 如果您发现内容无法下载，请稍后再次尝试；或者到消费记录里找到下载记录反馈给我们.
- 下载后发现下载的内容跟说明不相乎，请到消费记录里找到下载记录反馈给我们，经确认后退回积分.
- 如下载前有疑问，可以通过点击"提供者"的名字，查看对方的联系方式，联系对方咨询.
