function bresenhamStep(x1,x2,y1,y2) dX=x2-x1; dY=y2-y1; deltaX=abs(dX*2); deltaY=abs(dY*2); grad = dY/dX; plot(x1,y1,'ro'); plot(x2,y2,'ro'); plot([x1,x2],[y1,y2],'b-'); drawnow(); %case1 if(deltaX>deltaY&&dX>0&&grad>0) D = deltaY - deltaX…