function Q5_code

    t = 0:0.01:600;
    k = 0.0073;
    A0 = 10;
    
    y = (A0^2*k*t)./(1 + A0*k*t);
    
    figure(1)
    plot(t,y,'LineWidth',3)
    xlabel('time')
    ylabel('concentration, mol/L')
    title('Quantity of A that has Reacted, where A+B -> AB, A0 = B0 = 10')
    grid on
end