Witam,
Mam równanie liniowe które wyglada dokladnie tak:
10^((-2)*ph1(i))*Ka + 10^(-ph1(i))*(-Ka*(C1(i)) - Kw) - Ka*Kw = 0
ph1, Ka oraz Kw sa dane ph1 zmienia sie co petle Ka i Kw to stałe, chce obliczyć C1(i) dla każdego ph1, które mam. Probowałem się oprzeć na rozwiązaniu podanym tutaj: http://4programmers.net/Forum/Newbie/188001-matlab_-_rozwiazywanie_rownania_z_jedna_niewiadoma, ale nie bardzo chce to działać po moich przekształceniach.
Kod tego wyglada tak:
for i=1:(length(ph1))
sym(C1(i));
f = 10^((-2)*ph1(i))*Ka + 10^(-ph1(i))*(-Ka*(C1(i)) - Kw) - Ka*Kw;
[C1(i)] = solve(f,C1(i));
end