simulation for j =2, k, l