gnuplot

Вот пример прямоугольного
set yrange [-1:3]
f(x) = x<=0 ? 0 : x>0 && x <=1 ? 1 : 0
plot f(x)

А чтобы периодический?

Вот пример периодического прямоугольного:
gnuplot> a = 5 //длина промежутка между импульсами
gnuplot> b = 3 //длина импульса
gnuplot> c = 1 //высота импульса
gnuplot> set xrange [-4*(a+b) : 4*(a+b)]
gnuplot> set yrange [-c : 3*c]
gnuplot> f(x) = x < ceil(x/(a+b))*(a+b) && x > floor(x/(a+b))*(a+b)+a ? c:0
gnuplot> plot f(x)