% task3a
%
%Makes a plot of the relationship of discharge versus channel slope for
%values of slope varying from 0 to 0.2 (about 12 degrees), all other
% variables constant. Answers the question, How does discharge vary with channel slope?
%
%Author:  Wendy Bohon with significant help from David Raleigh, Feb. 26,
%2009.


function task3a(width , depth , roughness)


% increment slope by 0.01 from 0 to .2
slopes = [0:.01:.2]

% call task1 function for velocity/discharge
[discharges velocity hydralicRadius] = task1(width , depth , roughness , slopes , 0);

%This plots shows how discharge varies with varying channel slope.
figure(2)
plot(slopes , discharges , '^b-');
title('Discharge v/s Varying Channel Slope')
xlabel('Slope')
ylabel('Discharge')
print('dischargeSlope.jpeg')