We present a fast one dimensional hybrid method to efficiently simulate
extensive air showers up to the highest observed energies. Based on
precalculated pion showers and a bootstrap technique, our method predicts the
average shower profile, the number of muons at detector level above several
energy thresholds as well as the fluctuations of the electromagnetic and
hadronic components of the shower. We study the main characteristics of
proton-induced air showers up to ultra-high energy, comparing the predictions
of three different hadronic interaction models: SIBYLL 1.7, SIBYLL 2.1 and
QGSjet98. The influence of the hadronic interaction models on the shower
evolution, in particular the elongation rate, is discussed and the
applicability of analytical approximations is investigated.Comment: 17 pages (Latex), 15 postscript figures, 7 table