clc
clear
% Inglis formula
Basin_name = {'A';'B';'C';'D';'E'};
A = [12,15,28,91,4]; % area in Square kilometer
Qp = 124*A./sqrt(A+10.24)

bar(Qp)
set(gca,'Xticklabel',Basin_name)

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