La simulación de los últimos estadios de la vida de las estrellas necesita de métodos hidrodinámicos multidimensionales debido a la pérdida de simetría esférica de los sistemas y a la frecuente aparición de inestabilidades hidrodinámicas. Hasta fechas recientes casi todos los métodos hidrodinámicos utilizados, tanto de tipo Lagrangiano como Euleriano eran de naturaleza explícita. En este trabajo nos hemos propuesto la construcción de un método Lagrangiano implícito, que permitirá abordar sistemas en quasi equilibrio. A continuación se describen los fundamentos físicos y matemáticos del código construido, las dificultades encontradas y la propuesta de trabajo futuro.