Вот такой вопрос: зачем искать изолинии численными методами, если ваша функция задана аналитически?
Ну, а насчет рисования графиков — есть же уйма библиотек (MathGL, gnuplot и т.п.), упрощающих этот процесс…
Как понимаю, матрицу точек можно получать разными способами.
Тем же Кригингом, к примеру.
Вполне добротная подборка. Как альтернативу, можете в комменты ссылок подбросить, если посчитаете полезным.
Рисование сеточных графиков трехмерных функций и изолиний к ним