let mk_div a b = 
  mk_mult a (mk_expt (-1) b)