An algorithm for synthesis of a class of multivariable positive real functions