function [y] = S(x)
y = cumsum([0.0, x]);