первая интерполяционная формула ньютона

[Sturvi]

#include <cstdlib>
#include <iostream>
#include <math.h>

using namespace std;

int main(int argc, char *argv[])
{
    float x[6], y[6], x1=1.5;
    int n=3,i;
    for (i=0; i<=5; i++)
        {
              cin>>y[i];
              x[i]=i;
        }
    int  If, l;
    double h1, zn;
    float S=y[0], S1=0, S2=1, ch, h;
    h=x[1]-x[0];
    for (i=1; i<=n; i++)
        {
              ch=y[i]-y[i-1]; // дельта y0
              If*=i; // i!
              h1=pow(h,i); //  h^n
              zn=If*h1;   // 1!*h^n
              S1=ch/zn;  // первая часть
              S2*=x1-x[i-1];  // вычисление второй части
              S+=S1+S2;  // прибовляем в общую сумму
              S1=0;
        }
    cout<<S;
              
    system("PAUSE");
    return EXIT_SUCCESS;
}

Подскажите в чем ошибка? что сделано не так... у меня не выдает нужный ответ

[Sturvi]

S+=S1+S2;

кроме этой строчки...

[sever-42]

попробуй под отладчиком раза 2 пройтись
вроде все так, ктоме S+=S1*S2;
If*=i; // i! тут if = if * i, if чему равно?

[Nazgulko]

if не было проинициализированно перед тем как ты его стал использовать
+ как уже написали надо S+=S1*S2 а у тебя S1+S2