Разработка и отладка циклических алгоритмов и программ

#include "stdafx.h"
#include <iostream>
#include <stdio.h>
#include <math.h>
 
using namespace std;
int fackt(int k)
{
    int f=1;
    for(int i=1;i<=k;i++)
    {
        f*=i;
    }
    return f;
}
double Sx(double x,int n)
{
    double sum=0;
    for(int k=0;k<=n;k++)
    {
        sum+=(pow(-1.0,k))*((2*k*k+1)/fackt(2*k))*(pow(x,2*k));
    }
    return sum;
}
 
int main()
{
    setlocale(LC_ALL,"Russian");
    int i,n,ity=0;
    double a,b,h,x;
    cout<<"Введите значение а(с какого значения начинает изменяться х):\n";
    cin>>a;
    cout<<"Введите значение b(до какого значения будет изменяться х):\n";
    cin>>b;
    cout<<"Введите значение h(с каким шагом будет изменяться х):\n";
    cin>>h;
    cout<<"Введите значение n(нужное для вычесления функций):\n";
    cin>>n;
    ity=(b-a)/h;
    double* Y=new double[ity];
    double* S=new double[ity];
    x=a;
    for(i=0;i<=ity;i++)
    {
        Y[i]=(1-x*x/2)*cos(x)-(x/2*sin(x));
        S[i]=Sx(x,n);
        x+=h;
    }
    x=a;
    for(i=0;i<=ity;i++)
    {
        cout<<"x: "<<x<<" Y(x)= "<<Y[i]<<" S(x)= "<<S[i]<<" S(x)+ |Y(x)-S(x)|= "<<S[i]+abs(Y[i]-S[i])<<"\n";
        x+=h;
    }
    system("PAUSE > void");
    return 0;
}