#include <stdio.h>
#include <math.h>

#define f(x) ((x)*(x))

int main(void)
{
  int i,n;
  double x,fx,d,u,dx,sum;

  n=1000;      // 分割数
  d=0;         // 積分の下限
  u=4;         // 積分の上限
  dx=(u-d)/n;  // 分割の１区間

  sum=0;       // sumはΣの値の計算
  for(i=0;i<n;i++)  // 分割数だけ回す
    {
      x=d+i*dx;
      sum=sum+f(x);
    }
  printf("解析解：        %15.12f\n",(u*u*u-d*d*d)/3);
  printf("数値積分(短冊)：%15.12f\n",sum*dx);
  printf("数値積分(台形)：%15.12f\n",(sum+(f(u)-f(d))/2)*dx);
  return 0;
}