www.gusucode.com > symbolic 源码程序 matlab案例代码 > symbolic/PlotPuiseuxSeriesApproximationExample.m

    %% Plot Puiseux Series Approximation
%% 1
syms x
f = exp(x)/x;
s6 = series(f, x)
s7 = series(f, x, 'Order', 7)
s8 = series(f, x, 'Order', 8)
%% 2
fplot([s6 s7 s8 f])
legend('approximation up to O(x^6)','approximation up to O(x^7)',...
            'approximation up to O(x^8)','exp(x)/x','Location', 'Best')
title('Puiseux Series Expansion')