u=RF(ps2(u,k),qs2(u,k),rs2(u,k));