用matlab求解哈雷彗星在1986年2月9日到达了近日点（最接近太阳的点,取太阳为原点）,那时它的位置和速度分别为位置

这是从哪弄来的题目?你确定 运动方程没错?我做出的结果怎么不对啊,在近日点是圆弧,后面好像是一条直线,离太阳越来越远,不再回来
dt=0.01;u=4*pi^2;
t(1)=0;
x(1)=0.325514;y(1)=-0.459460;z(1)=0.166229;
vx(1)=-9.096111;vy(1)=-6.916686;vz(1)=-1.305721;
x(2)=vx(1)*dt+x(1);y(2)=vy(1)*dt+y(1);z(2)=vz(1)*dt+z(1);
r(1)=sqrt(x(1)^2+y(1)^2+z(1)^2);
i=1;
while t