- Codice:
function f = portate(lv)
global x y z F
f(1) = F*z(1) - lv(1)*x(1) - lv(2)*y(1);
f(2) = F - lv(1) - lv(2);
end
che "chiamo" con
- Codice:
soluzioni = fsolve(@portate, x0)
lv(1) e lv(2) hanno dei significati fisici particolare ed è per questo che non posso permettermi di confonderle. In particolare lv(1) è la portata di un liquido e lv(2) è la portata di un vapore. Quando chiamo fsolve che mi restituisce il vettore "soluzioni" come faccio a capire se soluzione(1) = lv(1) o solzione(1) = lv(2) ?? Mi sono riuscito a spiegare?
Grazie in anticipo