matlab code for simulation in laser ablation

ed on the heat conduction equation: \[ \rho C_p \frac{\partial T}{\partial t} = k \nabla^2 T + Q \] Where: \(\rho\): density \(C_p\): specific heat capacity \(k\): thermal conductivity \(Q\): heat source term from lase