function [y] = sample(a,b,n)
y = (a : (b-a)/n : b);