## Not run:
# plot(c(0, 1), c(0, 1), type="n",
# xlab = "", ylab = "", main = "Segment Intersection")
# grid()
# for (i in 1:20) {
# s1 <- matrix(runif(4), 2, 2)
# s2 <- matrix(runif(4), 2, 2)
# if (segm_intersect(s1, s2)) {
# clr <- "red"
# p1 <- s1[1, ]; p2 <- s1[2, ]; p3 <- s2[1, ]; p4 <- s2[2, ]
# A <- cbind(p2 - p1, p4 - p3)
# b <- (p3 - p1)
# a <- solve(A, b)
# points((p1 + a[1]*(p2-p1))[1], (p1 + a[1]*(p2-p1))[2], pch = 19, col = "blue")
# } else
# clr <- "darkred"
# lines(s1[,1], s1[, 2], col = clr)
# lines(s2[,1], s2[, 2], col = clr)
# }## End(Not run)
Run the code above in your browser using DataLab