We describe our method of determining solar cosmic-ray flux and spectrum in the past, based on the comparison of different cosmogenic isotopes. For the period, AD 1781–1950, we have detected several intervals with a high probability of powerful solar flares.