Un Método de Penalización Interior de C0 para la Ecuación Cristalina del Campo de Fase