Problem with ploting the solution

Make sure to change the boundary condition

to this

bord2 = 'near(x[0], 1)'

as I pointed out in the original post