function v = velocity(A)% velocity.m Compute velocity profile profile for a rxnmodel object%% INPUT/OUTPUT PARAMETERS%        A : rxnmodel object%% Copyright (c) by Raymond A. Adomaitis, 1998-2004unpack(A,'var')unpack(A,'param')v = vo*Ro*T/Tamb/RSF;