Includes all functionality sans sampling via Gaussian distributions

p_zt=h*(z^k)/g?