clear all;

global a alphdot m;

a = 1.019;

V_0 = 0.5;

x_0 = 2.7;

alphdot = 0.9;

m=1.2;

%solve the ODE equations

resource = [];
out1=[];
out2=[];

[t,x] = ode15s('ODE_resources',[0.01:0.1:7],[V_0; x_0]);
    
%out1=[out1 x(:,1)];       
%out2=[out2 x(:,2)];

% plot v vs t

t=t+1.5;

hold on;
plot(t,x(:,1));

xlabel('time (hours)'); ylabel ('Growth of bacteria in volume');
title('Graph to show the growth of bacteria');
axis([0 7 0 7]);