for i=1:10
% solve the system of ODEs that are specified in bistability.m
[t,y] = ode45('bistabilitymin', [t0 tfinal], y0);
EE(i) = y(:,3);
end;