t = linspace(0, 10, 1000);
A0 = 1;%начальная амплитуда(cм)
k = 9.9;%коэффициент упругости пружины(cН/cм)
ro = 11300/(10^-6);%плотностьсвинца(кг/cм^3)
m = 1;%масса свинца(кг)
nu = 35;%вязкость глицерина(cН*с/cм^2)
V = (m/ro);%объём свинцового шара(cм^3)
r = (((3*V)/(4*pi))^(1/3));%радиус свинцового шара(м)
b = 6*pi*r*nu;%коэффициент сопротивления глиценрина
alpha = b/(2*m);%коэффициент затухания
w = sqrt(k/m-alpha^2);
A = A0*exp(-alpha*t);
x = A.*cos(w*t);
figure("Name", "График затухающих колебаний", "NumberTitle", "off");
hold on;
xlabel("t, с");
ylabel("x(t),A(t), см");
grid minor; 
set(gca, 'xtick', t(1):0.5:t(end), 'xlim', [t(1) t(end)]);
set(gca, 'ytick', linspace(-A0, A0, 21), 'ylim', [-A0 A0]);
title("График затухающих колебаний");
plot(t, x);
plot(t, A);
plot(t, -A);
legend("x(t)", "A(t)", "-A(t)", -1);


