ptozp_subst

[f  s  [p q] v] ptozp_subst g   (ox_asir function)
It returns ptozp(subst(f,s,p/q))
Example 1:  [ (x Dx - s) (s) [2 3] (x,s)] ptozp_subst 
   ===>     3 x Dx - 2