Polynomial synthetic division: Difference between revisions
Content added Content deleted
Walterpachl (talk | contribs) m (add REXX) |
Walterpachl (talk | contribs) (→{{header|REXX}}: extended to support order of divisor >1) |
||
Line 142: | Line 142: | ||
=={{header|REXX}}== |
=={{header|REXX}}== |
||
<lang rexx>/* REXX |
<lang rexx>/* REXX Polynomial Division */ |
||
/* 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_dd '1 -12 0 -42' |
||
Call set_dr '1 -3' |
Call set_dr '1 -3' |
||
q.0=0 |
q.0=0 |
||
⚫ | |||
do While dd.0> |
do While dd.0>=dr.0 |
||
q=dd.1/dr.1 |
q=dd.1/dr.1 |
||
Do j=1 To |
Do j=1 To dr.0 |
||
dd.j=dd.j-q*dr.j |
dd.j=dd.j-q*dr.j |
||
End |
End |
||
Line 154: | Line 158: | ||
Call shift_dd |
Call shift_dd |
||
End |
End |
||
⚫ | |||
⚫ | |||
⚫ | |||
Exit |
Exit |
||
Line 199: | Line 202: | ||
list=list','q.i |
list=list','q.i |
||
End |
End |
||
Return list']' |
Return list']' |
||
mk_list_dd: |
|||
list='['dd.1'' |
|||
Do i=2 To dd.0 |
|||
list=list','dd.i |
|||
End |
|||
Return list']' |
|||
</lang> |
|||
{{out}} |
{{out}} |
||
<pre>[1,-12,0,-42] / [1,-3] |
<pre>[1,-12,0,-42] / [1,-3] |