/* Trapezoidal rule */
#include <stdio.h>

double f(double x)
 {return 4.0/(1.0+x*x);}

int main()
{
 int i, n ;
 double a=0.0, b=1.0 , h, s=0.0, x;

 printf("Enter number of partitions = ");
 scanf("%d", &n) ;

 h = (b-a)/n ;

 for (i=1;i<=n-1;i++) s = s + f(a + i*h);

  s=h/2*(f(a)+f(b))+ h* s;

  printf("%20.12f\n", s) ;
  return 0;
}

y=[6.28345, 3.22246, 5.39914, 3.60865, 3.65644, 4.9384, 6.71493, 4.66608, 4.95818, 4.87514, 3.65405, 5.58718, 4.40297, 3.56216, 4.26301, 7.17059, 3.8773, 3.65189, 2.62594, 6.38777, 5.10231, 4.97363, 4.92931, 6.88091, 4.02733, 4.91927, 5.35723, 5.22676, 4.93593, 3.10154, 3.80014, 4.06031, 4.0938, 2.32529, 2.33921, 2.54891, 1.34026, 2.67622, 0.942549, 4.21709, 2.67793, 1.7866, 3.44514, 0.0870244, -0.296949, -0.0974269, 3.39252, -0.730539, 1.95505, 2.28344, 0.691118, 1.63496, 3.37392, -0.796132, 3.18227, 2.27203, 2.54027, 1.34103, -0.0650801, 2.77979, 0.714989, -0.835475, 0.531536, 1.98343, -0.335429, 1.24667, 0.81676, 1.38371, 0.667464, 2.42286, 1.00915, -0.0767671, 2.68488, -0.8332, 2.4854, 1.54572, 1.3195, 4.73976, 0.990827, 1.99647, 3.02067, 2.3319, 2.73444, 1.29938, 4.91075, 1.16644, 1.53354, 4.70443, 3.47393, 3.30089, 3.79532, 2.69182, 4.33596, 1.89338, 3.20495, 2.26482, 4.02946, 4.99474, 2.44974, 2.19765, 3.03248, 2.99834, 5.27905, 4.484, 1.15131, 4.72466, 4.10129, 4.63075, 2.25385, 3.92437, 4.09418, 3.91403, 1.342, 1.18493, 2.21738, 2.8371, 1.28084, 4.20303, 3.24356, 3.70767, 2.99487, 1.75703, 2.94221, 3.75199, 1.57677, 2.51171, 5.17309, 1.45874, 3.7486, 2.98649, 4.31283, 3.06652, 2.98886, 1.06716, 3.79495, 2.95823, 3.86032, 2.09225, 2.44123, 2.30055, 4.37755, 3.12176, 3.90007, 1.82306, 2.71626, 4.70279, 2.9145, 1.78829, 5.05247, 5.75262, 4.72955, 5.63518, 4.70264, 2.93391, 5.29481, 5.67629, 3.21819, 3.58788, 6.99283, 6.66301, 6.87711, 1.91899, 5.12603, 3.32968, 5.6337, 3.56705, 3.56726, 3.27185, 5.77811, 7.05581, 6.19436, 1.93323, 6.78286, 5.59857, 5.15317, 5.77036, 5.11903, 3.79379, 5.02595, 3.97256, 4.76766, 4.2986, 5.9101, 3.75019, 4.45336, 3.41256, 2.3897, 4.65437, 4.19994, 1.80178, 1.92235, 2.58043, 1.8049, 0.727705, 2.23352, 3.20678, 1.99824, -0.100039, -0.398669, 1.0254, 1.06878, 0.779708, 0.592887, 1.63836, 0.022458, 2.04382, 0.924611, 1.87949, -0.219661, 2.75017, 2.97187, 0.445748, 0.558413, 0.162202, 2.50927, 1.02211, 2.84669, 3.30699, 0.326872, 3.81569, 0.825398, 1.34227, 2.13528, 3.29024, -0.0106859, 3.46575, 1.8493, 2.85042, 1.55254, 2.14178, 3.87732, 2.93549, 3.72309, 2.3297, 0.492662, 2.01528, 3.55222, 3.57302, 2.80209, 3.99803, 2.50128, 3.18386, 2.12652, 3.82451, 3.95133, 1.21219, 3.24286, 4.94487, 3.599, 3.23689, 0.814642, 3.30503, 4.32353, 2.1707, 2.38114, 2.56006];