clc
clear
% Fanning formula
Basin_name = {'A';'B';'C';'D';'E'};
A = [12,15,28,91,4]; % area in Square kilometer
Qp = 2.64*A.^0.8

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

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