MATLABの無名関数内で匿名関数を作成しよう

Question

こんにちは、これは私にとって非常に問題なのかもしれませんが、私は固執していますし、しばらくの間これをやろうとしています。

私はこのような何かをしようとしている：

a=0.5; 
ODE= @(x,y) ((-2*(a^2)*x*y)/(x^2+y^2)^2)/(1-((a^2)*(x^2-y^2)/(x^2+y^2)^2)); 


%The classical RK4 solution as a function 
k1 = ODE(x,y); 
k2 = ODE(x+0.5*dx,y+k1*0.5*dx); 
k3 = ODE(x+0.5*dx,y+k2*0.5*dx); 
k4 = ODE(x+dx,y+k3*dx); 
rk4= @(x,y,dx) y + 1/6*(k1 + 2*(k2+k3) + k4)*dx; 

だから私は、最初の無名関数としてODEを定義し、その後、私はODEからKさんとKさんから最後にRK4を定義しようとしています。

しかし、これは機能しません。どのような提案をするか？

Answer 1

​​〜k4は、機能として定義または使用されていません。 k1 = ODE(x,y)を宣言すると、ベクトルxとyが定義されていないODE関数を評価しようとします。あなたはおそらく何を言っているのですか：

k1 = @(x,y,dx) (ODE(x,y)); 
k2 = @(x,y,dx) (ODE(x+0.5*dx,y+k1(x,y,dx)*0.5*dx)); 
k3 = @(x,y,dx) (ODE(x+0.5*dx,y+k2(x,y,dx)*0.5*dx)); 
k4 = @(x,y,dx) (ODE(x+dx,y+k3(x,y,dx)*dx)); 
rk4= @(x,y,dx) (y + 1/6*(k1(x,y,dx) + \ 
       2*(k2(x,y,dx)+k3(x,y,dx)) + \ 
       k4(x,y,dx))*dx);