import numpy as np in MC2004 recipe