В чем ошибка в вычислении интеграла?

[IVANN98]

Определенный интеграл sin(x) a=0, b=pi/2 (пределы) В ЧЕМ ОШИБКА КОДА??????

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

double f( double x)
{return sin(x);}
int i;
double a=0,b,d,S1,S2,e=0.00001,n, pi=3.14159265;
int main()
{

printf("n=");
scanf("%lf",&n);
/*printf("e=");
scanf("%lf",&e); */
b=pi/2;
d=(b-a)/n;
printf("d=%lf\n",d); //////
for (i=1;i<=n-1;i++)
{
a+=d;
S1+=((f(a)+f(a+d))/2)*d;
}
S2=0;
printf("S1=%lf\n",S1); //////

do
{
S2=S1;
n++;
d=(b-a)/n;
for (i=1;i<=n-1;i++)
{
a+=d;
S1+=((f(a)+f(a+d))/2)*d;
}

}
while(fabs(S1-S2)>e);
printf("d=%lf\n",d);
printf("integral=%lf",S1);

return 0;

}

[WorldMaster]

А в чем ошибка кода?

[IVANN98]

в результате не получается единица

[WorldMaster]

А сколько получается?

[IVANN98]

РЕЗУЛЬТАТ

n=20
d=0.078540
S1=0.996405

[IVANN98]

Результат

n=20
d=0.078540
S1=0.996405

[IVANN98]

точнее вот результат

n=20
d=0.078540
S1=0.996405
integral=1.074870

[olej.tsil]

Сообщение от IVANN98

integral=1.074870

В пределах той точности ... на которую вас хватило.