function dr=velocidades(JRR,edQ,MDH)

for i=1:1:size(MDH,3)
    dR(:,:,i)=JRR(:,:,i)*edQ;
end

dr=dR;