Feynman integrals of the banana topology are essential to understanding Feynman integrals' functional space and scattering amplitudes beyond multiple polylogarithms. I will focus on equal-mass cases at three and four loops in this talk. Recently we managed to obtain epsilon factorized differential equations via a similar basis as the two-loop sunrise. Based on the two-, three- and four-loop data, some interesting patterns emerge and may extend to all loops.