function?a（x0y0xendyend）
dx=（xend-x0）;
dy=（yend-y0）;
p=[];
p（1）=2*dy-dx;
twoDy=2*dy;
twoDyMinusDx=2*（dy-dx）;
x=[];
y=[];
x（1）=x0;
y（1）=y0;
for?i=2:20
????x（i）=x（i-1）+1;
????if?p（i-1）<0
????????p（i）=p（i-1）+twoDy;
???