function dotq = one_dim_link (t, q)
    global J P
    theta = q(1);
    omega = q(2);
    dottheta = omega;
    dotomega = (1/J)*(-P*cos(theta) + torque(t));
    dotq = [ dottheta; dotomega ];
end
