syms x;p=((31.5*x)^97/factorial(97))/sum(((31.5*x).^(0:97))./factorial(0:97));x=1:0.1:20;plot(x,subs(p),'ro');
clear;clcsyms xp=((31.5*x)^97/97)/sum(((31.5*x).^(0:97))./factorial(0:97));x=01:0.1:5;plot(x,subs(x),'*')