第一题:
#include<stdio.h> #include<math.h>
#include<float.h>//为了调用FLT_EPSILON,防止出现分母为零的情况。//
#define N 1000//保证有足够多的运算次数能达到比较精确的结果,且便于修改。//
#define pre 0.0000000000001//使循环能够在达到某一误差精度后停下来。//
double main() {
double p[N],m,err[N],relerr[N]; int k,n;
p[0]=0.879;//给定一个合理的初始值以此来进行迭代。// for(k=0;k<N;k++) {
m=p[k]-cos(p[k]);
p[k+1]=pow(p[k],m);//进行迭代的方程,由书上给出。// err[k]=fabs(p[k+1]-p[k]);//对绝对误差的运算。// relerr[k]=err[k]/(fabs(p[k+1])+FLT_EPSILON);//对相对误差的运算。//
n=k;//把k赋给n,使下面输出时能够在合适的时候停下来。
百度搜索“77cn”或“免费范文网”即可找到本站免费阅读全部范文。收藏本站方便下次阅读,免费范文网,提供经典小说综合文库非线性方程的数值计算方法实验(10)在线全文阅读。
相关推荐: