function b = basis(torque, rate)

  power = torque.*rate;

  b = [
    ones(size(torque)), ...
    torque, ...
    power, ...
    torque.^2, ...
    power.^2];

end
