Use max():
x=[-10:0.02:10]
y=exp(-(x-1).^2./2)+3.*exp(-(x-2).^2./2)
plot(x,y, 'b-', 'LineWidth', 2)
grid on;
% Find the max
[yMax, indexAtMax] = max(y)
% Put a circle around the max.
hold on;
plot(x(indexAtMax), yMax, 'ro', 'MarkerSize', 12, 'LineWidth', 2);
% Draw a line from the x-axis up to the max.
line([x(indexAtMax), x(indexAtMax)], [0, yMax], 'Color', 'r', 'LineWidth', 2);