fixed bug in the computation of the basis functions
computation of Sigma appears to be wrong (in comparison with matlab implementation)
improved output. should think about evaluating the relative distance!