Вычислить значение функции разложением в ряд Тейлора

double mycos(double x,double eps)
{
    double s,an;
    int n;
    n=0;
    an=1;
    s=0;
    while(fabs(an)>eps)
    {
         s+=an;
         n++;
         an*=-x*x/(2.*n-1.0)/(2.0*n); 
    }
    return s;
}
double mysin(double x,double eps)
{
    double s,an;
    int n;
    n=1;
    an=x;
    s=0;
    while(fabs(an)>eps)
    {
         s+=an;
         n++;
         an*=-x*x/(2.*n-1.0)/(2.0*n-2.0); 
    }
    return s;
}
double mytan(double x,double eps)
{
    return mysin(x,eps)/mycos(x,eps);
}
double myf(double x,double eps)
{
    return mytan(x)+mytan(x/2.0);
}
int main()
{
     double X,E;
     cin>>X>>E;
     cout<<myf(X,E);
     return 0;
}

#include <iostream>
#include <cmath>
using namespace std;
 
double mycos(double x,double eps)
{
    double s,an;
    int n;
    n=0;
    an=1;
    s=0;
    while(fabs(an)>eps)
    {
         s+=an;
         n++;
         an*=-x*x/(2.*n-1.0)/(2.0*n); 
    }
    return s;
}
 
double mysin(double x,double eps)
{
    double s,an;
    int n;
    n=1;
    an=x;
    s=0;
    while(fabs(an)>eps)
    {
         s+=an;
         n++;
         an*=-x*x/(2.*n-1.0)/(2.0*n-2.0); 
    }
    return s;
}
 
double mytan(double S1,double S2)
{
    return S1/S2;
}
 
int main()
{
     double x, eps;
     cin>>x>>eps;
     cout <<"  tan="<<tan(x)<<endl;
     cout<<"mytan="<<mytan(mysin(x,eps),mycos(x,eps))<<endl;
     cout <<"y="<<1/(mytan(mysin(x,eps),mycos(x,eps))+mytan(mysin(x/2,eps),mycos(x/2,eps)))<<endl;
system("pause");
return 0;
}

#include <iostream>
#include <iomanip> 
#define _USE_MATH_DEFINES
#include <math.h>
#include <cstdlib>
using namespace std;
double mycos(double x, double eps)
{
    double s, an;
    int n;
    n = 0;
    an = 1;
    s = 0;
    while (fabs(an)>eps)
    {
        s += an;
        n++;
        an *= -x * x / (2.*n - 1.0) / (2.0*n);
    }
    return s;
}
double mysin(double x, double eps)
{
    double s, an;
    int n;
    n = 1;
    an = x;
    s = 0;
    while (fabs(an)>eps)
    {
        s += an;
        n++;
        an *= -x * x / (2.*n - 1.0) / (2.0*n - 2.0);
    }
    return s;
}
double mytan(double x, double eps)
{
    return mysin(x, eps) / mycos(x, eps);
}
double myf(double x, double eps)
{
    return mytan(x,eps) + mytan(x / 2.0,eps);
}
int main()
{
    const double E = 1e-8;
    for (double x = 0; x < M_PI_2; x += M_PI * 0.05)
    {
        cout << setw(10)<<x << setw(14) << myf(x, E) << setw(14) << tan(x) + tan(x / 2.)<<endl;
    }
    system("pause");
    return 0;
}