Code covered by the BSD License
% Combustion rate function function [x,dx] = sz_comb(th,xr,ths,thb) if th < ths x = xr; dx = 0; else x = xr + (1 - xr) * (1 - exp(-5 * ((th - ths)/thb)^3)); dx = (1-xr) * (15/thb) * (((th-ths)/thb)^2)*exp(-5*((th-ths)/thb)^3); end end
Contact us