Абсолютный ассемблер. Вычисление факториала - Assembler

Формулировка задачи:

Добрый вечер, помогите пожалуйста написать программу для вычисления факториала для одно-адресной вычислительной машины на абсолютном ассемблере. У меня имеется пример программы(для вычисления A=B+C*D ) на нем:

       BEGIN 1000 ;Задает начальный адрес программы
FIRST  LDA   KC ;Загрузка в регистр A<-C
       MUL   KD    ;A<-A*D
       ADD   KB    ;A<-A+B
       STA   PA    ;Вывод результата в переменную PA
       HLT           ;Остановка процессора
PA     RESW   1 ;Резервирует заданное количество слов для данных   
KB     WORD   1   
KC     WORD   2  ;Формирует константу занимающее 1 слово
KD     WORD   3
       END FIRST ;Задает первую исполняемую  команду программы

Решение задачи: «Абсолютный ассемблер. Вычисление факториала»

Листинг программы

lda 40; загружаем аргумент
factorial     ;рассчитываем факториал сорока
STA   PA    ;Вывод результата в переменную PA

Объяснение кода листинга программы

lda 40 - загрузка аргумента (в данном случае число 40) в регистр A
factorial - вычисление факториала числа, находящегося в регистре A
STA PA - результат вычисления факториала записывается в переменную PA