clc
clear
% Mayer formula
Basin_name = {'A';'B';'C';'D';'E'};
A = [12,15,28,91,4]; % area in Square kilometer
Qp = 175*A(A>=10).^0.5

bar(Qp,0.5)
set(gca,'Xticklabel',Basin_name(A>=10))

xlabel('Sub-basin Names')
ylabel('peak discharge(m^3/s)')
title('Mayer method for peak discharge estimation')