jisi.m

来自「自己编的四节龙格琨塔程序」· M 代码 · 共 7 行

%Creat the function for jisi
%Let x(1)=x, x(2)=y

function dx=jisi(t,x,a,b,c)
s=sqrt((c-x(1))^2+(a*t-x(2))^2);
dx=[b*(c-x(1))/s;b*(a*t-x(2))/s];