In the talk, I will present an efficient algorithm to compute the decomposition of the Jacobians of modular curves of arbitrary level, using modular symbols. This is obtained by working intrinsically with the curve, unlike previous methods. I will also discuss the possible consequences for deriving equations of modular curves.