Superconductivity in the two-dimensional Hubbard model: Gutzwiller wave function solution