Lực hấp dẫn trong lòng Trái Đất

------------------------------------------------------------------------------
Nếu đã có giao điểm, bạn werty thử test với hàm y =a.x^n +b với b là đoạn OK , a và n tìm ra được nhờ vào các điều kiện:
y(x1) = y1 và y''''(x1)= tg(@).
Không biết còn hàm nào thoả mãn các điều kiện trên không?

Được haidelft sửa chữa / chuyển vào 12:34 ngày 20/10/2006

Còn một điều kiện biên nữa là tại vị trí trên mặt đất thì y`` bằng vô cùng, do đó tớ nghĩ hàm đa thức (n nguyên) không phải là nghiệm. Tớ đã thay hàm y = ax^n + b vào phương trình vi phân (6) thì thấy không thỏa.
Tớ post đồ thị y và y` của quỹ đạo nhanh nhất mà tớ tìm được (ứng với trường hợp alpha = pi / 4) để mọi người cùng tham khảo:


Được werty98 sửa chữa / chuyển vào 18:36 ngày 20/10/2006

Post luôn các file mô phỏng trên Matlab:
http://users.tpg.com.au/adslcvn/misc/diff_func1.m
http://users.tpg.com.au/adslcvn/misc/simgravity2.m