clear all
close all
clc

% using inline functions
f=inline('A*cos(B*t)+C*t+D','t','A','B','C','D');
best=[1 2 3 4];
quad(@(t)f(t,best(1),best(2),best(3),best(4)),0,10);

% using ode45 with options and parameters
tspan=[0 10];
y0=[1 0 0];
options=odeset('RelTol',1e-3);
params=[1 2 3 4];
[T Y]=ode45('myfun2',tspan,y0,options,params);

colors=['r','g','m','k'];
figure(1)
hold on
for i=1:3
    plot(T,Y(:,i),colors(i))
end

init=[.2 .4 0];
k1=1;
k2=.3;
k3=.001;
k4=0;
opt=[];
[T2 Y2]=ode45('chem',tspan,init,opt,k1,k2,k3,k4);
figure(2)
hold on
for i=1:3
    plot(T2,Y2(:,i),colors(i))
end