%Rungen funktion interpolointi kuutiollisella splinillä
% Tiedosto rungesplini.m   HA 18.8.04
close all; clear
f=inline('1./(1+25*x.^2)','x')  
x=linspace(-1,1,100);                                
plot(x,f(x),'k');  hold on
for n=5:2:21                   % Piirretään n:n arvoilla 5,7,9,...,21 
    X=linspace(-1,1,n+1);                                          
    s=spline(X,f(X),x);
    plot(x,s)
end


