Можно ли одной операцией в octave, matlab, python найти вектор с центральными разностями dx. Например, на octave элементарные действия:

t=[2*pi:-2*pi/n:pi/n];
x=cos(t);
y=sin(t);

dx(1)=(x(2)-x(n))/2.;
for i=2:n-1
dx(i)=x(i+1)-x(i-1)/2.;
end
dx(n)=(x(1)-x(n-1))/2.;

интересует сделать короткий просто-читаемый код с матричными операциями, а не поэлементными.

P.S. в C++ с использованием % на две строчки код короче, а в octave операции целочисленного деления не нашел