f(x) = sin(x^2)
g(x) = cos(x+sin(x)*3)
h(x) = f(g(x))

g = plot(f, (x,0,3)) + plot(g, (x,0,3), color='red') + plot(h, (x,0,3), color='green')
g += points([(a,f(a)) for a in [0,.1,..,3]], pointsize=15, color='black')
g

︠c9188bd6-76ad-4cb2-9806-da3e036463d5︠