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