Некорректный подсчет функции (язык C)

В лабе задали подсчитать такую функцию двумя циклами для заданного пользователем N.

Однако при любом N, даже 1 выводиться, как я понял, позитивная бесконечность. В чем может быть ошибка и как ее исправить?

#include <stdio.h>
#include <stdlib.h>
#include <math.h>
int main()
{
double P=1;
double S=0;
int n;
int i;
int j;
printf("Please insert N \n");
scanf("%d", &n);
for(i=1;i<=n;i++)
{
    for(j=1;j<=i;j++)
    {
    S=S+(4*j/3-1);
    }

P=P*((i+sqrt(i))/S);
S=0;
}
printf("%lf \n",P);

return 0;
}

Источник: https://ru.stackoverflow.com/questions/728162/%D0%9D%D0%B5%D0%BA%D0%BE%D1%80%D1%80%D0%B5%D0%BA%D1%82%D0%BD%D1%8B%D0%B9-%D0%BF%D0%BE%D0%B4%D1%81%D1%87%D0%B5%D1%82-%D1%84%D1%83%D0%BD%D0%BA%D1%86%D0%B8%D0%B8-%D1%8F%D0%B7%D1%8B%D0%BA-c