pastebin

pastebin

Paste Search Dynamic
Recent pastes
trapezoidalRule
  1. #include <stdio.h>
  2. #include <math.h>
  3.  
  4. double f(double x) {
  5.     return 1.0 / (1.0 + x * x);
  6. }
  7.  
  8. double trapezoidalRule(double a, double b, int n) {
  9.     double h = (b - a) / n;
  10.     double sum = f(a) + f(b);
  11.  
  12.     for (int i = 1; i < n; i++)
  13.     {
  14.         double x = a + i * h;
  15.         sum += 2.0 * f(x);
  16.     }
  17.  
  18.     return (h / 2.0) * sum;
  19. }
  20.  
  21. int main() {
  22.     double a = 0.0; // Lower limit of integration
  23.     double b = 6.0; // Upper limit of integration
  24.     int n = 6; // Number of intervals
  25.  
  26.     double integral = trapezoidalRule(a, b, n);
  27.  
  28.     printf("Approximate integral: %.6fn", integral);
  29.  
  30.     return 0;
  31. }
Parsed in 0.009 seconds