We present an inverse method to engineer uniform-intensity focal fields with arbitrary shape. Amplitude, phase, and polarization states, as adjustable parameters, are used to seek the desired focal fields in the non-iterative computational procedure. Our method can be applied to the cases with low and moderate numerical aperture (NA), in which case the feasibility and validity of our approach have been demonstrated in theory, simulation and experiment, respectively. For the case of higher NA, simulated results based on the Richards-Wolf diffraction integral are shown. We also made some discussions on the experiments with the higher NA. Our method should have wide applications in optical micro machining, optical trapping and so on.