% Shock and Detonation Toolboox
% http://www.galcit.caltech.edu/EDL/public/cantera/html/SD_Toolbox/
%
% Generate plots and output files for a ZND detonation with the shock front
% traveling at speed U.

clear;clc;
display('demo_ZNDshk')

P1 = 100000; T1 = 300; U = 2000; 
q = 'H2:2 O2:1 N2:3.76';
mech = 'h2air_highT.cti';  

[gas] = znd_shk(1, U, P1, T1, q, mech, 'h2air');