加快运动学问题的依存循环-Matlab

Question

我将使用这个简单的，独立的函数，应该将其作为输入：初始位置，力，时间步长和迭代次数；并返回输出：时间和位置向量。

function [time,pos] = example(pos0, force, N, dt) 

clc;

pos = zeros(N+1,3);
vel = zeros(N,3);
pos(1,:)=pos0;


for i = 1:N

        time(i+1) = i*dt;

        % solver finds the velocities given forces and previous positions
        vel(i,:) = force./(norm(pos(i,:)+1));

        % evolution calculates new positions using new velocities
        pos(i+1,:) = pos(i,:) + vel(i,:)*dt;

end

end

鉴于并行计算不是一种选择（循环项取决于先前的迭代），如何提高for循环的速度？ 矢量化将是一种选择吗？我该怎么办？