
disp('Setting up vars for a Long Fin problem...');
syms Tx Tb Tinf;
syms h k t w A p r x;
syms Qdot;

m = sqrt((h*p)/(k*A));

condition = Qdot == sqrt(h*p*k*A)*(Tb-Tinf)
eqn = (Tx-Tinf)/(Tb-Tinf)==exp(-m*x)