Вычислить факториал для числа в формате двойного слова

.model tiny
.code
.486
org 100h
.startup
mov ax,0789h
mov dx,0FAd7h
 
mov di,offset _str
 
mov cx,4
hi:xor  si,si
   shld si,dx,4
   shl  dx,4
   add  si,offset hextable 
   movsb
loop hi
 
mov cx,4
lo:xor  si,si
   shld si,ax,4
   shl  ax,4
   add  si,offset hextable 
   movsb
loop lo
 
mov dx,offset _str
mov ah,9
int 21h
 
xor ax,ax
int 16h
ret
 
hextable db '0123456789ABCDEF'
_str     db 9 dup('$')
end

.model tiny
.code
.286
org 100h
start:
 
std
mov x,len
sub x,2
lea di,buf
add di,x
mov si,di
mov word ptr [di],1
mov bx,1
lodsw
@1:
mul bx
stosw
@8:
mov cx,dx
mov ax,[si]
mul bx
cmp ax,0
jnz @f
mov [si],ax
add [si],cx
inc bx
dec y
jz @4
lea si,buf
add si,x
@9:
lodsw
cmp ax,0
je @9
mov di,si
add di,2
jmp @1
 
@@:
stosw
add [si],cx
jnc @2
adc dx,0
@2:
sub si,2
jmp @8
 
@4:
cld
lea bx,xtable
lea si,buf
lea di,buf1
mov dx,len+2
mov x,dx
@6:
sub x,2
lodsw
cmp ax,0
jz @6
add x,2
sub si,2
mov cx,16
lodsw
sub x,2
@@:
test ax,0f000h
jnz @5
sub cx,4
rol ax,4
jmp @b
 
@7:
lodsw
sub x,2
jz exit
mov cx,16
@5:
rol ax,4
push ax
and ax,0fh
xlatb
stosb
pop ax
and ax,0fff0h
sub cx,4
jnz @5
jmp @7
 
exit:
lea dx,buf1
mov ah,9
int 21h
xor ax,ax
int 16h
 
ret
buf1 db 200 dup ('$')
buf db 80 dup (?)
len equ $-buf
x dw ?
y db 40
xtable db '0123456789ABCDEF'
end start

.586p
.model flat
.code
include windows.inc
includelib user32.lib
extrn _imp__MessageBoxA@16:dword
n equ 73
;19!=121645100408832000
;20!=2432902008176640000
;21!=51090942171709440000
;22!=1124000727777607680000
;23!=25852016738884976640000
;24!=620448401733239439360000 
;25!=15511210043330985984000000
;26!=403291461126605635584000000
;27!=10888869450418352160768000000
;28!=304888344611713860501504000000
;29!=8841761993739701954543616000000
;30!=2,6525285981219105863630848e+32
;31!=8,22283865417792281772556288e+33
;32!=2,6313083693369353016721801216e+35
;33!=8,68331761881188649551819440128e+36
;34!=2,9523279903960414084761860964352e+38
;35!=1,0333147966386144929666651337523e+40
;36!=3,7199332678990121746799944815084e+41
;37!=1,3763753091226345046315979581581e+43
;38!=5,2302261746660111176000722410007e+44
;39!=2,0397882081197443358640281739903e+46
;40!=8,1591528324789773434561126959612e+47
;41!=3,3452526613163807108170062053441e+49
;42!=1,4050061177528798985431426062445e+51
;43!=6,0415263063373835637355132068514e+52
;44!=2,6582715747884487680436258110146e+54
;45!=1,1962222086548019456196316149566e+56
;46!=5,5026221598120889498503054288003e+57
;47!=2,5862324151116818064296435515361e+59
;48!=1,2413915592536072670862289047373e+61
;49!=6,082818640342675608722521633213e+62
;50!=3,0414093201713378043612608166065e+64
;51!=1,5511187532873822802242430164693e+66
;52!=8,0658175170943878571660636856404e+67
;53!=4,2748832840600255642980137533894e+69
;54!=2,3084369733924138047209274268303e+71
;55!=1,2696403353658275925965100847567e+73
;56!=7,1099858780486345185404564746372e+74
;57!=4,0526919504877216755680601905432e+76
;58!=2,3505613312828785718294749105151e+78
;59!=1,3868311854568983573793901972039e+80
;60!=8,3209871127413901442763411832234e+81
;61!=5,0758021387722479880085681217663e+83
;62!=3,1469973260387937525653122354951e+85
;63!=1,9826083154044400641161467083619e+87
;64!=1,2688693218588416410343338933516e+89
;65!=8,2476505920824706667231703067855e+90
;66!=5,4434493907744306400372924024784e+92
;67!=3,6471110918188685288249859096605e+94
;68!=2,4800355424368305996009904185692e+96
;69!=1,7112245242814131137246833888127e+98
;70!=1,1978571669969891796072783721689e+100
;71!=8,5047858856786231752116764423993e+101
;72!=6,1234458376886086861524070385275e+103
;73!=4,4701154615126843408912571381251e+105
start:  mov edi,offset result+len;edi - указатель на последний байт результата
    mov esi,n
    mov dword ptr [edi-4],esi
    sub edi,8
    mov ebp,1
    dec esi
@1: mov ecx,ebp
    xor ebx,ebx
@@: mov eax,dword ptr [edi+4*ecx]
    mul esi
    add eax,ebx
    adc edx,0
    mov ebx,edx
    mov dword ptr [edi+4*ecx],eax
    loop @b;.untilcxz
    je @2;  .if !ZERO?
    mov dword ptr [edi],edx
    cmp dword ptr [edi+4*ebp],0
    je @f
    inc ebp
@@: sub edi,4
@2: dec esi
    jne @1;.until ZERO?
; big-endian --> litle-endian
    mov edi,offset result
    mov cl,len/4;длина результата в двордах
@@: mov eax,[edi]
    bswap eax
    stosd
    loop @b
; перевожу hex->dec
    mov esi,offset terminator-1;
@3: mov cl,len;длина в байтах значения делимого
    mov edi,offset result
    xor eax,eax
@@: mov al,[edi];делимое
    div ten
    stosb
    loop @b
    or [esi],ah;остаток
    dec esi
    cmp dword ptr [edi-4],0;cmp result+36,0
    jne @3
    cmp dword ptr [edi-8],0;cmp result+32,0
    jne @3
    cmp dword ptr [edi-12],0;cmp result+32,0
    jne @3
; вывожу на экран
    inc esi
    push MB_OK + MB_ICONASTERISK
    push offset mesbox_title   
    push esi;offset mesbox_text с поправкой на ведущие нули
    push 0
    call _imp__MessageBoxA@16
    ret;выходим из программы
result  db 44 dup (0)
len = $-result
ten db 10
mesbox_text db 112 dup ('0')
terminator db 0
mesbox_title db 'Факториал',0
end start

.model tiny
.286
.code
org 100h
.startup
mov si,offset result1+len-2
mov di,offset result2+len-2
mov dx,1
std
calc:mov  bx,dx
     mov  cx,len-2
     _mul:pusha
          _add:lodsb
               add al,[di]
               aaa
               stosb
               adc byte ptr [di],0
          loop _add
          popa
          dec  bx
     jnz _mul
     inc  dx
     xchg si,di
     inc  cx
     pusha
     rep  movsb
     popa 
     cmp  dx,n
jnz calc
 
lea  dx,[si-len+2]
bcdtoasc:add byte ptr [si],'0'
         dec si
loop bcdtoasc
         
mov  ah,9
int  21h
 
xor  ax,ax
int  16h
ret
 
n=29
result1 db 30 dup(0),1,'$'
result2 db 30 dup(0),1,'$'
len = $-result2
end