function [j]=jerk(t)
    j = 4*%pi*%pi*cos(2*%pi*t);
endfunction

function [a]=acceleration(t)
    a = 2*%pi*sin(2*%pi*t);
    endfunction

function [v] = vitesse(t)
    v = 1 - cos(2*%pi*t);
endfunction

function [x] = position(t)
    x = (t - 0.5*(1/%pi)*sin(2*%pi*t))  ;
endfunction

t = 0:0.005:1;

j = jerk(t);
a = acceleration(t);
v = vitesse(t);
x = position(t);

scf(0);
clf;

subplot(4, 1, 1)
plot(t, x)
xtitle("position", "t/T", "$x/x_0$")

subplot(4, 1, 2)
plot(t, v)
xtitle("vitesse", "t/T", "$v/(x_0/\mathrm{T})$")

subplot(4, 1, 3)
plot(t, a)
xtitle("accélération", "t/T", "$a/(x_0/\mathrm{T}^2)$")

subplot(4, 1, 4)
xtitle("à-coup", "t/T", "$j /(x_0/\mathrm{T}^3)$")
plot(t, j)