function [pp]=qspliner(x,y,ypn)
n = length(x);
h = zeros(n,1);
d = zeros(n,1);
b= zeros(n,1);
s2= ones(n,1);
t =1;
for i=1:n-1
d(i) = 2.*((y(i+1)-y(i))./(x(i+1)-x(i)));%rechte seite der gleichung
h(i+1) = x(i+1) - x(i);%gitterweiten
end
d(n)= ypn; 
M = diag(ones(n-1,1),1) + eye(n,n);
s1 = M\d; %erste ableitungen an den knoten
while t<n
s2(t+1) = ((s1(t+1) - s1(t))./(x(t+1) - x(t)));%2. ableiteungen an den knoten
t=t+1;
end
a=s2(2:n)./2;% koeffizient an x^2(anpassung an mkpp)
b=s1-h.*s2;%koeffizient an x
c=y'-h.*s1+1/2.*(h.^2).*s2;;%koeffizient an x^0
pp = mkpp(x,[a b(2:n) c(2:n)]);
end