Estimados, me han hecho notar (con razón) que del código presentado por mí no queda claro cuál sería bp. Dice "bp=???;" pero luego no hay más referencias a ese valor.
El bp se usa para calcular tau2 (que también está comentado), y corresponde al xB en la base actual (bp=xB). Recuerden que si x=[xB; 0], entonces B*xB=b.
Como asumimos además (caso no-degenerado) que xB>0, los tau2 se calculan buscando el paso que debo dar para anular la i-ésima variable de la base, siempre que eso sea posible. Es durante esos cálculos que vuelve a aparecer bp.
Mañana también será clase de consulta.
Saludos.-