Our model generalizes the differential D(E) and integral D(>E) spectra of cosmic rays (CR) during the 11-year solar cycle. The empirical model takes into account galactic (GCR) and anomalous cosmic rays (ACR) heliospheric modulation by four coefficients. The calculated integral spectra in the outer planets are on the basis of mean gradients: for GCR – 3%/AU and 7%/AU for anomalous protons. The obtained integral proton spectra are compared with experimental data, the CRÈME96 model for the Earth and theoretical results of 2D stochastic model. The proposed analytical model gives practical possibility for investigation of experimental data from measurements of galactic cosmic rays and their anomalous component.