d(CO)/d(t) = (a+b*sin(pi*t/6)+F-v0*CO)/V # v0 = 1.67*10^12 #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)=72

transform the following polymath code into matlab code