yp_direction = np.matrix(np.zeros([G(0.001).shape[0],
                                               len(Poles_G)]))