HowarthsTransformation.m file provides a framework for solving the problem.
y3y5_0 = HowarthsTransformation(rhofun, miufun, hw, M, Pr, Gamma, y3y5_0guess)
The default values of the parameters are:
rhofun = @(h) h^(-1);
miufun = @(h) h^(2/3);
hw = 2; M = 0; Pr = 0.7;
Gamma = 1.4; y3y5_0guess = [0.1;2];
as described in the link:
https://en.wikipedia.org/wiki/Blasius_boundary_layer#Howarth_transformation
Example 1
y3y5_0 = HowarthsTransformation() %using only default values;

Example 2
y3y5_0 = HowarthsTransformation(@(h)1/(2+tanh(h)), @(h)h^0.7, 5, 0.7, 1.4);