Ошибка в программе на С

[Silken Lili]

Дано задание : удовлетворяют ли элементы последовательности данному рекурентному соотношению: с1*a(i-1)+c2*ai+c3*a(i-1)=b
Вроде написала, но есть ошибка, из-за которой программа не запускается.
Помогите, пожалуйста, с ошибкой =__=


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

int main(void)
{
double x, y, z, b, c1, c2, c3, sum;
FILE *fin=NULL;
if ( (fin=fopen("input.txt", "r"))==NULL)
{printf("didnt open"); return -1;}

printf ("write c1, c2, c3 and b\n");
scanf ("%lf%lf%lf%lf", c1, c2, c3, b);

while (fscanf(fin, "%lf%lf%lf", &x, &y, &z)==1)
sum=c1*x+c2*y+c3*z;

if (sum==b)
printf("you are right!");
if (sum!=b)
printf("You are wrong!");

return 0;
}

[BDA]

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

int
main(void)
{
    double a1, a2, b, c1, c2, c3, num;
    FILE *fin=NULL;
    if ((fin = fopen("input.txt", "r")) == NULL) {
        printf("Cant open\n");
        return -1;
    }
    printf("Input c1, c2, c3 and b:\n");
    scanf("%lf%lf%lf%lf", c1, c2, c3, b);  
    if (fscanf(fin, "%lf%lf", &a1, &a2) < 2) {
        printf("Count of elements too small\n");
        return -1;
    }
    while(c1 * a1 + c2 * a2 + c3 * a1 == b && fscanf("%lf", &num) > 0) {
        a1 = a2;
        a2 = num;
    }
    if (c1 * a1 + c2 * a2 + c3 * a1 == b) {
        printf("True\n");
    } else {
        printf("False\n");
    }
    return 0;
}

Не проверял.

[Silken Lili]

мм, при компилировании он показывает "comparing floating point with == and != is unsafe". у меня уже такое не раз было, так и не поняла почему.
и еще, не могли бы вы объяснить вот эту часть :

while(c1 * a1 + c2 * a2 + c3 * a1 == b && fscanf("%lf", &num) > 0) {
a1 = a2;
a2 = num;
}

[ViktorR]

"comparing floating point with == and != is unsafe". у меня уже такое не раз было, так и не поняла почему.

Так понимаю, что сравнение == или != проходит надежно для целых чисел.
Вещественные надо сравнивать с каким то эпсилон, типа:
if (a-b) <= eps.




Как-то так ...

[BDA]

ViktorR прав. Нужно ввести какое-либо eps.

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

#define EPS 0.000001

int
main(void)
{
    double a1, a2, b, c1, c2, c3, num;
    FILE *fin=NULL;
    if ((fin = fopen("input.txt", "r")) == NULL) {
        printf("Cant open\n");
        return -1;
    }
    printf("Input c1, c2, c3 and b:\n");
    scanf("%lf%lf%lf%lf", c1, c2, c3, b);  
    if (fscanf(fin, "%lf%lf", &a1, &a2) < 2) {
        printf("Count of elements too small\n");
        return -1;
    }
    while(fabs(c1 * a1 + c2 * a2 + c3 * a1 - b) < EPS && fscanf("%lf", &num) > 0) {
        a1 = a2;
        a2 = num;
    }
    if (fabs(c1 * a1 + c2 * a2 + c3 * a1 - b) < EPS) {
        printf("True\n");
    } else {
        printf("False\n");
    }
    return 0;
}

Сначала считываем 2 элемента.
Затем заходим в условие while.
while(fabs(c1 * a1 + c2 * a2 + c3 * a1 - b) < EPS && fscanf("%lf", &num) > 0) {
a1 = a2;
a2 = num;
}
Сначала проверяется условие для считанных ранее a1 и a2 (первый и второй элементы последовательности).
Если оно истинно, тогда выполняется вторая часть условия, т.е. считывается еще 1 элемент (если возможно).
Если эти условия выполнены, то элементы "сдвигаются" (первый отбрасывается, второй становится на место первого, а ново считанный на место второго).
Затем цикл повторяется.
Он выполняется, пока есть элементы в файле и пока выполняется заданное условие на последовательность.

[Silken Lili]

спасибо большое, вроде все понятно)