This study investigates into the analysis of ferromagnetic thin films with a body centered cubic lattice and three spin layers, utilizing the solution of the fourth-order perturbed Heisenberg Hamiltonian equation with seven magnetic energy parameters. Spin-exchange interaction, magnetic dipole interaction, second-order magnetic anisotropy, fourth-order magnetic anisotropy, applied magnetic field, demagnetization energy, and stress-induced anisotropy were all taken into account. According to 3D plots, the minimum order of energy was observed when the second order magnetic anisotropy constant in the middle spin layer is less than those of the bottom and top spin layers. In all cases, the values of stress-induced anisotropy at the maxima of 3D plots are exactly the same when the values of the second-order magnetic anisotropy constants of the bottom, middle, and top spin layers are interchanged with each other. In 2D plots, the angle between consecutive magnetic easy and hard directions is approximately 90 degrees in all cases. Additionally, the magnetic easy and hard directions were observed to have exactly the same values when the second-order magnetic anisotropy constant of the spin layers’ changes. These results were compared with the results obtained using the second and third order perturbed Heisenberg Hamiltonian.