<DIV>A method for enhancing the stability and robustness of explicit schemes in computational fluid dynamics is presented. The method is based in reformulating explicit schemes in matrix form, which can then be modified gradually into semi or strongly implicit schemes. From the point of view of matrix-algebra, explicit numerical methods are special cases in which the global matrix of</DIV>
<DIV>coefficients is reduced to the identity matrix I. This extreme simplification leads to severe limitation of their stability range, hence of their robustness. In this paper it is shown that a condition, which is similar to the Courant–Friedrich–Levy (CFL) condition, can be obtained from the stability requirement of inversion of the coefficient matrix. This condition is shown to be relax-able, and that a class of methods that range from exp