| 查看: 698 | 回复: 0 | ||
[求助]
积分编程运行结果出现错误y1和y2=-1.#IND00是怎么回事
|
|
#include<stdio.h> #include <math.h> #define PI 3.14159 double b1,K; double integral(double(*fun)(double x),double a,double b,int n) { double s,h,z; int i; s=(fun(a)+fun(b))/2; h=(b-a)/n; /*积分步长*/ for(i=1;i<n;i++) s=s+fun(a+i*h); z=s*h; return z;/*返回积分值*/ } double f(double x) { return((K-b1*(pow(K,2)-pow(x,2))/2)/sqrt(pow(x,2)-pow((K-b1*(pow(K,2)-pow(x,2))/2),2))); /*修改此处可以改变被积函数*/ } main() { double y,y1,L,K,jd; printf("请分别输入L和K的值:\n"  ;scanf("%1f %1f",&L,&K); b1=4/(2*K+1); while(fabs(y1)>0.01) { b1=b1+0.001; y=integral(f,0.5,K,10000);/*修改此处可以改变积分上下限和步长*/ y1=L-2*y; jd=acos(K*(2-b1*K)+b1/4)*180/PI; } printf("b1=%f\n y=%lf\n y1=%lf\n jd=%lf\n",b1,y,y1,jd); } 看网上有人说是非法运算溢出什么的，但是我把b1的初值再加上0.001就没有为0的情况了，结果还是出现那种情况，求高手指点 |
回复此楼» 猜你喜欢
售一区SCI文章T0P，我:8O.551.O54,科目全,可十急
已经有6人回复
-
售一区SCI文章T0P，我:8O.551.O54,科目全,可十急
已经有6人回复
-
售SCI一区文章，我:8.O.551.O.5.4,科目全,可伽急
已经有7人回复
-
售SCI-T0P文章，我:8O.5.5.1.O.54,科目齐全,可+急
已经有11人回复
-
售SCI一区T0P文章，我:8.O.55.1.O54,科目全,可伽急
已经有11人回复
-
售SCI一区T0P文章，我:8O.55.1.O.5.4,科目齐全,可+急
已经有13人回复
-
售SCI一区T0P文章，我:8.O.55.1.O54,科目全,可伽急
已经有8人回复
-
售SCI一区T0P文章，我:8.O.55.1.O.5.4,科目全,可+急
已经有8人回复
-
售SCI一区T0P文章，我:8.O55.1.O.54,科目全,可十急
已经有14人回复
-
售SCI一区T0P文章，我:8O.55.1.O.5.4,科目齐全,可+急
已经有11人回复
点击这里搜索更多相关资源