matlabで次の微分方程式を解こうとしています。(これらは、フーフィー ポリアコフ モノポール ansatz の yang-mills-higgs ラグランジアンから得られた方程式です)。これは私の関数ファイルです。2 つの変数 h と k と、変数 t に対するそれらの導関数があります。私の x(1)=h、x(2)=k、x(3)=dh\dt、x(4)=dk\dt. すべての関数の初期値は 0 です。
function xprime = monopole( t,x )
%UNTITLED Summary of this function goes here
% Detailed explanation goes here
xprime(1)=x(3);
xprime(2)=x(4);
xprime(4)=(1/(t.^2)).*((x(2).^2)-1).*x(2) + 4.*(x(1).^2).*x(2);
xprime(3)=(2/(t.^2)).*(x(2).^2).*x(1)-(1-(x(1)).^2).*x(1)-(2/t).*x(3);
xprime=xprime(:);
end
次のコードを実行すると >
> t0=0;
>> tf=10;
>> x0=[0 0 0 0];
>> [t,s]=ode45(@monopole,[t0,tf],x0);
>> plot(t,s(:,1));
私は何も得ていません。グラフ ウィンドウが表示されますが、何も含まれていません。この方程式には解があるはずです。点線の曲線は、1 から始まる曲線が k で、0 から h が得られる曲線です。
私の間違いは何ですか?