V = {
let outside_right = nj.zeros(Math.round((numx-1)/2))
let barrier = nj.ones(Math.round(barrier_width/dx)).multiply(barrier_height)
let well = nj.zeros(Math.round(well_width/dx))
let outside_left = nj.ones(numx - Math.round(well_width/dx) - Math.round(barrier_width/dx) - Math.round((numx-1)/2) ).multiply(barrier_height)
let V = nj.concatenate(outside_left,well,barrier,outside_right)
return V.reshape(numx,1)
}