library(plotly)
horse <- matrix(nrow = 100, ncol = 100)
for (j in 1:100) {
    for (i in 1:100) {
        horse[i, j] <- (i-50)^2 / 2 - (j-50)^2 / 3
    }
}
plot_ly(z = ~horse, type = "surface")