for b=[ 20 50 100];
m=1000;
%b =50;
num=[1];
den=[m,b];
G=tf(num,den)
step(G); pause
hold on; end
legend( '20' ,'50', '100')

1111111111111111111111

for m=[ 300 1000 2000]
 %m=1000;
b=500;
%b =50;
num=[1];
den=[m,b];
G=tf(num,den)
step(G); pause
hold on; end
legend( '300' ,'1000', '2000')
2222222222222222222222222222222222222
for K=[.05 .01 .5 1 1.5]
b=0.1;
j= .01;
R=1;
L=0.5;
num=K;
den=[(J*L) ((J*R)+(L*b)) ((b*R)+K^2)];
G= tf(num,den);
step(G);pause
hold on
end
legend ('.05 ','.01', '.5', '1', '1.5');
33333333333333333333333333333333333333
for R= [10 30 100];
L= 50;
Cap= .05;
A = [0 1/Cap; -1/L -R/L];
B=[0 ;1/L];
C=[1 0];
D=[0];
step(A,B,C,D); pause
hold on
end