//MAC5749	Análise de Formas e Classificação
//EP4		Thin-plate splines para campos vetoriais
//Aluno		Robson Ribeiro Correia

//Este arquivo contém definições de funções para o cálculo
//de thin-plate splines

//Função Ro
//Entrada:	(xi,yi), (xj,yj), posições dos pontos no plano cartesiano
//Saída:	distância euclidiana entre os pontos

function [D]=Ro(xi,yi,xj,yj)
   D=sqrt((xi-xj)^2+(yi-yj)^2)
endfunction


//Função G
//Entrada:	saída da função Ro 
//		(nessa implementação, a distância euclidiana entre dois pontos)
//Saída:	0		se Ro=0
//		Ro^2 log(Ro)	caso contrário

function [J]=G(r)
   if r == 0 then
      J=0;
   else
      J=r^2*log(r);
   end
endfunction