In this paper we study the quantum entanglement of two identical qubits which interact with external electromagnetic elds and the thermostat and each other by time-dependent dipole-dipole interaction. We develop original method in path integral approach for numerical calculation of density matrix and Peres-Horodecki parameter (the measure of qubits entanglement). The results indicate on possibility of high-entanglement states exiting and long-time non-destructive control of them