line.m

关于质数的一些程序

function fun=line(x0,y0,x1,y1)
if(x0==x1) %vertical
    y=y0:0.01:y1;
    x=x0+(x1-x0)*(y-y0)/(y1-y0);
    plot(x,y);
else
    x=x0:0.1:x1;
    y=y0+(y1-y0)*(x-x0)/(x1-x0);
    plot(x,y);
end