El algoritmo tiene un bug, desde la época de mi código en excel, que da error cuando se pretenden modelar dos (o n) componentes iguales, es decir, con igual Hcr e igual m(Hcr). El problema aparece al querer hallar los Msi porque queda un denominador que se hace nulo.
Lo que hay que hacer es, meter algún if o case que diga que cuando las n componentes son iguales (mismo Hcri y mismo mi(Hcri), los Msi son todos iguales entre sí y se calculan asì:
Msi = Mtip / (n * mi(Htip))
donde n es el número de componentes. mi(Htip) se calcula como se estaba haciendo para el sist de ecuaciones a resolver con linsolve. Lo que no se usa, es Hx (puede decir cualquier cosa o nada que da lo mismo).
Creo que así no debería haber más problemas. Verificar.
Me queda la duda por ej para el caso de 3 componentes, qué pasa si 2 son iguales entre sí y la tercera no. Quedará para seguir pensando, pero al menos la corrección de este caso es útil.
El algoritmo tiene un bug, desde la época de mi código en excel, que da error cuando se pretenden modelar dos (o n) componentes iguales, es decir, con igual Hcr e igual m(Hcr). El problema aparece al querer hallar los Msi porque queda un denominador que se hace nulo. Lo que hay que hacer es, meter algún if o case que diga que cuando las n componentes son iguales (mismo Hcri y mismo mi(Hcri), los Msi son todos iguales entre sí y se calculan asì:
Msi = Mtip / (n * mi(Htip))
donde n es el número de componentes. mi(Htip) se calcula como se estaba haciendo para el sist de ecuaciones a resolver con linsolve. Lo que no se usa, es Hx (puede decir cualquier cosa o nada que da lo mismo).
Creo que así no debería haber más problemas. Verificar.
Me queda la duda por ej para el caso de 3 componentes, qué pasa si 2 son iguales entre sí y la tercera no. Quedará para seguir pensando, pero al menos la corrección de este caso es útil.