Sequential Implementation

float Trap(float local_a, float local_b, int local_n , float h) {

float integral; /* Store result in integral */

float x;

int i;

float f(float x); /* function we're integrating */

integral = (f(local_a) + f(local_b))/2.0;

x = local_a;

for (i = 1; i <= local_n-1; i++) {

x = x + h;

integral = integral + f(x);

}

integral = integral*h;

return integral;

} /* Trap */

Previous slide

Next slide

Back to first slide

View graphic version