Проблемы с реализацией программы

#include "stdafx.h"
#include<iostream>
#include<cmath>
 
using namespace std;
 
#define min 0.1//минимальное значение
#define max 0.2// максимальное значение
#define step 0.01 // шаг
 
double Sin (double x, double eps)// приближенное вычисление синуса
{
    double s = 1;
    double Rn = 1;
    double k = 1;
    do
    {
        Rn = Rn*(((-1)*x*x*x)/(2*k*(k+2)));
        s = s+Rn;
        k++;
    }
    while(fabs(Rn)>=eps);
 
    return s;
}
 
double Sh(double x, double eps)
{
    double s = 1;
    double Rn = 1;
    double k = 1;
    do
    {
        Rn = Rn*((x*x*x)/(2*k*(k+2)));
        s = s+Rn;
        k++;
    }
    while(fabs(Rn)>=eps);
    return s;
}
double square_root(double x, double eps)
{
    double f1 = 1.5;
    double f2 = 0;
    while (fabs(f1 - f2)>eps)
    {
        f2 = f1;
        f1 = 0.5*(f2+x/f2);
    }
    return f1;
 
}
 
int _tmain(int argc, _TCHAR* argv[])
{
    
    double eps = 0;
    printf("eps = ");
    scanf("%lf", &eps);
    double epsPhi = eps/(1.5);
    double epsG = eps/(1.08);
    double epsF = eps/3;
    int n = 1 + (int)((max-min)/step); // вычисление колва шагов
    double *a = new double[n];// массив для приближенных значений z
    double *b = new double[n]; // массив для точных значений z
    int k = 0;
    for (double x = min; x <= max; x = x+step)
    {
        a[k] = square_root(Sin(x+0.74,epsG),epsF)*Sh((0.8*(x*x)+0.1),epsPhi);
        b[k] = sqrt(sin(x+0.74))*sinh(0.8*pow(x,2)+0.1);
        printf("Z*(%.3f) = %.10f\n", x, a[k]);// приближенные значения
        printf("Z(%.3f) = %.10f\n", x, b[k]);//точные значения
        printf("Abs(Z-Z*) = %e\n\n", fabs(a[k]-b[k])); // вычисление и вывод погрешности
        k++;
 
    }
 
    return 0;
}