%% Question 10

F = input('Enter a value for F:\n');
a = [0.246 0.250 0.643]';
b = [0.806 1.000 1.213]';

x = linspace(0, 25);
X = repmat(x, [3 1]);

A = repmat(a, [1 length(x)]);
B = repmat(b, [1 length(x)]);

Y = F * exp(-A .* X) .* sin(B .* X);
plot(X', Y')
xlabel('position');
ylabel('wave amplitude');
legend('ethyl alcohol', 'water', 'oil');
title('wave amplitude versus position');