% building example
% 2 body mass system

global k m A w

w = 55;     % rad/sec
k = 10^8;   % N/m
m = 25000;  % kg
A = 0.01;   % m
L = 2;      % m

IC = [ 0 0 0 0 ];

[t,y]=ode45('building', [0 10], IC');   
subplot(2,1,1)
plot(t,y(:,1),'r')
subplot(2,1,2)
plot(t,y(:,2),'r')


xlabel('t, sec')
ylabel('y')

