Fout gevonden, nu 3867.13330195116.
Dan zitten we dicht in elkaars buurt.
Code: Selecteer alles
from sympy import *
import numpy as np
from matplotlib import pyplot as plt
m, gewicht_lucht, gewicht_par, Fcoulomb_lucht, Fcoulomb_par, rho_par, V, rho_metaal, epsilon_r_par, k, g, l, theta_lucht, theta_par, Q = symbols('m gewicht_lucht gewicht_par Fcoulomb_lucht Fcoulomb_par rho_par V rho_metaal epsilon_r_par k g l theta_lucht theta_par Q')
eq1 = Eq(g, 10) #N/kg
eq2 = Eq(k, 8.988*10**9) #N/kg
eq3 = Eq(epsilon_r_par, 2) #N/kg
eq4 = Eq(rho_par, 800) #N/kg
eq5 = Eq(theta_lucht, 70*pi/180) #N/kg
eq6 = Eq(theta_par, 60*pi/180) #N/kg
eq7 = Eq(Fcoulomb_lucht, k*Q**2/(l*sin(theta_lucht))**2)
eq8 = Eq(Fcoulomb_par, (k/epsilon_r_par)*Q**2/(l*sin(theta_lucht))**2)
eq9 = Eq(gewicht_lucht, V*rho_metaal*g) #N/kg
eq10 = Eq(gewicht_par, V*(rho_metaal-rho_par)*g) #N/kg
eq11 = Eq(tan(theta_lucht), Fcoulomb_lucht/gewicht_lucht) #N/kg
eq12 = Eq(tan(theta_par), Fcoulomb_par/gewicht_par) #N/kg
solve([eq1,eq2,eq3,eq4,eq5,eq6,eq7,eq8, eq9, eq10, eq11, eq12], dict=True)
[{Fcoulomb_lucht: 10178679290.909*Q**2/l**2,
Fcoulomb_par: 5089339645.45448*Q**2/l**2,
V: 97656.0488955274*Q**2/l**2,
epsilon_r_par: 2.00000000000000,
g: 9.81000000000000,
gewicht_lucht: 3704736286.03268*Q**2/l**2,
gewicht_par: 2938331614.30058*Q**2/l**2,
k: 8988000000.00000,
rho_metaal: 3867.13330195116,
rho_par: 800.000000000000,
theta_lucht: 1.22173047639603,
theta_par: 1.04719755119660}]