You can manually generate wavefunction files for N, N+1 and N-1 states, and then use Multiwfn to plot Fukui function, please check the very detailed example in Section 4.5.4 of Multiwfn manual. A even more easier way is directly using the conceptual density functional theory (CDFT) module of Multiwfn to fully automate this process, this module has been systematically introduced in this book chapter: Chapter Realization of Conceptual Density Functional Theory and Info...
, corresponding usage example is given in Section 4.22 of Multiwfn manual.