function [ x, gamma, omega0, sq ] = dho_position( m, c, k, x0, t )
%UNTITLED3 Summary of this function goes here
%   Detailed explanation goes here

gamma=c/(2*m); omega0=sqrt(k/m);
zeta=c/(2*sqrt(m*k)); phi=acos(sqrt(1-zeta^2));
sq=sqrt(1-zeta^2);

x= x0/sq * exp(-gamma*t) .* cos(omega0*t*sq - phi);

end

