New algorithms, combining asymptotic numerical method (ANM) and method of fundamental solutions, are proposed to compute bifurcation points on branch solutions of a nonlinear bi-harmonic problem. Three methods, mainly based on asymptotic developments framework, are then proposed. The first one consists in exploiting the ANM step accumulation close to the bifurcation points on a solution branch, the second method allows the introduction of an indicator that vanishes at the bifurcation points, and finally the first real root of the Padé approximant denominator represents the third bifurcation indicator. Two numerical examples are considered to analyze the robustness of these algorithms. KEYWORDS asymptotic numerical method, bifurcation branch, bifurcation indicator, bi-harmonic problem, method of fundamental solutions 1 Numer Methods Partial Differential Eq. 2019;35:2091-2102. wileyonlinelibrary.com/journal/num