clc
clear
% Dickens formula
A = 150; % area in Square kilometer
C1 = 14:28;% for mountain areas
Qp_mountain = C1*A.^0.75;% peak discharge in mountain areas

subplot(2,1,1)
plot(C1,Qp_mountain,'-*r')
xlabel('C  coefficient  range for mountain areas')
ylabel('peak discharge(m^3/s)')
title(['Dickens method for peak discharge estimation in mountain area',...
    ' Watershed Area: ',num2str(A)])

C2 = 2.8:.1:5.6;% for plain areas
Qp_Plain = C2*A.^0.75;%  peak discharge in plain areas
subplot(2,1,2)
plot(C2,Qp_Plain,'-*b')
xlabel('C  coefficient  range for plain areas')
ylabel('peak discharge(m^3/s)')
title(['Dickens method for peak discharge estimation in plain area',...
    ' Watershed Area: ',num2str(A)])
