In this paper, first, we will try to introduce the gravitational domain wall as a physical system. In the second step, we also introduce the Hun differential equation as a mathematical tools. We factorize the known Heun’s equation as form of operators P +, P − and P0. Then we compare the differential equation of gravitational domain wall with corresponding Hun equation. In that case the above-mentioned operators can be obtained for the gravitational system by the comparing process. Finally, we employ such operators and achieve the corresponding symmetry algebra with the usual commutation relation of operators to each other. Here, by having such operators, we investigate the stability of system