ContourPlot[((1+((1.95023*^-19)x2)/((7.04882*^-6)-x2)^4)sh(.1y))/(cosh(.1y)+sqrt(1-(((1.95023*^-19)x2)/((7.04882*^-6)-x2)^4)(sh(.1y))^2))+((x2)-1.80002*^-6)/((x2)-7.04882*^-6), {x,-1,1},{y,-1,1}]