Question

菜鸟C语言求大神改错（定积分基本数值算法）

a=0 b=3.1415926 N=1000问题1：运行程序时后两种没有输入a、b、N的提示问题2：第三个函数的运行值是2.7的一半左右，为什么不是2.7左右？程序哪里写错了？谢谢大神的指导！！！#include <stdio.h>#include <stdlib.h>#include <math.h>void main(){ void I1(); void I2(); void I3(); I1(); I2(); I3();}void I1(){ double f(double x); int N,i; double a,b,deltax,I1=0; printf("a,b,N\n"); scanf("%lf%lf%d",&a,&b,&N); deltax=(b-a)/N; for(i=0;i<N;i++) I1=I1+deltax*f(a+deltax*i); printf("\n输出I1值\n"); printf("%lf\n",I1);}void I2(){ double f(double x); int N,i; double a,b,deltax,I2=0; printf("a,b,N\n"); scanf("%lf%lf%d",&a,&b,&N); deltax=(b-a)/N; I2=I2+0.5*deltax*(f(a)+f(b)); for(i=1;i<N;i++) I2=I2+deltax*f(a+deltax*i); printf("\n输出I2值\n"); printf("%lf\n",I2);}void I3(){ double f(double x); int N,i; double a,b,deltax,I3=0; printf("a,b,N\n"); scanf("%lf%lf%d",&a,&b,&N); deltax=(b-a)/N; I3=I3+1/3*deltax*(f(a)+f(b)); for(i=1;(2*i-1)<N;i++) I3=I3+4/3*deltax*f(a+deltax*(2*i-1)); for(i=1;(2*i)<N;i++) I3=I3+2/3*deltax*f(a+deltax*(2*i)); printf("\n输出I3值\n"); printf("%lf\n",I3);}double f(double x){ double y; y=cos(1/(1+x*x)); return (y);}

Answer 1

#include <stdio.h>
#include <stdlib.h>
#include <math.h>
 main()
{ void I1();
 void I2();
 void I3();
 I1();
 I2();
 I3();
}
void I1()
{ double f(double x);
 int N,i;
 double a,b,deltax,I1=0;
 printf("a,b,N\n");
 scanf("%lf%lf%d",&a,&b,&N);
 deltax=(b-a)/N;
 for(i=0;i<N;i++)
 I1=I1+deltax*f(a+deltax*i);
 printf("\n输出I1值\n");
 printf("%lf\n",I1);
}
void I2()
{ double f(double x);
 int N,i;
 double a,b,deltax,I2=0;
 printf("a,b,N\n");
 scanf("%lf%lf%d",&a,&b,&N);
 deltax=(b-a)/N;
 I2=I2+0.5*deltax*(f(a)+f(b));
 for(i=1;i<N;i++)
 I2=I2+deltax*f(a+deltax*i);
 printf("\n输出I2值\n");
 printf("%lf\n",I2);
}
void I3()
{ double f(double x);
 int N,i;
 double a,b,deltax,I3=0;
 printf("a,b,N\n");
 scanf("%lf%lf%d",&a,&b,&N);
 deltax=(b-a)/N;
 I3=I3+1.0/3*deltax*(f(a)+f(b));            //有改动
 for(i=1;(2*i-1)<N;i++) 
   I3=I3+4.0/3*deltax*f(a+deltax*(2*i-1));            //有改动
 for(i=1;(2*i)<N;i++)
   I3=I3+2.0/3*deltax*f(a+deltax*(2*i));            //有改动
 printf("\n输出I3值\n");
 printf("%lf\n",I3);
}
double f(double x)
{ double y;
 y=cos(1/(1+x*x));
 return (y);
}

追问

大神，为啥加个".0"就好了啊？？？

追答

两个整数的商，仍然是一个整数。所以，3/4=0而3.0/4=0.75