clear all
close all
%Inserimento dati
theta0=0;        %Frazione di ricoprimento a t=0
time=[0 2];      % tempi
ka=0.076;        % atm-1s-1
kd=0.095;        %s-1
p=[10,7.5,5];
%Isoterme di adsorbimento a diverse p
f1=@(t,theta) ka*p(1)*(1-theta)-kd*theta;
[t1,theta1]=ode45(f1,time,theta0);
plot(t1,theta1,'b')
axis([0 2 0 1]);
xlabel('tempo / s'),ylabel('theta')
hold on

f2=@(t,theta) ka*p(2)*(1-theta)-kd*theta;
[t2,theta2]=ode45(f2,time,theta0);
plot(t2,theta2,'r');

f3=@(t,theta) ka*p(3)*(1-theta)-kd*theta;
[t3,theta3]=ode45(f3,time,theta0);
plot(t3,theta3,'g');
legend('p=10 atm', 'p=7.5 atm', 'p=5 atm')



