Kampuchiec

Доброго времени суток помогите найти ошибку в программе

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




int main()
{
float x, s, f, x_nach, x_kon, delta_x, eps, temp;
int n, fact;

printf("\nVvod ishodnyh dannyh\n");
printf("\nx_nach=");
scanf("%f", &x_nach);
printf("\nx_kon=");
scanf("%f", &x_kon);
printf("\ndelta_x=");
scanf("%f", &delta_x);
printf("\neps=");
scanf("%f", &eps);

printf("\nDlya x=%5.3f do %5.3f s shagom $5.3f i tochnostyi %5.3f\n", x_nach, x_kon, delta_x, eps);

printf("\n x n s f");
printf("\n--------------------------------------");

x = x_nach;

do
{
temp = 1.0;
s = 0;
n = 1;
f = exp(-x);
fact=1;
while (fabs(temp) > eps)
{
temp=pow(-1,n)*(pow(x,n)/fact);
s+=temp;
n++;
fact+=n*(n-1);

}


printf("\n%5.3f %d %0.6f %0.6f", x, n, s, f);

x += delta_x;
}
while (x <= x_kon);

printf("\n--------------------------------------");

return 0;
}

Миниатюры

 

Vladimir_S

Это еще что за ахинея? Разберитесь, для начала, что такое факториал и как его правильно считать.

Попробуйте так:

do
{
s = 1;
n = 1;
f = exp(-x);
fact = 1;
while (fabs(temp) > eps)
{
temp=pow(-1,n)*(pow(x,n)/fact);
s+=temp;
n++;
fact*= n;

}
...

И да, fact лучше задать в формате float, иначе он может вылететь за разрядную сетку.

Kampuchiec

Спасибо за исправление, но что то не так почему то программа дальше n=1 не идет что делать?

Vladimir_S

А какое eps вводите?

Kampuchiec

eps=0.01 (10 символов)

Vladimir_S

Ой, да, забыл спросить - при каком х проверяли?

Kampuchiec

От 1 до 10 с шагом в 1, при смене на начальное 0.1 и конечным 1 и с шагом 0.1, такой же результат

Vladimir_S

Попробуйте так:
1. Задайте х=1.
2. Внутри цикла организуйте вывод на экран промежуточных значений n, temp, s и покажите, что выйдет.

Kampuchiec

Выводит вот такую картину

Миниатюры

 

Vladimir_S

Нет, не то.
Еще раз.
1. Уберите цикл по х, оставьте только х=1.
2. В цикле по n организуйте вывод указанных параметров.
3. Покажите текст программы.

К сожалению, СИ-транслятора у меня нет, поэтому я переписал Ваш код на Паскаль так, чтобы было максимально близко. Вот что получилось (стандартной функции Pow у меня нет, поэтому она введена вручную):
Вот что получилось у меня:

Видно, как s постепенно сходится к значению 1/е = 0.367879, что свидетельствует о правильности алгоритма. Вывод: Вы где-то накосячили с вводом программы. Поэтому и прошу выложить Ваш листинг.
Да, Вам еще следует перед циклом while добавить temp = 1.