Вы можете использовать растровую диаграмму, чтобы прояснить, что функция увеличивается только от нуля. Вы также можете иметь отрицательные значения t и x, хотя это также будет симметрично с этой стороны.
foo <- function(x, t) x ^ 2 + 4 * sqrt(t * x) + t ^ 2
df <- data.frame(
x = rep(seq(0, 40, by = 1), times = 41),
t = rep(seq(0, 40, by = 1), each = 41)
)
df$z <- foo(df$x, df$t)
library(ggplot2)
ggplot(df) +
geom_raster(aes(x = x, y = t, fill = z)) +
scale_fill_viridis_c()