use constant_iterator for dx, use predefined multiply functor for y*dx