clear clc

X=[25,36,49,64,81]; Y=[5,6,7,8,9];
x=[25:0.4:81];
y=cspline(X,Y,x); 

plot(X,Y,'r*',x,y)
figure

plot(x,y-sqrt(x))