  We show that there is an inconsistency in the class of solutions obtained in
Phys. Rev. D {\bf 95}, 084037 (2017). This inconsistency is due to the
approximate relation between lagrangian density and its derivative for
Non-Linear Electrodynamics. We present an algorithm to obtain new classes of
solutions.
