Bresenham画直线MATLAB实现

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;