clear all
clc

TEnd = input('Please give a end time: ')
k1= 1;
k2= 2;
Ca(1)= 5;
Cb0= 0;
Cc0= 0;
T = 0:2:TEnd
[T, Ca]= ode(- k1*Ca(T) ,[0 T],5);
plot(T, Ca(T))