Abstract. This paper deals with the numerical approximation of mild solutions of elliptic-parabolic equations, relying on the existence results of Bénilan and Wittbold (1996). We introduce a new and simple algorithm based on Halpern's iteration for nonexpansive operators (Bauschke, 1996;Halpern, 1967;Lions, 1977), which is shown to be convergent in the degenerate case, and compare it with existing schemes (Jäger and Kačur, 1995;Kačur, 1999).Mathematics Subject Classification. 65M12, 35K65, 35K55, 65N22.