for n = 1 : npassos;

    //Calculando condições de contorno
    Q(1) = Q(3);
    Q(N+3) = Q(N+1);

    //Calculo do Upwind de Primeira Ordem
    for i = 2 : N + 2
        Qmais1(i) = Q(i) - C*(Q(i) - Q(i - 1));
    end

    //Atualizando t e u
    t = t + dt;
    Q = Qmais1;
    plot(x,Q,'bo-');
    sleep(dt);

end