# HG changeset patch
# User huffman
# Date 1234995424 28800
# Node ID 17ddfd0c350697fed49c0dc571fbc635849d4fed
# Parent  666f5f72dbb544f2f29cd1c960212922b6dde268
composition of polynomials

diff -r 666f5f72dbb5 -r 17ddfd0c3506 src/HOL/Polynomial.thy
--- a/src/HOL/Polynomial.thy	Wed Feb 18 12:24:06 2009 -0800
+++ b/src/HOL/Polynomial.thy	Wed Feb 18 14:17:04 2009 -0800
@@ -1120,6 +1120,11 @@
 
 subsection {* Synthetic division *}
 
+text {*
+  Synthetic division is simply division by the
+  linear polynomial @{term "x - c"}.
+*}
+
 definition
   synthetic_divmod :: "'a::comm_semiring_0 poly \<Rightarrow> 'a \<Rightarrow> 'a poly \<times> 'a"
 where [code del]:
@@ -1246,6 +1251,32 @@
   by (simp add: expand_fun_eq)
 
 
+subsection {* Composition of polynomials *}
+
+definition
+  pcompose :: "'a::comm_semiring_0 poly \<Rightarrow> 'a poly \<Rightarrow> 'a poly"
+where
+  "pcompose p q = poly_rec 0 (\<lambda>a _ c. [:a:] + q * c) p"
+
+lemma pcompose_0 [simp]: "pcompose 0 q = 0"
+  unfolding pcompose_def by (simp add: poly_rec_0)
+
+lemma pcompose_pCons:
+  "pcompose (pCons a p) q = [:a:] + q * pcompose p q"
+  unfolding pcompose_def by (simp add: poly_rec_pCons)
+
+lemma poly_pcompose: "poly (pcompose p q) x = poly p (poly q x)"
+  by (induct p) (simp_all add: pcompose_pCons)
+
+lemma degree_pcompose_le:
+  "degree (pcompose p q) \<le> degree p * degree q"
+apply (induct p, simp)
+apply (simp add: pcompose_pCons, clarify)
+apply (rule degree_add_le, simp)
+apply (rule order_trans [OF degree_mult_le], simp)
+done
+
+
 subsection {* Order of polynomial roots *}
 
 definition