clear all theta = 5; %Declaring the constants n=8; t = [1:1000]; %Time span nutrient = [0.01: 0.01: 10]; for t=1:1000 %Generating the data y(t) = (nutrient(t))^n / ((theta)^n + (nutrient(t))^n); end; plot(nutrient,y) xlabel('time'); ylabel('Growth'); title('Graph to show the growth of bacteria')