示例#1
0
文件: h1pt2.c 项目: tcsiwula/c_code
int main(void) {
    double  area, area2;       /* Store result in area       */
    double  a, b;       /* Left and right endpoints   */
    int     n;          /* Number of trapezoids       */
    double  h;          /* Trapezoid base width       */
    
    printf("Enter a, b, and n\n");
    scanf("%lf", &a);
    scanf("%lf", &b);
    scanf("%d", &n);
    
    h = (b-a)/n;
    area = Trap(a, b, n, h);
    
    if(n%2 == 0){
        area2 = Simp(a, b, n, h);
        printf("With n = %d parabolas, our estimate\n", n);
        printf("of the area from %f to %f = %.15f\n",
               a, b, area2);
        double abs = area - area2;
        
        int num = fabs(abs);
        printf("The absolute vale of the difference between %f and %f is %d.", area, area2, num);
    
       
    }
    printf("With n = %d trapezoids, our estimate\n", n);
    printf("of the area from %f to %f = %.15f\n",
           a, b, area);
    
    main();
    return 0;
}  /* main */
示例#2
0
int main() {
	
	float delta, b, x, N;
	float Res1, nElements;
	
	N = 10;
	delta = 0.1;
	x = 3;	
	
	Res1 = 0.0;
	
	for(b = 0; b < x; b += delta) {
		printf("%f\t%f\n",b,Res1);
		Res1 += Simp(b,b+delta,N);
	}
		
	return 0;
	
}