Mercurial > hg > octave-lyh

comparison scripts/polynomial/mkpp.m @ 6746:a8105a726e68

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
[project @ 2007-06-19 08:18:34 by dbateman]
author dbateman
date Tue, 19 Jun 2007 08:18:34 +0000
parents 01036667884a
children 93c65f2a5668
comparison
equal deleted inserted replaced
6745:ca6668e47518 6746:a8105a726e68
62 %! t=[sin(x),cos(x)]; 62 %! t=[sin(x),cos(x)];
63 %! m=diff(t)./(x(2)-x(1)); 63 %! m=diff(t)./(x(2)-x(1));
64 %! b=t(1:4,:); 64 %! b=t(1:4,:);
65 %! pp = mkpp(x, [m(:),b(:)]); 65 %! pp = mkpp(x, [m(:),b(:)]);
66 %! xi=linspace(0,pi,50); 66 %! xi=linspace(0,pi,50);
67 %! plot(x,t,"x;control;",xi,ppval(pp,xi),";interp;"); 67 %! plot(x,t,"x",xi,ppval(pp,xi));
68 %! legend("control","interp");