sympyで(微分)「2019 神戸大学 前期日程 数学（理科系）1.(1)」 をやってみたい

読売新聞
https://www.yomiuri.co.jp/nyushi/sokuho/k_mondaitokaitou/kobe/mondai/img/kobe_zenki_sugaku_ri_mon.pdf
https://sokuho.yozemi.ac.jp/sokuho/k_mondaitokaitou/13/kaitou/img/kobe_zenki_sugaku_ri_kai.pdf
wolframalpha
log(x)/x
https://www.wolframalpha.com/input/?i=log(x)%2Fx
max{log(x)/x} = 1/e at x = e

fullscript.py

from sympy import *
x = symbols('x', real=True)
f = symbols('f', cls=Function)
f=log(x)/x
eq=f.diff(x,1)
v = solve(eq)
print("v=",v[0])
print(f.subs(x,v[0]))
# v= E
# exp(-1)