Fully comment poly.cl and add new tests.cl with some testing

poly.cl

                        (term-addition (car poly) y))
              (cdr poly)) (car poly))))))))
 
+;; Map a single polynomial onto each term of another polynomial
 (defun map-onto-poly (single-poly poly)
   (map 'list #'(lambda (x) (multiply-singles single-poly x)) poly))
 
+;; Multiply two singular polynomials together
 (defun multiply-singles (single-poly1 single-poly2)
   (cond
     ((equal (variable-symbol single-poly1) (variable-symbol single-poly2))
                         (exponent single-poly2)))
             (* (multiplier single-poly1) (multiplier single-poly2))))))
 
+;; Multiply two polynomials together, see the map-poly will be recursed
+;; and the poly will remain constant
 (defun multiply-terms (map-poly poly)
   (cond
     ((equal map-poly nil) nil)
     (T (remove nil (append (map-onto-poly (car map-poly) poly) 
                (multiply-terms (cdr map-poly) poly))))))
-
-; (poly+ '(((x 2) 3)) '(((y 2) 4)))
-; (poly+ '(((x 2) 3) ((y 2) 3)) '(((y 2) 4)))

tests.cl

+(require "poly.cl")
+
+(defun run-tests ()
+  (if (equal 
+        (poly+ '(((x 2) 5)) '(((x 2) 3)))
+        '(((x 2) 8)))
+    (format t "5x^2 + 3x^2 ... passed - simple addition~%")
+    (format t "5x^2 + 3x^2 ... failed - simple addition~%"))
+  (if (equal
+        (poly+ '(((x 2) 5) ((y 2) 6)) '(((x 2) 3) ((y 2) 3) ((z 3) 4)))
+        '(((z 3) 4) ((x 2) 8) ((y 2) 9)))
+    (format t "5x^2+6y^2 + 3x^2+3y^2+4z^3 ... passed - complex addition~%")
+    (format t "5x^2+6y^2 + 3x^2+3y^2+4z^3 ... failed - complex addition~%"))
+  (if (equal 
+        (poly- '(((x 2) 5)) '(((x 2) 3)))
+        '(((x 2) 2)))
+    (format t "5x^2 - 3x^2 ... passed - simple subtraction~%")
+    (format t "5x^2 - 3x^2 ... failed - simple subtraction~%"))
+  (if (equal
+        (poly- '(((x 2) 5) ((y 2) 6)) '(((x 2) 3) ((y 2) 3) ((z 3) 4)))
+        '(((z 3) -4) ((x 2) 2) ((y 2) 3)))
+    (format t "5x^2+6y^2 - 3x^2+3y^2+4z^3 ... passed - complex subtraction~%")
+    (format t "5x^2+6y^2 - 3x^2+3y^2+4z^3 ... failed - complex subtraction~%"))
+  )