A numerical assessment of the free energy function for fractional-order relaxation