Inductor Calculator

from math import *


density = 11340.0 #lead

breakdown_amps = 4e8 #1e9 with future tech

mu0 = 1.2566370614e-6

def inductor( area_wire, diameter_coil, turns ):
    current = area_wire * breakdown_amps

    #radius_wire = (area_wire/pi)**0.5 #circular
    radius_wire = (area_wire ** 0.5)*0.5 #square
    radius_coil = 0.5*diameter_coil - 0.5*radius_wire
    #square coil
    #length = 4.0*(area_coil**0.5) * turns
    #circular coil
    length = 2.0*pi*radius_coil * turns
    area_coil = pi * radius_coil*radius_coil

    mass = length * area_wire * density

    energy = mu0 * turns*turns * area_coil * current*current / (2.0*length)

    print("Mass:   %f\nEnergy: %f\nLength: %f,\nCurrent: %f"%(mass,energy,length,current))

inductor( 0.01*0.01, 0.1, 10 )
inductor( 0.0002, 1.0, 500 )