Для функции вычислить ее приближенное значение с точностью ε=10^-3 в точке х=1,5

// задание3.cpp 
//
#include "stdafx.h"
#include <iostream>
#include <iomanip>
#include <math.h>
using namespace std;
 
int _tmain(int argc,_TCHAR* argv[])
{ 
    setlocale(LC_CTYPE,"Russian");
    float eps=0.001;
    //Определяем необходимые переменные х
    //k-номер члена ряда
    //y-переменная для хранения текущего значения суммы ряда.
    //y1-переменная для хранения предыдущего значения суммы ряда.
    //y0-переменная для хранения значения функции вычисленного
    //на основе стандатрных математических функций.
    double x,y,y1,y0;
    int k=0;
 
    //Вводим переменную х
    cout<<"Введите х=";cin>>x;
 
    y=1;//Задаем начальное значение у
 
    //Основной цикл для расчета суммы ряда
    do
    {
        //сохраняем предыдущие значение суммы ряда
        y1=y;
 
        //Вычисляем текущее значение суммы ряда
        k++;
        //kf-переменная для хранения факториала
        int kf=1;
        //Вложенный цикл для расчета факториала
        for(int m=1;m<=(k*2+1);m++)
        kf*=m;
        //накопление суммы ряда
        y+=pow(-1.,k)*(pow(x,2*k+1)/kf);
    }while(fabs(y-y1) > eps);
        //Количество членов ряда
         k=k+1;
 
    //Расчет точного значения для проверки
    y0=sin(x);
    cout<<"Результат расчета с "<<9<<"членами ряда:"<<setw(9)<<setprecision(5)<<y<<endl;
    cout<<"Точное значение ряда:"<<setw(9)<<setprecision(5)<<y0<<endl;
 
    return 0;
}

if(k%2==0)
y+=p;
else y-=p;