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

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

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

  sum=0;
  for(i=0;i<n;i++)  // 分割数だけ回す
    {
      x=d+i*dx;
      fx=x*sin(x);
      sum=sum+fx*dx;
      if(i%100==0)
	{
	  printf("%4d %8.5f %8.5f %8.5f\r",i,x,fx,sum);
	}
    }
  printf("数値積分結果 %15.13f  (真値 %15.13f)\n",sum,M_PI);
  return 0;
}