n = 4;
x_data = linspace ( 0.0, 1.0, n );
y_data = exp ( x_data );

x = linspace ( -0.5, 1.5, 25 );
y = exp ( x );
y2 = zeros ( size ( x ) );
for i = 1 : n
  y2 = y2 + y_data(i) * ell(i,x_data,x);
end

plot ( x, y, 'b-', ...
       x, y2, 'r-', ...
       x_data, y_data, 'k.', 'markersize', 25 );

print ( '-dpng', 'exp_any.png' );

