Polynomial synthetic division: Difference between revisions
→{{header|REXX}}: extended to support order of divisor >1
Walterpachl (talk | contribs) m (add REXX) |
Walterpachl (talk | contribs) (→{{header|REXX}}: extended to support order of divisor >1) |
||
Line 142:
=={{header|REXX}}==
<lang rexx>/* REXX
/* extended to support order of divisor >1 */
call set_dd '1 0 0 0 -1'
Call set_dr '1 1 1 1'
Call set_dd '1 -12 0 -42'
Call set_dr '1 -3'
q.0=0
Say list_dd '/' list_dr▼
do While dd.0>
q=dd.1/dr.1
Do j=1 To
dd.j=dd.j-q*dr.j
End
Line 154 ⟶ 158:
Call shift_dd
End
▲Say list_dd '/' list_dr
▲say 'Quotient:' mk_list_q() 'Remainder:' dd.1
Exit
Line 199 ⟶ 202:
list=list','q.i
End
Return list']'
mk_list_dd:
list='['dd.1''
Do i=2 To dd.0
list=list','dd.i
End
Return list']'
</lang>
{{out}}
<pre>[1,-12,0,-42] / [1,-3]
|