Usage example for Lagrangeeval related to Figure 3.1
0001 % Usage example for Lagrangeeval 0002 % related to Figure 3.1 0003 0004 x = (-1:0.01:1)'; 0005 [~,xn] = size(x); 0006 y = ones(xn)./(ones(xn)+25*x.^2); 0007 plot(x,y,'k') % plot Runge's example 0008 hold 0009 n =16; 0010 inter = 2/n; 0011 nodes = (-1:inter:1)'; 0012 values = ones(n+1,1)./(ones(n+1,1)+25*nodes.^2); 0013 [y] = Lagrangeeval(nodes, values, x ); 0014 plot(x,y) % plot the interpolant 0015 plot(nodes, values,'o') % plot the data