Dismiss Notice
Join Physics Forums Today!
The friendliest, high quality science and math community on the planet! Everyone who loves science is here!

MATLAB Rate of Reaction in matlab

  1. Apr 3, 2013 #1
    Can Anybody help to code this in matlab

    A->B --- A is Sbbstrate here and B is Product
    Y[/a]=stiochimetric
    Initial values are A[0] and B[0]
    Q= reaction quotient

    R=active if Q<k ----- R is reaction
    R=Stops if Q=k
    t=? ------ t is the time

    d[A][/dt]= -kA Conc_A(t,A[0],K) and it return A
    A=A[0] e[-k][t-t[0]]

    And for B
    ΔB = YΔA
    ΔB=B-B[0] = (A[0]-A)Y Conc_B (B[0],Y,A[0],A) gets A and Return B
    B=B[0] + (A[0]-A)Y

    t=[0.001 0.002 0.003...... 1ms ]

    A=Conc_A(t)
    B=Conc_B(t)

    Plot(t,A,B)
     
  2. jcsd
  3. Apr 4, 2013 #2
    Please do reply ?
     
  4. Apr 5, 2013 #3
    look up RK45
     
  5. Apr 9, 2013 #4
    I think the Matlab command is ode45 for an RK4-type algorithm.
     
Share this great discussion with others via Reddit, Google+, Twitter, or Facebook




Loading...