En Matlab/Octave el algoritmo queda como sigue.