  In this work, we explore the dynamics of entanglement of an isolated quantum
system consisting of two time-dependent, coupled harmonic oscillators. Through
the use of a numerical method that relies on the estimation of the system's
Wigner representation by a specific Gaussian function, we investigate the time
evolution of the entanglement entropy after an instant quench in the inherent
parameters of the system. Besides, from the comparison of the results obtained
from the analytical expression for the time-dependent von Neumann entropy with
the numerically computed entropy data, the effectiveness of the numerical
method is tested for a variety of angular frequency combinations. Also, we
analyze how the entropy of entanglement change as a function of time.
