d=0:2.5:10;
t=(0:30:60)';
c=[95,14,0,0,0;88,48,32,12,6;67,64,54,48,41];
di=0:2:10;
ti=(0:20:60)';
ci=interp2(d,t,c,di,ti)
surf(di,ti,ci);