// Zenbakizko Kalkulua Ingeniari Kimikoentzat // OCW-EHU // http://ocw.ehu.es/irakaskuntza-teknikoak/zenbakizko-kalkulua-ingeniari-kimikoentzat/Course_listing // 4. gaia: Ekuazio Diferentzial Arruntak // Ariketa 4_5 // Asier Aranzabal function cprima=tanke_sistema(t,c) cprima=zeros(5,1); cprima(1)=-0.12*c(1)+0.02*c(3)+1; cprima(2)=0.15*c(1)-0.15*c(2); cprima(3)=0.025*c(2)-0.225*c(3)+4; cprima(4)=0.1*c(3)-0.1375*c(4)+0.025*c(5); cprima(5)=0.03*c(1)+0.01*c(2)-0.04*c(5); endfunction t=[0:2:100]'; c=ode([0;0;0;0;0],0,t,tanke_sistema); [t,c'] plot2d(t,c'); hl=legend(['c1';'c2';'c3';'c4';'c5']); xtitle('4_4ariketa','denbora,s','kontzentrazioa, mol/l'); emaitza=Euler(0,[0 0 0 0 0],2,100,tanke_sistema)