Anonymous user
Polynomial long division: Difference between revisions
m
→{{header|Racket}}
Line 1,406:
(if (zero? pi) d i)))
(define (lead p) (vector-ref p (deg p)))
(define (mono
(define (poly*cx^n c n p) (vector-append (make-vector n 0) (for/vector ([pi p]) (* c pi))))
(define (poly+ p q) (poly/lin 1 p 1 q))
Line 1,427:
[(< N D) (values 0 n)]
[else (define c (/ (lead n) (lead d)))
(define
(define
(define-values (q1 r1) (poly/ r d))
(values (poly+ q q1) r1)]))
| |||