Matlab or Octave Script for this problem.

% This octave/Matlab script demonstrates the state space model of the 

J=0.01;
b=0.1;

K=0.01; R=1; L=0.5; A1=[[0,1,0];[0,-b/J,K/J];[0,-K/L,-R/L]]; B1=[0,0,1/L]'; C1=[0,1,0]; % The output is the angular velocity omega. sys1=ss(A1,B1,C1);  % Note that theta cannot be controlled. % So we use only the theta dot, and i_a as state variables. A2=[[-b/J, K/J];[-K/L,-R/L]]; B2=[0,1/L]'; C2=[1,0]; sys2=ss(A2,B2,C2); figure(1) clf;