// Zenbakizko Kalkulua Ingenieri Kimikoentzat
// OCW-EHU
// http://ocw.ehu.es/irakaskuntza-teknikoak/zenbakizko-kalkulua-ingeniari-kimikoentzat/Course_listing
// 3. gaia: Ekuazio sistemak
// Ariketa 3.6
// Asier Aranzabal

clc

t=[0 0.0002 0.0012 0.0062 0.0313 0.2463 0.6332 1.1257 1.8355 2.4681 3.1247 3.78 4.5991 5.2544 5.9097 6.5651 7.4358 8.2594 9.2257 10];
alfa=[5.5000 5.4989 5.49282 5.46268 5.31395 4.2543 3.06531 2.36388 2.10041 2.14127 2.2496 2.3536 2.4512 2.5036 2.5384 2.5614 2.5792 2.5883 2.5943 2.5969];
beta_=[0.6818 0.6817 0.6811 0.6781 0.6638 0.5832 0.5312 0.5222 0.5351 0.5489 0.5602 0.5683 0.5747 0.5779 0.5799 0.5812 0.5822 0.5827 0.583 0.5832];

for i=1:length(beta_);
A=[5.5 2.6;0.9*beta_(i) 1-beta_(i)];
b=[alfa(i); -(0.75-2*beta_(i))];
x=gauss(A,b);

CA(i)=x(1);
CS(i)=x(2);
CR(i)=1-CA(i)-CS(i);
end

t=t'
CA=CA'
CR=CR'
CS=CS'

plot2d(t,[CA CR CS])
hl=legend(['CA';'CR';'CS']);
xtitle('3_6ariketa','denbora,s','kontzentrazioa, mol/l')