% main.m — click "Run" to execute this in your browser.

x = linspace(0, 4*pi, 400);
y = waves(x);

fprintf('peak value:  %.4f\n', max(y));
fprintf('zero crossings: %d\n', sum(abs(diff(sign(y))) > 0));

plot(x, y);
title('damped sine: e^{-x/10} sin(x)');
xlabel('x');
ylabel('y');
grid on;