mod_G_Gd = np.zeros([np.shape(G(0.0001))[0], len(w)])