double precision function areteideale( xyz, direction ) double precision xyz(3), direction(3) areteideale = 10 return end