//mass of air flow(in kg/min) ma=5; //mass of fuel flow(in kg/min) mf=0.5; //specific gravity of fuel Pf=0.75; //pressure of air at inlet(in bar) p1=1; //temperature at inlet(in K) T1=300; //velocity coefficient Vc=0.8; //ratio of pressure drop of metering orifice to choke pd=0.8; //specific heat a constant pressure Cp=1005; //gas constant R=0.287; //coefficient of discharge of fuel Cf=0.6; //adiabatic constant y=1.4; //velocity at throat(in m/s) C2=100;