The problem of the synthesis of a homogeneous Lyapunov function for an asymptotically stable homogeneous system is studied. First, for systems with nonnegative degree of homogeneity, several expressions of homogeneous Lyapunov functions are derived, which depend explicitly on the supremum or the integral (over finite or infinite intervals of time) of the system solutions. Second, a numeric procedure is proposed, which ensures the construction of a homogeneous Lyapunov function. The analytical results are illustrated by simulations.