clear all
close all
clc
 
a=8;b=15 
jstxb=[];jstyb=[];zscxb=[];zscyb=[];
d=0.1;dt=0.05;N=1/dt;t=0;
 
 
jstx=200;jsty=0;zscx=0;zscy=0;
while (sqrt((jstx-zscx)^2+(jsty-zscy)^2)>=30)
    
    t=t+dt
    jstx=jstx-b*dt*jstx/sqrt(jstx^2+(a*t-jsty)^2);
    jstxb=[jstxb,jstx];
    jsty=jsty+b*dt*(a*t-jsty)/sqrt(jstx^2+(a*t-jsty)^2);
    jstyb=[jstyb,jsty];
    zscy=a*t;
    zscyb=[zscyb,zscy];
end
figure(1);
zscxb=zeros(length(zscyb));
plot(jstxb,jstyb,'o',zscxb,zscyb,'.');
hold on
i=0;
jsty
pause
while (sqrt((jstx-zscx)^2+(jsty-zscy)^2)>=d)
    if mod(i,20)==0
    a=0.5*a
    b=1.1*b
    end
    t=t+dt;
    jstx=jstx-b*dt*jstx/sqrt(jstx^2+(a*t-jsty)^2);
    jstxb=[jstxb,jstx];
    jsty=jsty+b*dt*(a*t-jsty)/sqrt(jstx^2+(a*t-jsty)^2);
    jstyb=[jstyb,jsty];
    zscy=a*t;
    zscyb=[zscyb,zscy];
    i=i+1;
    jsty
end
zscxb=zeros(length(zscyb));
figure(2);
plot(jstxb,jstyb,'o',zscxb,zscyb,'.');
size(zscyb)