data(triMesh)
meshP = triMesh$MeshPts$p
BoundConx = BoundCony = numeric(NROW(meshP))
BoundConx[1:NROW(meshP)] = BoundCony[1:NROW(meshP)] = 1
BoundConx[c(10, 11, 12)] = BoundCony[c(10, 11, 12)] = 0
bound = ApplyBC.2d(meshP, BoundConx, BoundCony)
Run the code above in your browser using DataLab