clc;
clear all;
disp("displacement thickness")
disp("velocity distribution = u/U=y/d")
disp("i) the displacement thickness, d*:")
disp("d*=intergrate((1-u/U,U,0,d)")
disp("d*=d/2")
disp("ii) the momentum thickness, th")
disp("th = integrate(u(1-u/U)/U,u,0,d)")
disp("th=d/6")
disp("iii) The energy thickness de")
disp("de=integrate (u/U*(1-(u/U)^2),U,0,d")
disp("de=d/4")
disp("iv) The value of d*/th")
disp("d*/th =3")