function dx = Rrobot(t,x); % R robot arm dx = zeros(2,1); dx(1) = x(2); dx(2) = -135*x(2)-90*x(1)+30*pi;