// Zenbakizko Kalkulua Ingenieri Kimikoentzat
// OCW-EHU
// http://ocw.ehu.es/irakaskuntza-teknikoak/zenbakizko-kalkulua-ingeniari-kimikoentzat/Course_listing
// 5. gaia: Integrazioa
// Ariketa 5.4
// Asier Aranzabal


clc
t=0:2:22;
C=[.244 .201 .331 .145 .143 .102 .206 .334 .552 .438 .145 .114];
Q=[64 83 62 191 165 153 139 67 30 39 99 197]*60;

tcrit=12.5:.25:15.5;
Ccrit=interp1(t,C,tcrit,'spline');
Qcrit=interp1(t,Q,tcrit,'spline');
mcrit=Ccrit.*Qcrit;
masa_cont=inttrap(tcrit,mcrit);

printf('La masa total de contamintante derramado entre las 12:30 y las 15:30 asciente a %f g',masa_cont)

//emaitza 61.898671 g

//adierazpen grafikoa
ti=0:.25:22;
Ci=interp1(t,C,ti,'spline');
Qi=interp1(t,Q,ti,'spline');
plot(t,C.*Q,'o')
plot(ti,Ci.*Qi,'-')
plot(tcrit,Ccrit.*Qcrit,'-r')
xtitle('5_4ariketa','denbora, h','kutsatsaile kupura, g')