>>341
<p>解答例1</p></H3>

<p>
<table border BgColor="#aaffff"><td><font size="+1"><pre><font color=black>#include <stdio.h>

int main(void)
{
double a,b,sekibun,exact<font color=red>,x,dx</font>;
int n<font color=red>,i</font>;

printf("x^2をaからbまで台形則を用いて数値積分します。\n");
printf("aとbを入力してください。\n");
scanf("%lf %lf",&a,&b);
printf("分割数nを入力してください\n");
scanf("%d",&n);
<font color=red>dx=(b-a)/n;
x=a;
sekibun=0.0;
for(i=0; i<n; i++)
{
sekibun += (x*x+(x+dx)*(x+dx))*dx/2.0;
x += dx;
}</font>

exact=(b*b*b-a*a*a)/3.0;

printf("数値積分値は%20.15lfです。\n",sekibun);
printf("厳密値は%20.15lfです。\n",exact);
printf("誤差は%20.15lfです。\n",exact-sekibun);
return 0;
}