Вопрос Помогите пожалуйста найти ошибку в коде программы. C++

[danil201317]

Дан ряд S = sum(1 / (k^2 + p * k + q)), необходимо найти частичную сумму n элементов. ПОдскажите пожалуйста где ошибся, ибо в программе при n = 1, выдает, что S = 0, а при n = 2, выдает ответ соответствующий при n = 1.

#include
#include
#include
using namespace std;

int main()
{
setlocale(LC_ALL, "rus");

double S = 0;
double S1 = 0;
float S2 = 0;
float S3 = 0;
float a = 0;
double k = 0;
double q = 0;
double p = 0;
int N = 26;
int i = 1;

q = 0.1 * i + sin(N*i / 5);
p = 2 * sin(N*i / 5);

cout << "q = " << q << endl;
cout << "p = " << p << endl;

for (int n = 1; n < 3; n++) {
cout << "При n = " << n << endl << "S = " << S1 << endl;
S1 += (1 / (pow(n,2) + p * n + q));
}

system("pause");
return 0;
}