  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.
