boom | Python Fiddle

from __future__ import division
from scipy.optimize import fsolve
from math import exp
from numpy import linspace
from numpy import array 
import numpy
from matplotlib import pyplot as plot 
from scipy.integrate import odeint 
import matplotlib.pyplot as plt
from math import sqrt
from math import log10
from math import pi

# Input

tubes = 200
R = 8.314
Fa0 = (400/3600)/200 # kmol/s
P0 = 200 # kPa
T0 = 423.15 # K
Tu = 298.15 # K
Ca0 = (P0/(R*T0)/1000)
Q = Fa0/Ca0
Qtube = Q/200

# Parameters
d =  0.03
k0 = 700
E = 40000
Wt = 700 # kg
K = 2 # kPa
Cp = array([500,250],dtype=float)
Hr = -47500 # J/mol
U = 200 # W/m^2 K
rho = 800 # kg/m^3
e = 1
a = (pi/4)*d**2

#Te = (T0 + ((-Hr)/(F0[0]*Cp[0] + F0[1]*Cp[1]))*x*F0[0])

def prf(var,w):
    x,T,P = var
    return[((1/Fa0)*k0*exp(-E/(R*T))*(1-x))*(T0*P)/(Qtube*(1 + e*x)*T*P0),
           #(-K*P0*T*(1 + e*x))/(P*T0),
           -K*(P0/P)*(T/T0)*(1 - e*x)/(rho*a),
            (U*(4/d)/rho)*(((Tu - T) + (-Hr)*(k0*exp(-E/(R*T))*Ca0*((1 + x)/(1 + e*x))*(T0*P)/(T*P0)))/(Fa0*Cp[0]))]
            
Int = array([0,T0,P0],dtype=float)
wspan = linspace(0,Wt,100)
Ci = odeint(prf,Int,wspan)

plt.plot(wspan, Ci[:,0])
plt.show()
print Ci[:,0]

Browser Version Not Supported

Due to Python Fiddle's reliance on advanced JavaScript techniques, older browsers might have problems running it correctly. Please download the latest version of your favourite browser.
Chrome 10+ Firefox 4+ Safari 5+
IE 10+

Let me try anyway!

Python Fiddle

Python Cloud IDE