Przedstawiono algorytm obliczeń polegający na połączeniu metody siatek z interpolacją wyników przy przechodzeniu do siatek o większej gęstości. Postępowanie takie zapewnia dużą dokłądność wyników przy małej liczbie węzłów sieci i może pełnić rolę sprawdzaniu w stosunku do metod wariacyjnych.