In this paper, we study the existence and the numerical estimates of solutions for a specific types of fractional differential equations. The nonlinear part of the problem, however, presupposes certain hypotheses. Particularly, for exact localization of the parameter, the existence of a non-zero solution is established, which requires the sublinearity of the nonlinear part at the origin and infinity. We also take into consideration several theoretical and numerical examples. One of the main novelties of this paper is to use the variational method to investigate the properties of solutions of boundary values problems involving Atangana–Baleanu fractional derivative for the first time.