interpolation after delaunay triangulation

is there a function that performs interpolation of data in a grid created by delaunay triangulation? for example: X = [1,3,6,2,2]; Y = [2,4,1,6,8]; Z = [12,15,10,8,20];

after performing DT = delaunaytri(X,Y,Z);

i would like to calculate the Z values at the next 3 points: x =[4,1,5]; y = [1,3,2];

thank you,