Application of low intensity electric fields to interfere with tumor growth is being increasingly recognized as a promising new cancer treatment modality. Intratumoral modulation therapy (IMT) is a developing technology that uses multiple electrodes implanted within or adjacent tumor regions to deliver electric fields to treat cancer. In this study, the determination of optimal IMT parameters was cast as a mathematical optimization problem, and electrode configurations, programming, optimization, and maximum treatable tumor size were evaluated in the simplest and easiest to understand spherical tumor model. The establishment of electrode placement and programming rules to maximize electric field tumor coverage designed specifically for IMT is the first step in developing an effective IMT treatment planning system. Finite element method electric field computer simulations for tumor models with 2 to 7 implanted electrodes were performed to quantify the electric field over time with various parameters, including number of electrodes (2 to 7), number of contacts per electrode (1 to 3), location within tumor volume, and input waveform with relative phase shift between 0 and 2π radians. Homogeneous tissue specific conductivity and dielectric values were assigned to the spherical tumor and surrounding tissue volume. In order to achieve the goal of covering the tumor volume with a uniform threshold of 1V/cm electric field, a custom least square objective function was used to maximize the tumor volume covered by 1V/cm time averaged field, while maximizing the electric field in voxels receiving less than this threshold. An additional term in the objective function was investigated with a weighted tissue sparing term, to minimize the field to surrounding tissues. The positions of the electrodes were also optimized to maximize target coverage with the fewest number of electrodes. The complexity of this optimization problem including its non-convexity, the presence of many local minima, and the computational load associated with these stochastic based optimizations led to the use of a custom pattern search algorithm. Optimization parameters were bounded between 0 and 2π radians for phase shift, and anywhere within the tumor volume for location. The robustness of the pattern search method was then evaluated with 50 random initial parameter values. The optimization algorithm was successfully implemented, and for 2 to 4 electrodes, equally spaced relative phase shifts and electrodes placed equidistant from each other was optimal. For 5 electrodes, up to 2.5cm diameter tumors with 2.0V, and 4.1cm with 4.0V could be treated with the optimal configuration of a centrally placed electrode and 4 surrounding electrodes. The use of 7 electrodes allow for 3.4cm diameter coverage at 2.0V and 5.5cm at 4.0V. The evaluation of the optimization method using 50 random initial parameter values found the method to be robust in finding the optimal solution. This study has established a robust optimization method for temporally optimizing electric field tumor coverage for IMT, with the adaptability to optimize a variety of parameters including geometrical and relative phase shift configurations.