Aim of this study was to present a P-wave model, based on a linear combination of Gaussian functions, to quantify morphological aspects of Pwave in patients prone to atrial fibrillation. Five minutes ECG recordings were performed in 25 patients with permanent dual chamber pacemakers set at 40/min in order to have spontaneous beats. ECG signals were acquired using a 32-lead mapping system for high-resolution biopotential measurement (ActiveTwo, Biosemi, The Netherlands, sample frequency 2 kHz, 24 bit resolution). Four healthy subjects were also recorded as a control group. Up to 8 Gaussian models have been computed for each averaged P-wave extracted from every lead. The P-wave morphology is then evaluated by the following parameters: best model orders @ degrees of freedom adjusted R-square (AdjRsq) =97.5%; minimum (sigmamin) and maximum (sigmamax) standard deviation of the Gaussians included in the model, number of relative maxima and minima (max+min), and zeroes of the fit. Significant differences in the best model order were obtained between the control group and patients group. Accordingly, the number of relative maxima and minima was higher in the patient group. These parameters might all be markers of the fractionated electrical activity that characterizes paroxysmal AF patients in sinus rhythm.