It seems to be a matter of the order of commands.
x = 1:60;
y = 1/3.*x;
plot(x,y)
grid on
axis([0 60 0 20])
axis equal
will return
what you don't want, as it screws up your limits.
So rather use:
axis equal
axis([0 60 0 20])
and it is alright: