clc delta=0.6; // mm delta1=delta/3; theta=2/15*delta; disp("Displacement thickness =") disp(delta1) disp("mm") disp("Momentum thickness =") disp(theta) disp("mm")