Предыдущая страница реферата | 2  3  4  5  6  7  8  9  10  11  12 | Следующая страница реферата



/*Открытие файла-источника*/

while((fp=fopen("data3.xls","w"))==NULL)

{

puts("Error!!! Can't open file nInput name of filen");

}

/*Ввод количества итераций*/

printf("Enter number of Itteration N= ");

scanf("%d",&N);

/*Вычисление шага интегрирования*/

h=(a+b)/n;

printf("Step=%.3fn",h);

/*******Вычисление интеграла методом трапеций*******/

for(j=1;j<=N;j++)

{

h=(a+b)/(int(pow(2,j-1))*n);

Sum[j]=0;

for(i=0;i<=(int(pow(2,j-1))*n);i++)

{

x=a+i*h;

if(i==0)

t=1;

else

t=2;

y=t*(h/2)*(sin(2*x));

Sum[j]=Sum[j]+y;

}

if (j>1)

{

coa=(Sum[j]-Sum[j-1])/Sum[j-1];

printf("Criterion of accuracy=%.5f Number of iteration=%dn",coa,j);

fprintf(fp,"%.7ft",coa);

fprintf(fp,"%dtn",j);

}

}

printf("The sum by a method of trapezes=%.7fn",Sum[1]);

fprintf(fp,"The sum by a method of trapezes=%.7fn",Sum[1]);

/*******Вычисление интеграла методом Симпсона*******/

for(j=1;j<=N;j++)

{

h=(a+b)/(int(pow(2,j-1))*n);

Sum[j]=0;

for(i=0;i<=(int(pow(2,j-1))*n);i++)

{

x=a+i*h;

if(i==0||i==n)

t=1;

Рекомендуем скачать другие рефераты по теме: скачать реферат бесплатно без регистрации, менеджмент.

Предыдущая страница реферата | 2  3  4  5  6  7  8  9  10  11  12 | Следующая страница реферата