The paper presents an application of the homotopy analysis method for solving the one-phase fractional inverse Stefan design problem. The problem was to determine the temperature distribution in the domain and functions describing the temperature and the heat flux on one of the considered area boundaries. It was demonstrated that if the series constructed for the method is convergent then its sum is a solution of the considered equation. The sufficient condition of this convergence was also presented as well as the error of the approximate solution estimation. The paper also includes the example presenting the application of the described method. The obtained results show the usefulness of the proposed method. The method is stable for the input data disturbances and converges quickly. The big advantage of this method is the fact that it does not require discretization of the area and the solution is a continuous function.