out <= addr[15:8] * addr[7:0];

  114 * Initialisation des valeurs : X(0)=0, Y(0)=R...
  115
  116          LDA   #0
                              ===== Page 3 =====
 
  117          STA   YFRAC
  118          STA   XFRAC
  119          STA   X          X(0) = 0.
  120          STA   X+1
  121          LDA   R          Y(0) = R.
  122          STA   Y
  123          LDA   R+1
  124          STA   Y+1
  125
  126          JSR   DOTXYC
  127
  128 * Pas d'iteration si (R = 0) ou (R < 2 en COL140).
  129
  130          LDA   R+1
  131          BNE   VARX
  132          LDA   R
  133          BEQ   RETBAS
  134          LDX   GRAFWRK
  135          CPX   #5
  136          BNE   VARX
  137          CMP   #1
  138          BNE   VARX
  139 RETBAS   RTS              Retour au BASIC.
  140
  141 * Si COL140, XYBIS <- Y/2,
  142 * si BW560,  XYBIS <- Y*2,
  143 * sinon:     XYBIS <- Y.
  144
  145 VARX     LDA   Y
  146          STA   XYBIS
  147          LDA   Y+1
  148          STA   XYBIS+1
  149          LDX   GRAFWRK
  150          CPX   #5
  151          BEQ   VARXA
  152          CPX   #9
  153          BNE   VARX1
  154          ASL   XYBIS
  155          ROL   XYBIS+1
  156          JMP   VARX1
  157 VARXA    LSR   XYBIS+1
  158          ROR   XYBIS
  159
  160 * (XFRAC,X) <- (XFRAC,X) + (XYBIS,XYBIS+1).
  161
  162 * Si une carry est generee de cette addition 16 bits,
  163 * alors modification de X+1...
  164
  165 VARX1    LDA   XFRAC
  166          CLC
  167          ADC   XYBIS      Addition poids faibles...
  168          STA   XFRAC
  169          LDA   X
  170          TAX
  171          ADC   XYBIS+1     et forts.
  172          STA   X
  173          BCC   VARX2      Si pas de report,
                              ===== Page 4 =====
 
  174          INC   X+1         sinon modification de X+1
  175          JMP   VARX3       et allumage d'un nouveau point.
  176 VARX2    CPX   X          Y-a-t-il eu modification de X ?
  177          BEQ   VARY       Non.
  178 VARX3    JSR   DOTXYC
  179
  180 * Si COL140, XYBIS <- X*2,
  181 * si BW560,  XYBIS <- X/2,
  182 * sinon:     XYBIS ,- X.
  183
  184 VARY     LDA   X
  185          STA   XYBIS
  186          LDA   X+1
  187          STA   XYBIS+1
  188          LDX   GRAFWRK
  189          CPX   #5
  190          BEQ   VARYA
  191          CPX   #9
  192          BNE   VARY1
  193          LSR   XYBIS+1
  194          ROR   XYBIS
  195          JMP   VARY1
  196 VARYA    ASL   XYBIS
  197          ROL   XYBIS+1
  198
  199 * (Y,YFRAC) <- (Y,YFRAC) - (XYBIS+1,XYBIS).
  200
  201 * Si report apres cette soustraction 16 bits, alors
  202 * modification de Y+1.
  203
  204 VARY1    LDA   YFRAC      Soustraction poids faible.
  205          SEC
  206          SBC   XYBIS
  207          STA   YFRAC
  208          LDA   Y          Soustraction poids fort.
  209          TAY
  210          SBC   XYBIS+1
  211          STA   Y
  212          BCS   VARY2      Si pas de report.
  213          DEC   Y+1        Modification de Y+1
  214          JMP   VARY3      d'ou allumage du point.
  215 VARY2    CPY   Y          Y-a-t-il eu modification de Y ?
  216          BNE   VARY3      Oui: alors allumage.
  217
  218 VARY4    JMP   VARX       Nouveau point.
  219
  220 VARY3    JSR   DOTXYC     Allumage.
  221
  222          LDA   Y          Est-ce le dernier ? (Y=0).
  223          ORA   Y+1
  224          BNE   VARY4      Non: nouveau point.
  225
  226          RTS              Retour au BASIC.
  227
  228 * DOTXYC est la routine qui assure l'allumage des 4 points
  229 * de coordonnees : (CX+X,CY+Y),(CX-X,CY+Y),
  230 *                  (CX-X,CY-Y),(CX+X,CY-Y).
                              ===== Page 5 =====
 
  231
  232 DOTXYC   LDA   X
  233          CLC
  234          ADC   CX
  235          STA   ABS1
  236          PHA
  237          LDA   X+1
  238          ADC   CX+1
  239          STA   ABS1+1
  240          PHA
  241          LDA   Y
  242          CLC
  243          ADC   CY
  244          STA   ORD1
  245          LDA   Y+1
  246          ADC   CY+1
  247          STA   ORD1+1
  248
  249          JSR   DOTCALL    Allumage de (CX+X,CY+Y).
  250
  251          LDA   CX
  252          SEC
  253          SBC   X
  254          STA   ABS1
  255          LDA   CX+1
  256          SBC   X+1
  257          STA   ABS1+1
  258
  259          JSR   DOTCALL    Allumage de (CX-X,CY+Y).
  260
  261          LDA   CY
  262          SEC
  263          SBC   Y
  264          STA   ORD1
  265          LDA   CY+1
  266          SBC   Y+1
  267          STA   ORD1+1
  268
  269          JSR   DOTCALL    Allumage de (CX-X,CY-Y).
  270
  271          PLA
  272          STA   ABS1+1
  273          PLA
  274          STA   ABS1
  275
  276          JMP   DOTCALL    Allumage de (CX+X,CY-Y).

Legend is:
() mutibyte operation (i.e. bind, add, substract, compare, either 16 or 24bits)
[] to designate pixel coordinates

XFRAC <-0
YFRAC <-0
(X, X+1) <- (0, 0)
(Y, Y+1 <- (R, R+1)

]LOOP:
(XYBIS, XYBIS+1) <- (Y,Y+1)
(XFRAC, X, X+1) <- (XFRAC, X, X + 1) + (XYBIS, XYBIS+1, 0) * 24bits add
If OldX (i.e. X before 24bits add)<> X OR X+1 changed (i.e. incremented) THEN CALL DotXYC

(XYBIS, XYBIS+1) <- (X,X+1)
(YFRAC, Y,Y+1) <- (YFRAC, Y, Y + 1) - (XYBIS, XYBIS+1,0) * 24bits substract
If OldY (i.e. Y before 24bits substract)<> Y OR Y+1 changed (i.e. decremented) THEN CALL DotXYC:IF (Y,Y+1)= (0,0) THEN END
JMP ]LOOP

DotXYC.. sets the four pixels in a symmetry around the [CX, CY] circle center: [CX-X,CY-Y], [CX-X,CY+Y], [CX+X,CY+Y] and [CX+X,CY-Y]..