We present a new method of calculating scalar propagator and vertex fu
nctions in the two-loop approximation, for arbitrary masses of particl
es, It is based on a double integral representation, suitable for nume
rical evaluation, Real and imaginary parts of the diagrams are calcula
ted separately, so that there is no need to use complex arithmetics in
the numerical program.