%load_ext sympy.interactive.ipythonprinting


import sympy
from sympy import var,Eq
var('a:z')#普通の小文字のaからzは全てsympyのシンボル（変数）として扱うことにしました。
eq1=Eq(a*x**2+b*x+c,0)
eq1


eq2=eq1.subs({a:1,b:0,c:1})# 各係数に値を代入。
from sympy import solve
ans=solve(eq2)
ans



def print_eq(eq):
    from IPython.display import display, Math
    from sympy import latex
    display(Math(latex(eq).replace('\imath','i')))
    
    

print_eq(ans)


print ans


print eq1
print eq2


from sympy import atan, pi
y=2*x+1
exp0=atan(y**2/(2*y+10)**2)
exp1=sympy.integrals.Integral(exp0,(x,-10.0,pi))
print_eq(exp1)
print(exp1)
print(exp0)



def exp0_func(x):
    from numpy import arctan             # numpyではarctan, sympyではatan
    return arctan((2*x + 1)**2/(4*x + 12)**2)
_x=np.arange(-10.0,np.pi,0.1)
plt.plot(_x,exp0_func(_x))



from scipy import integrate
integrate.quad(exp0_func,-10.0,np.pi)


exp1.evalf()


from sympy.plotting import plot
plot(exp0,(x,-10.0,pi))


import sympy.physics.units as unit
unit.find_unit("avogad")[1]


unit.avogadro_number


unit.find_unit("planc")


unit.planck


unit.planck.evalf()


unit.find_unit(unit.pressure)


3.5*unit.psi