Trapezoidal Method Of Numerical Integration

 #include<stdio.h>  
 float integrate(float);  
 void main()  
 {  
      float a,b,h,n,i,sum=0;  
      printf("\n\tINTEGRATION THROUGH TRAPEZOIDAL RULE");  
      printf("\n\t====================================");  
      printf("\n\nEnter the upper limit a ");  
      scanf("%f",&a);  
      printf("\n\nEnter the lower limit b ");  
      scanf("%f",&b);  
      printf("\n\nEnter the number of sub-intervals ");  
      scanf("%f",&n);  
      h=(b-a)/n;  
      sum=integrate(a)+integrate(b);  
      for(i=1;i<n;i++)  
           sum=sum+(2*integrate(a+i));  
      sum=sum*(h/2);  
      printf("\nThe value of integral is %f",sum);  
 }  
 float integrate(float n)  
 {  
      return 1/(1+(n*n));  
 }