Что не правильно(нельзя пользоваться библиотекой cmath)

#include <iostream>
#include <iomanip>
#define pi 3.1415962
#define e 2.71828
double Mcos(double x, int acc=7);
double Msqrt(double x);
double Mfact(double x);
double Mpow(double x , double y);
double Msin(double x , int acc=7);
double Mtan(double x , int acc=7);
double trig(double x);
double Mabs(double x);
double Mfunc(double x);
double Mmax(double a, double b);
 
int main(){
    setlocale (LC_ALL,"Russian");
    double y , a , b , h;
    std::cout<<"Левая граница а=";
    std::cin>>a;
    while (true){
    std::cout<<"Правая граница b=";
    std::cin>>b;
    if (a < b) break;
    else std::cout<<"Ошибка, правая граница больше левой!";
                }   
    std::cout<<"Введите шаг h=";
    std::cin>>h;
    int size=(int)((b-a)/h)+1;
    double Y[size][2] , x=a;
    for (int i=0 ; x < b ; i++){
        Y[i][0]=x;
        y=Mfunc(x);
        Y[i][1]=y;
        x+=h;   
    }
    double max_y=Y[0][01];
    for (int i=0 ; i<=size ; i++){
        std::cout<<std::setw(5)<<Y[i][0]<<"   "<<std::setw(5)<<Y[i][1]<<'\n';
    if (Mabs(max_y)<Mabs(Y[i][1])) max_y=Y[i][1] ;
    }
    std::cout<<"Максимaльное значение модуля функции="<<max_y<<'\n';
    return 0;   
}
 
double trig(double x)
{
    int t;
    t=x/(2*pi);
    x-=t*(2*pi);
    if (x>pi) x-=2*pi;
    if (x<pi) x+=2*pi;
}
 
double Mfact(double x)
{
 
    double res=1;
    for (int i=1 ; i<=x ; i++){
    res*=i; 
    }
    return res;
}
 
double Mpow(double x, double y ){
    double res=1;
    for (int i=1;i<=y;i++){
    res*=x; 
    }
    return res;
}
 
double Msin(double x, int acc)
{
    double res = 0;
    x = trig(x);
    for(int i = 0; i < acc; i++) res += Mpow((-1), i) * (Mpow(x, 2 * i + 1) / Mfact(2 * i + 1));
    return res;
 
}
 
double Mcos(double x, int acc)
{
    double res = 0;
    x = trig(x);
    for(int i = 0; i < acc; i++) res += Mpow((-1), i) * (Mpow(x, 2 * i) / Mfact(2 * i));
    return res;
}
 
double Mabs(double x){
    if (x<0) x=x*(-1);
    return x;
}
 
double Msqrt(double x)
{
    double res = x / 2;
    double t;      
    do {
        t = res;
        res = (t + (x / t)) / 2;
    } while ((t - res) != 0);     
    return res;
}
 
double Mtan(double x, int acc){
    double res = Msin(x, acc) / Mcos (x, acc);
    return res;
}
 
double Mfunc(double x){
    double y;
    if (Mabs(x-0.5)<1) y=Mmax(x,Mtan(x))+Msqrt(Mabs(x));
    if (Mabs(x-0.5)>=1) y=Mpow(e,x)*Mcos(x);
 
}
 
  double Mmax(double a, double b) {
    return (a > b) ? a : b ;
}