Computational algorithm for solving the diophantine equations 2^n ± α · 2^m + α^2 = x^2 - Szalay László (HOUSTON JOURNAL OF MATHEMATICS, 2020)