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

#include <vcl.h>
#include <iostream.h>
#include <conio.h>
#include <math.h>
#pragma hdrstop
     
//---------------------------------------------------------------------------
 
#pragma argsused
int main(int argc, char* argv[])
{   int n=0; int k=0;
float an=1;
float a=1;
float x;
 float sn=1;  float se=1;
    float y;
   for(x=1; x<=2; x+=0.1)
{
 while(n<=15)
    {
    an=an*x/(n+1);
    sn+=an;
    n++;
    }                                           //  n=n+1;
 
     while(a>0.0001)
   {
     a=a*x/(k+1);
     se+=a;
     k=k+1;
   }
 
   y=exp(x);
     cout<<"x="<<x<<"   "<<"sn="<<sn<<" "<<"se="<<se<<" "<<"y="<<y<<endl;
   }
 
   getch();
        return 0;
}

#include <iostream>
#include <string.h>
using namespace std;
 
int main()
{
float x = 1.0;
int n=1;//со второго слагаемого
 
float rec_elem = 1.0;
float summa=1.0;//первое слагаемое для n=0 всегда 1
while( n < 15){
rec_elem *= x/ n;
summa += rec_elem;
n++;
}
cout << summa << endl;
 
system("pause");
    return 0;
}

double myexp(double x,double eps)
{
    double an=1;
    double s=0;
    int n=0;
    while(fabs(an)>eps)
    {
         s+=an;
         an*=x/(n+1);
         n++;
    } 
    return s;
}
double myexp(double x,int nmax)
{
    double an=1;
    double s=0;
    int n=0;
    while(n<=nmax)
    {
         s+=an;
         an*=x/(n+1);
         n++;
    } 
    return s;
}
int main()
{   
     double a,b;
     const int k=10;
     cin>>a>>b; 
      double dx=(b-a)/k;
     cout<<"   x  \t  exp(x)\t  myexp(x)\n";
     for(double x=a;x<=b;x+=dx)
     {
         double y1=exp(x);
         double y2=myexp(x,1e-6);
         cout<<x<<"\t "<<y1<<"\t"<<y2<<endl;
     }
     cout<<"   x  \t  exp(x)\t  myexp(x)\n";
     for(double x=a;x<=b;x+=dx)
     {
         double y1=exp(x);
         double y2=myexp(x,5);
         cout<<x<<"\t "<<y1<<"\t"<<y2<<endl;
     }
     getch();
     return 0;
}