内部RAM 20H单元有两个非零的BCD数，请编写求两个BCD数的积并把积送入21H单元的程序

在8051单片机中，处理BCD码需要特殊的指令来转换和操作。首先，我们需要读取两个非零的BCD数存储在20H和21H单元，然后将它们相乘，并将结果保存到21H单元。这里是一个简单的示例程序，假设我们有`MOV A, #20H` 和 `MOV R1, #21H` 来准备源地址。

; 将20H单元内容转为十进制然后左移一位
MOVC A @A+PC       ; 指针自动加一，相当于 MOV A, BCD_20H
ADD A, #1         ; 左移一位用于乘法

; 循环乘法，因为BCD乘法通常采用二进制表示先按位相乘
MUL AB            ; A=A*BCD_20H
MOV R2, A        ; 结果暂存R2

; 移动到下一个BCD数
INC A             ; 指向下一个BCD数字
DJNZ A, MUL_LOOP  ; 当A变为0时跳出循环，否则继续

; 把结果存回21H单元，注意需要先清零21H单元
MOV R1, #21H      ; 目标地址
CLR C              ; 清除低4位（低四位可能会溢出）
MOV B, R2          ; 结果的高4位
ADC A, #0          ; 进位加0，确保高位无误地存入

; 更新地址并存入
MOVX @R1, A       ; 存储结果

MUL_LOOP:
    SJMP MUL_LOOP   ; 返回循环开始

; BCD_20H存放着待乘的第一个BCD数

请注意，以上代码是简化的示例，实际编程中可能需要更复杂的循环控制和错误检查。同时，如果内部RAM不够大，可能需要考虑临时变量或者内存分配。此外，对于8051的具体指令集和寄存器配置，可能需要参考具体的硬件手册进行调整。