The erosion rate of a streambank is calculated by the excessshear stress equation: e=k(a-c) where e is the erosion rate in m/s,k is the erodibility coefficient in m^3/ N s, a is the appliedshear stress (Pa) and c is the critical shear stress (Pa). Theapplied shear stress is calculated using the equation a=den*g*d*swhere den is the density of water (1000 kg/m^3), g is theacceleration due to gravity (9.8 m/s^2), d is the depth of water(m) and s is the bed slope (m/m). The critical shear stress is afunction of the silt and clay content of the streambank and iscalculated using the following equation:c=0.1+0.1779(SC%)+0.0028(SC%)2-2.34E-5(SC%)3 where SC% is thepercent of the streambank that is silt and clay. The erodibilitycoefficient is calculated using the equation k=0.20*c^-0.5. Developa program in MATLAB that calculates the erosion rate of astreambank. The user should enter the following values: depth ofwater, bed slope and percent of silt and clay. The program shouldwrite out “The erosion rate is:â€, the rate calculated and theunits. Test the program using the following: d=1, s=.001, SC=35 andd=0.1, s=.004, SC=50. Document the erosion rate. Paste all inputand output code into this Word document.