polynomial Division

; polynomialDivision
; Array[0] = Solution; Array[1] = Remainder
Func polynomialDivision($p, $q)
    Local $s[UBound($p,1)-UBound($q,1)+1], $i, $j
    For $i = UBound($p,1)-1 To UBound($q,1)-1 Step -1
        $s[$i-UBound($q,1)+1] = $p[$i] / $q[UBound($q,1)-1]
        For $j = UBound($q,1)-1 To 0 Step -1
            $p[$i-UBound($q,1)+1+$j] -= $s[$i-UBound($q,1)+1] * $q[$j]
        Next
    Next
    Local $a[2] = [$s, $p]
    Return $a
EndFunc