The study of liquid crystals is one of the active areas of research in physics. In this paper, the MLPG and direct MLPG (DMLPG) methods are used for the numerical study of the coupled nonlinear sine-Gordon equations in two dimensions arising from the modeling of some phenomena in the liquid crystals and superconductors. A comparison of numerical results between the MLPG and DMLPG methods highlight the precision and speed of the DMLPG method.