clear all

a = 1;
L = 0:0.01:(a/2);

Vol = zeros(1,length(L));

for i=1:length(L)
    Vol(i) = ((a-2*L(i))^2)*L(i);
    V2(i) = L(i)^3;
end


plot(L,Vol,'r-','linewidth',2)
xlabel('Comprimento L (m)')
ylabel('Volume da Caixa (m^3)')
grid on
hold on
plot(L,V2,'b--','linewidth',2)
hold off
legend('Volume','L^3','location','northwest')