We put forward a new method to calculate the leading order energy correlators in N=4 Yang-Mills theory. They are infrared-safe observables, and can be constructed from finite integrals. We aim to establish an integration by parts system and differential equations that exculde as many divergent integrals as possible. To avoid the divergent integrals, from algebra geometry point of view, we apply syzygy method to generate IBP system, and lift method to differential equations. For these kinds of atypical Feynman integrals, it's canonical differential equation matrix is block upper triangular. Our new method can work smoothly for three-point energy correlator. There are elliptic integrals for four-point energy correlator in arbitrary angle shape, our method can work well for unitary transcendental sector. We are going to check how many integrals are elliptic and if there is more complex curves will appear in four point energy correlator.