yz_direction = np.matrix(np.zeros([G(0.001).shape[0],
                                               len(Zeros_G)]))