| We introduce a new family of bivariate copulas which evolves according to the discrete process of heat equation. We prove the convergence of solution as well as the measure of dependence. Numerical experiments are also performed, which shows that our procedure works substantially well. |