Jawaban No 611.docx

%Jawaban no 611 (b2) global k1 k2 k3 Ca0=0.0199; Cr0=0; Cs0=0; k1=4.189*10^-3; k2=4.04*10^-3; k3=7.8*10^-3; t0=1; tf=1500; tspan=[t0 tf]; CaCrCs_vector=[Ca0 Cr0 Cs0]; [t,CaCrCs]=ode45(@persamaan3,tspan,CaCrCs_vector); t_Ca_Cr_Cs=[t CaCrCs] %Membuat grafik konsentrasi tiap komponen plot(t,CaCrCs(1:end,1),'b',t,CaCrCs(1:end,2),'g',t,CaCrCs(1:end,3),'r') title('Profil Konsentrasi Tiap Komponen') xlabel('waktu') ylabel('konsentrasi') legend('Ca','Cr','Cs') function dCaCrCsdt=persamaan3(t,CaCrCs_vector) global k1 k2 k3 Ca=CaCrCs_vector(1); Cr=CaCrCs_vector(2); Cs=CaCrCs_vector(3); dCadt=-k1*Ca; dCrdt=(k1*Ca)-(k2*Cr)-(k3*Cr*Cs); dCsdt=(k2*Cr)-(k3*Cr*Cs); dCaCrCsdt=[dCadt dCrdt dCsdt]'; end