Вычисление рядя C++

[alonix]

Необходимо вывести на экран таблицу значений функций Y(x) и ее разложения в ряд S(x) для x, изменяющегося от a до b с шагом h=(b-a)/10 - n количество элементов ряда, которые необходимо суммировать
Близость значений во всем диапазоне значений x укатывает на правильность вычисления S(X) и Y(X)
S(x)=1+cos(pi/4)*x/1!+...cos(n*pi/4)*x^n/n!
Y(x)=e^(x*cos(pi/4))*cos(x*sin(pi/4))
Для S(x) нахожу рекуррентную формулу

#include <vcl.h>
#include <iostream.h>
#include <iomanip.h>
#include <conio.h>
#include <math.h>

char* Rus(const char* text);



void main()
{
 double x1, x2, x, h;
 int N;
 cout << Rus("Введите начало интервала A: ");
 cin >> x1;
 cout << Rus("Введите конец интервала B: ");
 cin >> x2;
 cout << Rus("Введите количество элементов ряда N: ");
 cin >> N;
 h=(x2-x1)/10.0;
 x=x1;
 cout << "|-------|---------------|-------------|" << endl;
 cout << "|   x   |      s        |      y      |" << endl;
 cout << "|-------|---------------|-------------|" << endl;
 cout << setiosflags(ios::fixed);
 while (x<x2+0.001)
 {
  double a=1+cos(3.14/4);
  double s=a;
  for(int k=1; k<=N; k++)
   {
    a=a*cos((k+1)*(3.14/4))*x*k/k*cos(k*3.14/4);
    s+=a;
   }
   double y=exp(x*cos(3.14/4))*cos(x*sin(3.14/4));
   cout << setprecision(3) << "|" << setw(7) << x << "|" << setprecision(10)
        << setw(13) << s << "|" << setw(13) << y << "|"  << endl;
        x+=h;
  }
  getch();
}
char bufRus[256];
char* Rus(const char* text)
{
  CharToOem(text,bufRus);
  return bufRus;
}