کمک فوری در حل معادله دیفرانسیل مرتبه دو بروش رانگ کاتا نیشترم (نرم افزار متلب)

[Hamed_Saha]

سلام دوستان.
این پروزه ریاضی من هست که تا سه روز دیگه باید حلش رو از طریق متلب ارسال کنم, اگه کسی از دوستان بتونه کمکم کنه ممنون میشم.
من توی حل متلب از دستور ode45 که حل میکنم وقتی به گام 0.1 میرسم واگرا میشه و به 10 به توان 24 میرسه!
اینم بگم که خیلی وارد نیستم, کسی میتونه کمکم کنه فقط 3 روز دیگه فرصت دارم, اینم سوال:

مقدار y در معادله زير را به كمك روش رانگ – كوتا – نيشتروم درx=2 )با استفاده از برنامه نويسی در محيط (MATLAB بدست آوريد.

d[SUP]2[/SUP]y/dx[SUP]2[/SUP] = 83 dy/dx + 23√y[SUP]5[/SUP] + sin y

y(0) = 1
dy/dx (x=0) = 3
h= 0.05

اگه کسی بتونه کمکم کنه ممنون میشم تماس بگیره:

saha2alma@gmail.com
​

 

[حــامد]

اولا من کاهش مرتبه دادم یعنی 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])