Beispiel #1
0
def c6p19():
    # Impedance
    def z(r, c, l, w):
        return np.sqrt(r**(-2) + (w * c - (w * l)**(-1))**2)**(-1)

    # Constants
    r, c, l, z_desired = 225., 0.6 * 10**(-6), 0.5, 100.

    # Args for solution
    f = lambda w: z(r, c, l, w) - z_desired
    x0, x1, dx = 210, 215, 10**(-4)

    return {
        'newton_raphson': num.newton_raphson(f, x0),
        'secant_method': num.secant_method(f, x0, x1),
        'modified_secant': num.modified_secant(f, x1, dx)
    }
Beispiel #2
0
def c6p19():
	# Impedance
	def z(r, c, l, w):
		return np.sqrt( r**(-2) + ( w * c - (w * l)**(-1) )**2 )**(-1)

	# Constants
	r, c, l, z_desired = 225., 0.6 * 10**(-6), 0.5, 100.

	# Args for solution
	f = lambda w: z(r, c, l, w) - z_desired
	x0, x1, dx = 210, 215, 10**(-4)

	return {'newton_raphson': num.newton_raphson( f,x0 ), 'secant_method': num.secant_method( f, x0, x1 ), 'modified_secant': num.modified_secant( f,x1,dx) }