d(CO)/d(t) = (a+b*sin(pi*t/6)+F-v0*CO)/V # v0 = if(t<72) then (1.67E12) else (0) #ft3/h a = 35000 #lbmol CO/h b = 30000 #lbmol CO/h F = 3412.31 #lbmol CO/h V = 4*10^13 #ft3 pi = 3.14159265359 # t(0)=0 CO(0)=2E-10 t(f)=100

How do I turn this Polymath code into Matlab code? Thanks!