xdot=freeconv(t,x)

% Author: Housam Binous

% Free Convection Past a Vertical Flat Plate

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

% Email: [e-mail]

function xdot=freeconv(t,x)

global Pr

% derivation of the governing equations can be found in Deen's book
% Analysis of Transport Phenomena, OPU, 1998.

xdot(1)=x(2);
xdot(2)=x(3);
xdot(3)=-3*x(1)*x(3)+2*x(2)^2-x(4);
xdot(4)=x(5);
xdot(5)=-3*Pr*x(1)*x(5);

xdot=xdot';