substituiere x+y=z, 1+y'=z' , y'=z'-1
alles einsetzen:
z'-1=tan^2(z)
z'=tan^2(z)+1=1/cos^2(z) , separierbare DGL
dz*cos^2(z)=dx
integrieren, links mit partieller Integration
1/4(2z+sin(2z))=x+const
2z+sin(2z)=4x+const
Rücksubstituieren:
2(x+y)+sin(2(x+y))=4x+C
2y+sin(2(x+y))=2x+C
AWB einsetzen → C=π
Explizit nach der Funktion y kann man nicht auflösen.