function esthecis2(udat,Ddat,clr);

u = -udat(end):0.01:udat(end);
ud = [(-1).*fliplr(udat(2:end)) udat];
Dd = [fliplr(Ddat(2:end)) Ddat];
D = interp1(ud,Dd,u,'cubic');

hold on;
plot(u,D,clr);
grid on;