adaptive

'adaptive.m' allows to efficiently sample a function using a reduced number of points. It works by iteratively adding new points where needed. It is especially useful for functions which are computationally intensive (e.g. involve solving a differential equation).
Examples:
% using default options
[xi, yi] = adaptive(@(x) sin(x)*sin(2*x), [-20,20]);
% explicitly providing initial array an refinement option
f = @(x) exp(-x.^2/4).*sin(3*x);
[xi,yi] = adaptive(f, [ -5:5], 'minAngle',0.8*pi);