できませんでした。
オリジナル
wolframで
できませんでした。できるような気がしました。
微分方程式の解
y(x) = 1/8 (-x - cot^(-1)(1/3 tan(1/3 (9 x + π))))
?点検
sympyで
できませんでした。できるような気がしました。
from sympy import *
x = Symbol ('x',real=True)
y = Function('y')(x)
ans=dsolve(Eq(diff(y),cos(x+8*y)**2),ics={y.subs(x,pi/9):3*pi/36})
print("#",ans)
#plot (ans.rhs)
# Eq(x - atan(tan(x + 8*y(x))/3)/3 + pi*floor((-x - 8*y(x) - pi/2)/pi)/3 - atan(tan(2*pi/9)/3)/3 + 5*pi/9, 0)