Modelagem computacional de curvas

Estou trabalhando na produção da modelagem computacional de curvas, interessame neste caso
a interseção de variedades do tipo f(x,y,z) = 0.

Veja o artigo em andamento sobre este assunto na página de pré-prints do Curso de
Matemática em Sobral (procure o
link pré-prints).
É um trabalho em andamento que está descrito e faz referência a um programa
com que estou fazendo algumas experiências gráficas.

A idéia consiste em partir de uma condição inicial, um ponto P=(a,b,c) que pertença a interseçao e produzir uma sucessão de vetores tangentes a esta curva de módulo suficientemente pequeno para conseguir uma visão da curva.

Este é um item que deve me conduzir a modelar soluçoes de equaçoes diferencias do tipo

grad(u) = f

em que f é um par de equaçoes.modelagem da interseção de superficies

O programa mencionado acima produz este gráfico
onde o erro fica mais explicito.