f(x):=1/(x-3); wxdraw2d(grid=true, yrange=[-50,50], explicit(f(x),x,2,4), color=red, line_type=dots, parametric(3,t,t,-50,50), parametric(t,0,t,2,4) );
from here
