function xdot=mech_ode(t,x)
%-----------------------------------------------------------------
% Copyright 1998, by Bogumila and Zbigniew Mrozek, Krakow. For use
% with the book "MATLAB 5.x, SIMULINK 2.x, poradnik uzytkownika"
% (PLJ Warszawa, 1998, ISBN 83-7101-376-0).
%-----------------------------------------------------------------
% see page 196, chapter 8.5.5 of the book
%
% mech_ode.m zawiera model ukladu drgan II-go rzedu
% dx2/dt2 + beta*x + x = 0
% mech_ode jest wolane przez funkcje ode23 z M-pliku mechplot.m
% **** supporting M-file for mech_gui.m ******
global beta t2 % NIE oddzielac przecinkiem !
xdot=zeros(2,1); %zerowy wektor dwuelementowy
xdot(1)=x(2); % x(1)= x x(2)= dx//dt
xdot(2)= -1*x(1) -beta*x(2);