Con lắc đơn không tuyến tính hóa thỏa θ′′=−Lgsin(θ) với g=9.8, điều kiện đầu θ(0)=θ0,θ′(0)=ω0. Đặt ω=θ′, hệ bậc 1: θ′=ω, ω′=−Lgsin(θ). Dùng phương pháp Runge-Kutta bậc 2 (trung điểm) với bước h: k1=F(un),k2=F(un+2hk1),un+1=un+hk2. Cho L,θ0,ω0,h và số bước n, hãy tính θn và ωn.
Một dòng gồm các số L θ0 ω0 h n (n nguyên, θ0 tính bằng radian).
0.2≤L≤3, −1.5≤θ0≤1.5, −2≤ω0≤2, 0.001≤h≤0.05, 1≤n≤2000.
In θn và ωn trên một dòng, cách nhau khoảng trắng, mỗi số làm tròn 6 chữ số sau dấu chấm thập phân.
Ví dụ:
Đầu vào:
1 0.5 0 0.01 200
Đầu ra:
0.496515 0.180968
Đang tải editor...