Mixed model with quadratic term and relative 95%CI: PLOT

I performed a mixed model with a quadratic independent variable.

Now I would like to plot my model with this quadratic term (including 95% CI) and related experimental points.

I tried using the formula: plot(mdl) (https://it.mathworks.com/help/stats/linearmodel.plot.html) that sounds good because it plot the Scatter plot and add a plot of a linear regression model.

Nevertheless, it doesn't work in this mixed model.

How can I plot my data with a quadratic line?

I enclose my code

% S data set with 8 measures for each subject (random effect Sub)
% Performance: dependent variable
% Years2: quadratic term
% Year: simple term
lme1=fitlme(S, 'Performance~1+Years2+Years+(1|Sub)')

Actually, with this formula I was able to plot the regression line point by point and the relative scatter with a cycle but I'm not able to plot the 95%CI

Answer 1

[PRB]

Any suggestion?