f=derivation1(a)

% Author: Housam Binous

% Free Convection Past a Vertical Flat Plate

% National Institute of Applied Sciences and Technology, Tunis, TUNISIA

% Email: [e-mail]

function f=derivation1(a)

x0=[0 0 a(1) 1 a(2)];

[t,Y]=ode15s(@freeconv,[0 100],x0);

% we need to verify the outer fluid conditions
% this is a split boundary condition problem

f(1)=Y(end,2);
f(2)=Y(end,4)^2+Y(end,5)^2;