Code: Selecteer alles
from sympy import *
V, k, R, beta, g, m, V=symbols("V, k, R, beta, g, m, V")
V = k*R**2*( (6-4*cos(beta)-2*sin(beta))**0.5-1)**2/2+m*g*(R-R*sin(beta))
V = V.subs(g, 9.81).subs(R, 0.2).subs(k, 180)
dVdbeta=diff(V,beta).doit()
dVdbeta=dVdbeta.subs(beta, 33*pi/180)
dVdbeta.evalf()
result = solve(dVdbeta)
print(result)
Dat levert ook het juiste antwoord - ervan uitgaand dat de resultaten van wnvl1 en mij correct zijn.