import sympy as sp
from IPython.display import display

r,w = sp.symbols('r w',positive=True)
t = sp.symbols('t',real=True)
x = sp.Function('x')(t)
y = sp.Function('y')(t)
x = r*sp.cos(w*t)
y = r*sp.sin(w*t)
xy=sp.Matrix([x,y])
display("運動体の位置は",xy)
v = sp.Matrix([sp.diff(x,t),sp.diff(y,t)])
display("速度はベクトル",v)
vn = v.norm()
display("速さは",sp.simplify(vn))
a = sp.Matrix([sp.diff(x,t,2),sp.diff(y,t,2)])
display("加速度はベクトル",a)
an = a.norm()
display("加速度の大きさは",sp.simplify(an))