Xét phương trình vi phân y′=−ky+Asin(ωx) với điều kiện đầu y(0)=y0 (mô hình phân rã có ngoại lực dao động). Dùng phương pháp Euler cải tiến (Heun) với bước h: y~=yn+hf(xn,yn),yn+1=yn+2h[f(xn,yn)+f(xn+1,y~)]. Cho k,A,ω,y0,h và số bước n, hãy tính yn.
Một dòng gồm các số k A ω y0 h n (n nguyên).
0.1≤k≤5, −5≤A≤5, 0.1≤ω≤5, −5≤y0≤5, 0.001≤h≤0.1, 1≤n≤1000.
In ra yn, làm tròn 6 chữ số sau dấu chấm thập phân.
Ví dụ:
Đầu vào:
1 2 1 0 0.01 100
Đầu ra:
0.669052
Đang tải editor...