This paper presents an efficient and comprehensive MATLAB code to solve two-dimensional structural topology optimization problems, including minimum mean compliance, compliant mechanism synthesis and multi-load compliance problems. The Unsmooth Variational Topology Optimization (UNVARTOP) method, developed by the authors in a previous work, is used in the topology optimization code, based on the finite element method (FEM), to compute the sensitivity and update the topology. The paper also includes instructions to improve the bisection algorithm, modify the computation of the Lagrangian multiplier by using an Augmented Lagrangian to impose the constraint, implement heat conduction problems and extend the code to three-dimensional topology optimization problems. The code, intended for students and newcomers in topology optimization, is included as an appendix (Appendix A) and it can be downloaded from https://github.com/DanielYago together with supplementary material.
