sim('WoodandBerry1')
figure(1)
plot(t,xd+0.96,'r')
figure(2)
plot(t,xb+0.02,'g')
sim('WoodandBerry2')
figure(3)
plot(t,xd+0.96,'y')
figure(4)
plot(t,xb+0.02,'m')
