F=BIl simulation