function simulasi
clc, clear
yo = [0.1978495 0 0 1.582796];
tspan = [0:2:120];
[t,y] = ode23(@sub_simulasi,tspan,yo);
[t y]

plot (tspan,cA,'-*g', tspan,cMOH,'-+b', tspan,cME,'-*r',... 
tspan,cH2O,'-*c')
title ('grafik konsentrasi vs waktu');
xlabel (' t (menit)');
ylabel ('konsentrasi (mol/L)')
legend ('cA','cMOH','cME','cH2O')
%--------------------------------------------------------------------