function [model] = create_func

 model.odefun = @gas_phase_ODE;

end


function dxdt = gas_phase_ODE(x,t)
  
  k_ = 0.16;
  v = [-2 1];
  
  r = k_*x(1)^2;

  dxdt = v'*r;

end

