I want to prove how Sarah works and how I can create new models. I was thinking about a model that I saw in an article in order to get the results that it says. If I got the results that the article says, I could feel more confortable that I know how to use SARAH for my own problems.
In that article, they work with a model for an SU(5) gauge group with 5,5 ̅,10,10 ̅ (the 5 ̅ and 10 ̅ are the bars of 5 and 10 but I cannot wirte it well) and 24 chiral multiplets with multiplicities of (4, 7, 3, 0, 1) respectively, taking its contents to be one 24 of Higgs for the SU(5) breaking, 4(5+5 ̅) (four of the 5 and the 5 ̅) of Higgs some of which will be used for electroweak breaking and the remaining 3(5 ̅+10) which I identify as matter. Consistent with SU(5) invariance and R-parity conservation, they write a superpotential of the form

Which I don't understand how to write
Basically, I want to use the RGEs for the Yukawa couplings h_{t} from M_{G} (Grand unification scale) to M_{W} (electroweak scale)
