@xzczd I followed C.E. suggestion to input these conditions
@zhk First, do you have any reason to believe that your answer is wrong. If so, why? Second, you have the boundary condition du/dx = -u[x], why would you not have du/dx = NeumannValue[-u[x, y], x==1] instead of NeumannValue[u[x, y], x==1] ?