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 */
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; }