Entonces queremos encontrar al más pequeñox (por más pequeño, nos referimos al más pequeño‖) tal queAx=b, que también es el problema de mínimos cuadrados. def forward_subst(L,b): # T...Entonces queremos encontrar al más pequeñox (por más pequeño, nos referimos al más pequeño\|x\|^2) tal queAx=b, que también es el problema de mínimos cuadrados. def forward_subst(L,b): # This function solves $L x= b$ when $L$ is the lower-trigular matrix n = L.shape[0]; x = np.zeros(n); for i in range(n): x[i] = b[i] for j in range(i): x[i] = x[i] - L[i,j]*x[j] x[i] = x[i]/L[i,i] return np.matrix(x).T
