各位板友好，我的問題大概是這樣
我想利用4th order Runge-Kutta方法追蹤帶電粒子在靜電磁場中的軌跡，這個靜電磁場
無法用解析方法表示，所以我用Interpolation分別建立電場和磁場的內插函數，再來按照
4th order Runge-Kutta處理兩個couple的微分方程的方式寫一個函數rk4，然後用
NestWhile執行rk4
4th order Runge-Kutta
https://i.imgur.com/qSgdghB.png
其中efr、efz、Br和Bz都是電磁場的內插函數
https://i.imgur.com/raE1PUy.png
不過計算的粒子數量一多就需要耗費不少時間，我覺得應該是我的程式寫得不太好，想請
問各路高手該如何提升計算效率，compile rk4對計算時間沒有顯著影響，我目前覺得問
題應該是出在內插函數算太慢，這有辦法處理嗎? Interpolation後的函數應該不用再
compile吧? 還是說rk4本身就不應該那樣寫? 感恩感恩
compile rk4 (沒影響)
https://i.imgur.com/szV0baY.png