import sympy as sym
import numpy as np
from fractions import Fraction
from IPython.display import display
sym.init_printing()
t,x = sym.symbols("t x",positive=True)
y = abs(t-x)/(1+t**2)
f = sym.integrate(y,(t,0,1))
display("積分した結果をそのまま出力",f)
fdx=2*x*np.arctan(x)-np.log(1+x**2)-np.pi/4*x+1/2*np.log(2)
tana=sym.solve(fdx)
display(tana)