Runge-Kuttaの4次積分法を2次元行列（x、y）のmatlabに実装しようとしています。マトリックスには、各ポイントの高さの値（float）が含まれています。アイデアは、マトリックスにパーティクルを配置し、方向ベクトルuおよびvによって「プッシュ」されるときのパスを監視することです。
matlabに組み込まれている「gradient」関数を使用して、オイラーの積分手法を実装しました。ただし、RK4統合の場合、勾配によって与えられる所定の勾配ではなく、可変点で4つの勾配を取る必要があります。
Matlabの勾配関数自体にはこの種のパラメーターがないと思います（Matlabのドキュメントから理解していることから）。これにはどのようなアプローチを取るべきですか？
任意の点で勾配を計算する独自の勾配関数を作成する場合、それは直進関数ではないため、どのように作成するかはよくわかりません。このアプローチについて何か提案はありますか？
ありがとう。