The most computationally intensive portions in these m-files utilize the MATLAB Partial Differential Equation Toolbox to calculate self-consistent electronic wavefunctions, energies, densities, and band-bending diagrams for hexagonal and triangular core-shell nanowires. The flowchart depicted in Fig. summarizes the overall algorithmic processes within HADOKEN, which are described in extensive detail in the following paragraphs. The set_input_parameters.m and set_doping_density.m routines allow the user to input parameters specifying the desired core/shell side lengths, conduction band edge energies, effective masses, dielectric constants, mesh resolution, and doping density function, respectively. These quantities are then used by the main_scf_dirichlet.m and main_scf_neumann.m routines, which initiate the HADOKEN code with either Dirichlet or Neumann boundary conditions for V( x, ỹ). As sample input, the material parameters for GaN and Al 0.3 Ga 0.7 N have been provided in set_input_parameters.m in the global variables vector_of_V0, vector_of_masses, and vector_of_eps, which represent the conduction band discontinuities, effective masses, and dielectric constants of each of the nanowire regions, respectively. Researchers interested in other material compositions can simply replace these numerical values in the set_input_parameters.m file to enable self-consistent simulations for other materials.