how to plot figure of expression with `in RootOf`?

i use ilaplace to get the output signal but how to plot the output signal `xt`

codes are as this

syms t s
gs=(5*s + 4)/(s^4 + 4*s^3 + 2*s^2 + 3*s);
us= laplace(sin(t + pi/6), t, s);

xs = us*gs;
xt = ilaplace(xs, s, t)

and the xt is

 
xt =
 
sin(t)/4 - (9*cos(t))/4 + (5*symsum(exp(r3*t)/(3*r3^2 + 8*r3 + 2), r3 in RootOf(s3^3 + 4*s3^2 + 2*s3 + 3, s3)))/4 + (7*3^(1/2)*symsum((r3*exp(r3*t))/(3*r3^2 + 8*r3 + 2), r3 in RootOf(s3^3 + 4*s3^2 + 2*s3 + 3, s3)))/12 + (2*3^(1/2))/3 + (35*symsum((r3*exp(r3*t))/(3*r3^2 + 8*r3 + 2), r3 in RootOf(s3^3 + 4*s3^2 + 2*s3 + 3, s3)))/4 - (3^(1/2)*cos(t))/4 - (5*3^(1/2)*symsum((r3^2*exp(r3*t))/(3*r3^2 + 8*r3 + 2), r3 in RootOf(s3^3 + 4*s3^2 + 2*s3 + 3, s3)))/12 - (9*3^(1/2)*sin(t))/4 + (95*3^(1/2)*symsum(exp(r3*t)/(3*r3^2 + 8*r3 + 2), r3 in RootOf(s3^3 + 4*s3^2 + 2*s3 + 3, s3)))/12 + (9*symsum((r3^2*exp(r3*t))/(3*r3^2 + 8*r3 + 2), r3 in RootOf(s3^3 + 4*s3^2 + 2*s3 + 3, s3)))/4
 
>> 

so how to plot `xt` with `t`?

[Walter Roberson]

With difficulty, especially in releases before R2019a. You either need to know a fair bit about MuPAD programming, or you need to find someone willing to create the icky MuPAD code for you.

[dcydhb dcydhb]

but what is the relationship between MuPAD and the plot?

[dcydhb dcydhb]

i just want to know if there is a way to plot `xt-t` if t=0:0.1:50 , and thanks a lot!

Answer 1

[Walter Roberson]

Tt=0:0.1:50;
Y=double(subs(xt, t, Tt)) ;
plot(tt, Y)