Посчитать значение функции sin, используя разложение ее в ряд Тейлора (используя MMX)

#include "stdafx.h"
#include <iostream>
#include <stdlib.h>
#include <iomanip>
#include <windows.h>
#include <math.h>
 
 
using namespace std;
double x = 0.056739;
int f = 1;
double res;
int main()
{
    setlocale(LC_ALL, "Russian");
 
    __asm{
        lea esi, x
            lea edi, f
            fld qword ptr[esi]    //загрузить FPU x = 0.056739
            mov ecx, 5;           // счетчик цикла
            fldz;     //загрузить FPU sum = 0(помещаем в стек 0)
        
        cycle:
            fld qword ptr[esi]; //загружаем x = 0.056739
            fmul;               //st(0) = x*x
            fld qword ptr[esi]  //загружаем x 
            fmul                //st(0) = x*x*x
            fchs;               //меняем знак (st(0) = -x*x*x)
 
            fidiv WORD PTR[edi] //st(0) / f (целочисленное деление)
                mov ebx, f          //ebx = f
                inc ebx             //ebx = ebx + 1
                mov f, ebx          //f = f + 1
 
            fidiv WORD PTR[edi]     // st(0) / (f*(f + 1)) = -x ^ 3 / 6
            mov ebx, f              //ebx = f + 1
            inc ebx                 //ebx = ebx + 1
            mov f, ebx              //f = f + 2
            loop cycle;             //крутимся в цикле 5 раз
          fxch                  //в очередной раз меняем знак      (меняет местами содержимое st(0) и st(i))
            fstp res                 //res = x - x ^ 3 / 6 + x ^ 5 / 120 - x ^ 7 / 5040 + ...
    }
    cout << "sin:" << res << " Встроенный:" << sin(x) << '\n';
    cin.get();
 
    return(0);
}