Синус и косинус через Ряд Тейлора на внешнем ассемблере

#include <iostream>
#include <time.h>
#include <iomanip>
#include <math.h>
 
extern "C" double _cdecl Cos(double);
    using namespace std;
double Sin(double X)
{
       double eps = 0.01;
       double S = 0;
       double N =  1.;
       double D   = 2.;
       double edm =-1.;
       double ed  = 1.;
       double P   = X;
 
        double R;   
         int sw;
_asm {
        finit
        fld     X      // ; st(0)=x
        fstp    S       //;s=st(0)x
L1:     fld     N      //  ;n -- st(0)
        fmul    D      //  ;2*n --st(0)
        fst     R      //  ;st(0) -- R
        fmul    R     //  ;4*n*n -- st(0)
        fadd    R
        fst     R     //  ; знаменатель 4*n*n -- R
        fld     X    //   ;x --  st(0)
        fld     X    //   ;x -- st(1)
        fmul
        fld     edm    //   ; -1
        fmul           //  ;-x*x -- st(0)
 
fdiv    R
        fmul    P      //  ; элемент
        fst     P
fadd    S
        fstp    S
        fld     N
        fadd    ed     //  ;n=n+1
        fstp    N
        finit
        fld     P
        fabs            //;|p|
 
fcom    eps     //;сравнение
                        //;Eps>st(0) c3=0 c2=0 c0=1
                        //;Eps<st(0) c3=0 c2=0 c0=0
                        //;Eps=st(0) c3=1 c2=0 c0=0
        fstsw   sw     // ;слово состояния в основной процессор
        fwait          // ;Ожидание выполнения операции
        mov eax,sw
      and     ax,0700h   
// или    sahf  сохраняет  ah->регистр флагов процессора
//и потом          je     L1
 
        cmp     ah,0
        je     L1
        finit
        fld     S
        fld     X
        fsin
        fst     R
}                     // printf ("%lf %lf ",S,R);
return S;
 
}
//double Cos(double x)
//
//{
//  const double Pi= 3.1415926535897932384626433832795;
//  const double e = 0.0000000000000001;
//  double k=2,res=1,dpi=2*Pi,p=1,z=0;
//  _asm
//  {
//      finit
//      fld x
//;приведение х в интервал [-2*pi,2*pi]
//m1:
//      fcom dpi
//      fstsw ax
//      sahf
//      jb m2
//      fsub dpi
//      jmp m1
//m2:
//      fchs
//m3:
//      fcom dpi
//      fstsw ax
//      sahf
//      jb m4
//      fsub dpi
//      jmp m3
//m4:
//      fchs
//      fstp x
//;вычисление cos(x) с помощью ряда Тейлора
//m5:
//      fld p
//      fabs
//      fcom e
//      fstsw ax
//      sahf
//      jb mend
//      finit
//      fld p
//      fchs
//      fmul x
//      fmul x              ;-p*x*x - в st(0)
//      fld z
//      fld1
//      fadd                ;z+1 - в st(0), -p*x*x - в st(1)
//      fld z
//      fadd k              ;z+2 - в st(0), z+1 - st(1), -p*x*x - st(2)
//      fmul                ;(z+1)*(z+2) - st(0), -p*x*x - st(1)
//      fdiv                ;(-p*x*x)/((z+1)*(z+2))
//      fst p               ;записать значение st(0) в переменную p
//      fld res
//      fadd                ;прибавить к результату р
//      fstp res
//      fld z
//      fadd k              ;увеличить z на 2
//      fstp z
//      jmp m5
//mend:
//      nop
//  }
//  return res;
//}
void main()
{
    double a,b,h;
    double x,y,z,m;
    clock_t start,finish;
    int timems, times;
    cout<<"Input Min Value"<<endl;
    cin>>a;
    cout<<"Input Max Value"<<endl;
    cin>>b;
    cout<<"Input Shift"<<endl;
    cin>>h; 
    
    
    start=clock();
    cout.precision(7);
    cout<<"Min Value is "<<a<<" Max Value is "<<b<<" Shift is "<<h<<endl;
    for(double x=a;x<=b;x+=h)
    {
    double k = Sin(x);
    double m; /*= Cos(x);*/
    _asm
    {
        fld x
        call Cos
        fstp m
    }
    z=sqrt(1-k);    
    y=(k/m)/z;
    cout<<setw(10)<<k<<" "<<setw(10)<<m<<" "<<setw(10)<<y<<endl;
    }
    
    
finish=clock(); 
timems=(finish-start);
if (timems/1000)
{
    times=timems/1000;
    timems=timems%1000;
    cout<<"Elapsed "<<times<<" s "<<timems<<" ms."<<endl;
}
else
    cout<<"Elapsed "<<timems<<" ms."<<endl;
system("pause");
}