Skip to content

Instantly share code, notes, and snippets.

@lxm1117

lxm1117/ode system

Last active April 14, 2017 19:08
Show Gist options

  • Select an option

    Learn more about clone URLs
  • Save lxm1117/b039a742d5e6e2ed544dfb91966515df to your computer and use it in GitHub Desktop.

Select an option

Learn more about clone URLs
Save lxm1117/b039a742d5e6e2ed544dfb91966515df to your computer and use it in GitHub Desktop.
Download ZIP
matlab ode45 rxns
Raw
gistfile1.txt
rxn.m
function dydt = rxn( t,y )
%UNTITLED Summary of this function goes here
% Detailed explanation goes here
k1=260;k1_r=4000;
dydt1=-k1*y(1)*y(2)+k1_r*y(3);
dydt2=-k1*y(1)*y(2)+k1_r*y(3);
dydt3=k1*y(1)*y(2)-k1_r*y(3);
dydt=[dydt1;dydt2;dydt3];
end
////////
t=0:0.01:1;
y_init=[1000/602/0.125,800/602/0.125,0];
y_init=y_init';
[t,y]=ode45(@rxn,t,y_init)
////////
function dydt = rxn( t,y )
%UNTITLED Summary of this function goes here
% Detailed explanation goes here
k1=260;k1_r=4000;
k2=9.6;k2_r=70;
dydt1=-(k1+k2)*y(1)*y(2)+k1_r*y(3)+k2_r*y(4);
dydt2=-(k1+k2)*y(1)*y(2)+k1_r*y(3)+k2_r*y(4);
dydt3=k1*y(1)*y(2)-k1_r*y(3);
dydt4=k2*y(1)*y(2)-k2_r*y(4);
dydt=[dydt1;dydt2;dydt3;dydt4];
end
/////////
t=0:0.1:100;
y_init=[1000/602/0.125,800/602/0.125,0,0]'
[t,y]=ode45(@rxn,t,y_init);
Raw
ode system
%% no explicit expression
syms x(t) A(t) AB(t)
kf=1;kb=1;
ode1=diff(AB)==kf*A*x-kb*AB;
ode2=diff(A)==-kf*A*x+kb*AB;
ode3=diff(x)==-kf*A*x+kb*AB;
odes=[ode1;ode2;ode3];
Raw
rxn_odes with sin(t) input
function dydt = rxn( t,y )
%UNTITLED Summary of this function goes here
% Detailed expl anation goes here
% x(t) + A <-> AB
% y(1)=AB; y(2)=A; y(3)=x(t)
k1=260;k1_r=0;
dydt1=k1*y(2)*y(3)-k1_r*y(1);
dydt2=-k1*y(2)*y(3)+k1_r*y(1);
dydt3=-k1*y(2)*y(3)+k1_r*y(1)+10*sin(2*pi*t);
dydt=[dydt1;dydt2;dydt3];
end
///
Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment