enter image description hereI want to use gradient ascent pulse engineering (GRAPE) to optimize the CZ, however, my program cannot work properly, it has no ideal result. In fact, it is the reproduce of Fig. 1(d) in article “Time-Optimal Two- and Three-Qubit Gates for Rydberg Atoms_-Jandura_Pupillo_-2022”.enter image description here. Here is my matlab code, I cannot find why it is wrong: