We calculate the pressure (p), the energy density (epsilon) and the baryon density (n(B)) of QCD at finite temperatures (T) and chemical potentials (mu). The recently proposed overlap improving multi-parameter reweighting technique is used to determine observables at nonvanishing chemical potentials. Our results are obtained by studying n(f) =2+1 dynamical staggered quarks with semi-realistic masses on N-t = 4 lattices