Trapezoidal Rule

  Write a function of the following form:

double trapint(double f[],double r[],double dr[],int N)

Return value:

• Trapezoid rule approximation to .

Input:

• f[]: values of for ( ).
• r[]: values of .
• dr[]: values of .
• N: number of integration intervals.