资源描述：
《燃烧学第一次作业.doc》由会员上传分享，免费在线阅读，更多相关内容在应用文档-天天文库。
1、g=2.9对应系列二曲线programH2_Chained_Reactionimplicitnonerealtao0,CH0,tao,CH,tao1,CH1,dtao,taomax,CHmaxintegerii=0tao0=0CH0=0dtao=0.1tao1=tao0CH1=CH0taomax=100CHmax=2open(01,file='CH.dat')dowhile(abs(tao)<=taomax.and.abs(CH)<=CHmax)i=i+1tao=tao+dtaoCH=CH1+dtao*CHT(CH1)print*,i,tao,CHT(CH),CHwrite(01,*)tao
2、,CHT(CH),CHtao1=taoCH1=CHenddoclose(01)containsfunctionCHT(CH)implicitnonerealCHT,CHCHT=0,01+0.1*CHendfunctionendprogramH2_Chained_Reactiong=3.0对应系列四曲线programH2_Chained_Reactionimplicitnonerealtao0,CH0,tao,CH,tao1,CH1,dtao,taomax,CHmaxintegerii=0tao0=0CH0=0dtao=0.1tao1=tao0CH1=CH0taomax=100CHmax=2o
3、pen(01,file='CH.dat')dowhile(abs(tao)<=taomax.and.abs(CH)<=CHmax)i=i+1tao=tao+dtaoCH=CH1+dtao*CHT(CH1)print*,i,tao,CHT(CH),CHwrite(01,*)tao,CHT(CH),CHtao1=taoCH1=CHenddoclose(01)containsfunctionCHT(CH)implicitnonerealCHT,CHCHT=0.01endfunctionendprogramH2_Chained_Reactiong=3.1对应系列三曲线programH2_Chaine
4、d_Reactionimplicitnonerealtao0,CH0,tao,CH,tao1,CH1,dtao,taomax,CHmaxintegerii=0tao0=0CH0=0dtao=0.1tao1=tao0CH1=CH0taomax=100CHmax=2open(01,file='CH.dat')dowhile(abs(tao)<=taomax.and.abs(CH)<=CHmax)i=i+1tao=tao+dtaoCH=CH1+dtao*CHT(CH1)print*,i,tao,CHT(CH),CHwrite(01,*)tao,CHT(CH),CHtao1=taoCH1=CHend
5、doclose(01)containsfunctionCHT(CH)implicitnonerealCHT,CHCHT=0.01-0.1*CHendfunctionendprogramH2_Chained_Reactiont<30时，g=2.9；t>=30时，g=(50-CH)*0.1对应系列一曲线（前半段几乎覆盖了系列二曲线）programH2_Chained_Reactionimplicitnonereal::tao0,CH0,tao,CH,tao1,CH1,dtao,taomax,CHmaxinteger::ii=0tao0=0.CH0=0.dtao=0.1tao1=tao0CH1=C
6、H0taomax=100.CHmax=1.e30open(01,file='CH.dat')dowhile(abs(tao)<=taomax.and.abs(CH)<=CHmax)IF(tao<30)theni=i+1tao=tao+dtaoCH=CH1+dtao*CHT1(CH1)print*,i,tao,CHT1(CH),CHwrite(01,*)tao,CHT1(CH),CHtao1=taoCH1=CHELSEi=i+1tao=tao+dtaoCH=CH1+dtao*CHT2(CH1)print*,i,tao,CHT2(CH),CHwrite(01,*)tao,CHT2(CH),CHt
7、ao1=taoCH1=CHENDIFenddoclose(01)containsfunctionCHT1(CH)implicitnonerealCHT1,CHCHT1=0.01+0.1*CHendfunctionfunctionCHT2(CH)implicitnonerealCHT2,CHCHT2=0.01+CH*(0.1*CH-2)endfunctionendprogramH2_Chained_Reacti
