clc u=1.5; // Ns/m^2 v=0.5; // m/s H=0.02/2; // m t=-u*3*v/H; disp("The shear stress =") disp(t) disp("N/m^2") disp("It acts in the opposite direction to the flow.")