clc
clear
% Horton formula
Basin_name = {'A';'B';'C';'D';'E'};
A = [12,15,28,91,4]; % area in Square kilometer
T = [2,5,15,25,50,100];% Return periods
u = 1;
for t = T
    Qp(:,u) = 71.2*t.^0.25*A.^-0.5.*A;
    u = u+1;
end
bar(Qp)
set(gca,'Xticklabel',Basin_name)
legend((num2str(T')))
xlabel('Sub-basin Names')
ylabel('peak discharge(m^3/s)')
title('Horton method for peak discharge estimation in different return periods')