ilaplace 不能在xy平面画图 请大神指教！

Wkkkkkkkkk

程序：
syms s t
disp('象函数Lf')
Lf = (s+251.3)/(s*s*s*s+2*160*pi*s*s*s+45600*pi*pi*s*s+20000*160*pi*pi*pi*s+10000*25600*pi*pi*pi*pi)
disp('Lf的原函数f')
f = ilaplace(Lf)
subplot(1,2,1)
ezplot(Lf)
subplot(1,2,2)
ezplot(f)

提示：
>> ilaplace
象函数Lf

Lf =

(s + 2513/10)/(s^4 + 320*pi*s^3 + (3865934085194845*s^2)/8589934592 + (6658547215630757*s)/67108864 + 817126680320561/32768)

Lf的原函数f

f =

8589934592*symsum((root(s5^4 + 320*s5^3*pi + (3865934085194845*s5^2)/8589934592 + (6658547215630757*s5)/67108864 + 817126680320561/32768, s5, k)*exp(root(s5^4 + 320*s5^3*pi + (3865934085194845*s5^2)/8589934592 + (6658547215630757*s5)/67108864 + 817126680320561/32768, s5, k)*t))/(2*(3865934085194845*root(s5^4 + 320*s5^3*pi + (3865934085194845*s5^2)/8589934592 + (6658547215630757*s5)/67108864 + 817126680320561/32768, s5, k) + 4123168604160*root(s5^4 + 320*s5^3*pi + (3865934085194845*s5^2)/8589934592 + (6658547215630757*s5)/67108864 + 817126680320561/32768, s5, k)^2*pi + 17179869184*root(s5^4 + 320*s5^3*pi + (3865934085194845*s5^2)/8589934592 + (6658547215630757*s5)/67108864 + 817126680320561/32768, s5, k)^3 + 426147021800368448)), k, 1, 4) + (10793252814848*symsum(exp(root(s5^4 + 320*s5^3*pi + (3865934085194845*s5^2)/8589934592 + (6658547215630757*s5)/67108864 + 817126680320561/32768, s5, k)*t)/(2*(17179869184*root(s5^4 + 320*s5^3*pi + (3865934085194845*s5^2)/8589934592 + (6658547215630\\\r\n757*s5)/67108864 + 817126680320561/32768, s5, k)^3 + 4123168604160*pi*root(s5^4 + 320*s5^3*pi + (3865934085194845*s5^2)/8589934592 + (6658547215630757*s5)/67108864 + 817126680320561/32768, s5, k)^2 + 3865934085194845*root(s5^4 + 320*s5^3*pi + (3865934085194845*s5^2)/8589934592 + (6658547215630757*s5)/67108864 + 817126680320561/32768, s5, k) + 426147021800368448)), k, 1, 4))/5

错误使用 ezplot (line 162)
不能在 xy 平面中对 8589934592*symsum((root(s5^4 + 320*s5^3*pi + (3865934085194845*s5^2)/8589934592 +
(6658547215630757*s5)/67108864 + 817126680320561/32768, s5, k)*exp(root(s5^4 + 320*s5^3*pi +
(3865934085194845*s5^2)/8589934592 + (6658547215630757*s5)/67108864 + 817126680320561/32768, s5,
k)*t))/(2*(3865934085194845*root(s5^4 + 320*s5^3*pi + (3865934085194845*s5^2)/8589934592 +
(6658547215630757*s5)/67108864 + 817126680320561/32768, s5, k) + 4123168604160*root(s5^4 + 320*s5^3*pi +
(3865934085194845*s5^2)/8589934592 + (6658547215630757*s5)/67108864 + 817126680320561/32768, s5, k)^2*pi +
17179869184*root(s5^4 + 320*s5^3*pi + (3865934085194845*s5^2)/8589934592 + (6658547215630757*s5)/67108864 +
817126680320561/32768, s5, k)^3 + 426147021800368448)), k, 1, 4) + (10793252814848*symsum(exp(root(s5^4 + 320*s5^3*pi
+ (3865934085194845*s5^2)/8589934592 + (6658547215630757*s5)/67108864 + 817126680320561/32768, s5,
k)*t)/(2*(17179869184*root(s5^4 + 320*s5^3*pi + (3865934085194845*s5^2)/8589934592 + (6658547215630757*s5)/67108864 +
817126680320561/32768, s5, k)^3 + 4123168604160*pi*root(s5^4 + 320*s5^3*pi + (3865934085194845*s5^2)/8589934592 +
(6658547215630757*s5)/67108864 + 817126680320561/32768, s5, k)^2 + 3865934085194845*root(s5^4 + 320*s5^3*pi +
(3865934085194845*s5^2)/8589934592 + (6658547215630757*s5)/67108864 + 817126680320561/32768, s5, k) +
426147021800368448)), k, 1, 4))/5 绘图。

出错 sym/ezplot (line 61)
   h = ezplot(fhandle(f)); %#ok<EZPLT>

出错 ilaplace (line 9)
ezplot(f)

请问大神们应该怎么修改才能画出f的图像呢 谢谢！

daydayout

H=simplify(vpa(f))
ezplot(H)

替代原来的ezplot(f)