The recently developed formalism for the evaluation of nuclear form factors in neutrinoless double beta decay is applied to Ca-48, Ge-76, Se-82,Mo-100,Te-128 and Te-130 nuclei. Explicit analytical expressions that follows from this theoretical development, in the single-mode model for the decay of 48Ca, have been worked out. They are useful both for testing the full numerical calculations, and for analytically checking the consistency with other formalisms. Large configuration space calculations are compared with previous studies, where alternative formulations were used. Yet, besides using the G-matrix as residual interaction, we here use a simple delta-force. Attention is paid to the connected effects of the short range nuclear correlations and the finite nucleon size. Constraints on lepton number violating terms in the weak Hamiltonian (effective neutrino Majorana mass and effective right-handed current coupling strengths) are deduced. (C) 1999 Elsevier Science B,V.