F0=40000; UAS=8
zA0=0.05; zB0=0.2; P=101.3; DA=(1-1-1)/1
Cp=31.4; Q=-142400; Tw=673; RB=1100
pA=zA0*(1-x)*P/(1+DA*zA0*x)
pB=(zB0-zA0*x)*P/(1+DA*zA0*x)
k=22400*exp(-75370/8.314/T)
r=-k*pA*pB
x'=-r*RB/(zA0*F0)
T'=(Q*r*RB-UAS*(T-Tw))/(Cp*F0)
x # 0
T # 673
INTEGRAL V[0,2]  step 0.05
TREND x, T
output1 V,x,T step 0.05

計算結果例　F0=40000; UAS=8e6
---------------------------------------
   V          1:x          2:T         
---------------------------------------
 0            0            673.0000    
 0.10000000   0.2320226    712.8071    
 0.2000000    0.5260068    743.3337    
 0.3000000    0.7411176    745.1796    
 0.4000000    0.8430365    727.5061    
 0.5000000    0.8902307    709.4480    
 0.6000000    0.9160088    696.4883    
 0.7000000    0.9325017    688.1395    
 0.8000000    0.9442821    682.9569    
 0.9000000    0.9533046    679.7632    
 1.000000     0.9605165    677.7752