The idea is to approximate a surface part embedded into 3D non-plane contour by an arbitrary triangle grid.
Solving all three systems simultaneously one can obtain a new grid that will be the approximating minimal surface embedded into non-plane closed curve because of the minimum of the function \ \Pi where parameter \ j = 1,2,3.