x=[50; 150; 250; 350];
y=[-1845; -456; 1277; 2047];

% not-a-knot Bedingung wird bei spline() genutzt
s = spline(x,y,175);
pp = spline(x,y);

fprintf('c=%.6e\n',pp.coefs);
fprintf('s=%.6e\n',s);

% interp1 cubic
%s = interp1(x,y,'pchip');
%fprintf('pchip\n');
%fprintf('s=%.6e\n',s);