非线性方程实根的正割法.doc

山东师范大学数学科学学院实验报告实验课程： 非线性方程（组）的解法 实验项目： 非线性方程实根的正割法姓名 郭新国 学号：200708020244 班级： 二班 专业： 信计 指导教师： 朱爱玲 完成日期： 2010-4-20 实验目的:掌握解非线性方程实根的正割法的上机编程运算.实验内容: 问题1. 用正割法求x3-0.2x2-0.2x-1.2=0在区间 1,1.5内的实根（精确到1e-3）./ 用单点弦法求解 问题1中的方程，不妨设为x0 = 1.5, x1 = 1#include #include #include using namespace std;int main()int k;double x0,fx0,x,xk,fxk,err;double f(double x);k = 0;x0 = 1.5;fx0 = f(x0);xk = 1.0;err = 1.0;while ( err 0.001 )x = xk;fxk = f(xk);xk = xk - fxk*(xk - x0)/(fxk - fx0);err = fabs(x-xk);k+;cout K-th iteration: k endl; / 迭代次数cout err= setprecision(12) err endl; / 迭代序列最后两项之间的差cout x= setprecision(12) x endl; / 方程的解cout f(x)= setprecision(12) f(x) endl; / 方程中的函数f(x)在解x点的值return (0);double f(double x)return ( pow(x,3) - 0.2*pow(x,2) - 0.2*x - 1.2 );运行结果： 问题2. 用双点弦法解方程x=exp(-x)/ 用双点弦法求解 问题2中的方程x = exp(-x)，不妨设f(x) = x*exp(x)-1, 设x0 = 0, x1 = 1/ 1e-5为容许误差主要程序代码：#include #include #include using namespace std;int main() int k; double x,xi,xj,fxi,fxj,err; double f(double x); k = 0; xi = 0.0; xj = 1.0; err = fabs(xi - xj); while ( err 0.00001 ) / 迭代序列顺序 xi,xj,x. fxi = f(xi); fxj = f(xj); x = xj - fxj*(xj - xi)/(fxj - fxi); xi = xj; xj = x; err = fabs(xi - xj); k+; cout K-th iteration: k endl; / 迭代次数 cout err= setprecision(12) err endl; / 迭代序列最后两项之间的差 cout x= setprecision(12) x endl; / 方程的解 cout f(x)= setprecision(12) f(x) endl; / 方程中的函数f(x)在解x点