File:Loi de mouvement polynome 3 4 5.svg

function [j]=jerk(t)
    j = 60*(6*(t - 1).*t + 1);
endfunction

function [a]=acceleration(t)
    a = 60*((2*t - 3).*t + 1).*t;
    endfunction

function [v] = vitesse(t)
    v = 30*((t - 2).*t + 1).*t.*t;
endfunction

function [x] = position(t)
    x = ((6*t - 15).*t + 10).*t.*t.*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)