simulation of greisen–zatsepin–kuzmin limit