اولا من کاهش مرتبه دادم یعنی dx/dy رو برابر p گرفتم پس باید یک دستگاه معادلات را حل کنیم
p'=83p+px^2.5+sin(px)
p=y/x+c
پس ابتدا تابع زیر را تعریف میکنیم:
function udot=equ(x,u)
udot=[u(2)
83*u(2)+23*(u(1)*u(2))^(2.5)+sin(u(1)*u(2))];
سپس:
[x,u]=ode45(@equ,[0 2],[3;1])