매틀랩이나 c++로 연립2차 미분방정식을 풀어야 하는데

대략 x''= a(x'^2+y'^2)*cos(arctan(y'/x'))+b(x'^2+y'^2)*sin(arctan(y'/x'))

 y'' = c(x'^2+y'^2)*cos(arctan(y'/x'))+b(x'^2+y'^2)*sin(arctan(y'/x'))

이 식 코드랑 그래프를 구해야 되는데 어떻게 해야하는지 아시나요?