Hello John

This algorithm assumes smooth functions and the manipulation of the error terms has calculus background to them.

Another version of this algorithm could be obtained for continuous functions which are not smooth, and no calculus smooth operators needed.

In this case the error term will be minimized by a global minimizer e.g. Differential Evolution and that could run might fast on the new GPU servers.