# convert degree to radian; toradian(x_) := pi/180/degree*x;