编制一个 fu 函数文件，函数关系为 a(m,f) = f/m,再编制一个 M 文件，求出初速度 为 0 的直线运动，当 t=2s，4s，9s 时的速度、加速度、路程。

好的，我理解了你的要求。下面是 fu 函数的代码：

function [a] = fu(m, f)
% 计算加速度
a = f / m;
end

接下来是 M 文件的代码，其中包括了三个时间点 t=2s，4s，9s 的速度、加速度、路程的计算：

% 定义变量
m = 10; % 单位：kg
f = 100; % 单位：N
t = [2, 4, 9]; % 单位：s

% 调用函数计算加速度
a = fu(m, f);

% 计算速度
v = a .* t;

% 计算加速度
a_vec = repmat(a, 1, length(t));

% 计算路程
s = 0.5 * a .* t.^2;

% 输出结果
disp('t(s)   v(m/s)   a(m/s^2)   s(m)')
disp([t; v; a_vec; s])

请注意，以上代码中的“路程”指的是直线运动的位移。如果要计算直线运动的路径长度，需要使用积分的方法进行计算。