% Attention, utilise annee_version
f=@(t) cos(2*t).*t.^4;
g=@(x) sin(x)./x;
X0=0.5;

syms t x;
F=int(f(x),x,0,t);
disp('vérification');
disp(diff(F));
Gp=diff(g(x));
% Vieille syntaxe utilisant inline à deconseiller !!!!
% Ff=inline(char(vectorize(F)));
% Gf=inline(char(vectorize(Gp)));
% Vieille syntaxe utilisant @
if annee_version>=2019
    Ff=matlabFunction(F);
    Gf=matlabFunction(Gp);
else
    eval(['Ff=@(t)',vectorize(char(F)),';']);
    eval(['Gf=@(x)',vectorize(char(Gp)),';']);
end
% pas optimal mais developpement_limite(0,1,char(g(x))) marche pas !
% Attention, adaptée éventuellement à la version 2019 de matlab !
if annee_version>=2019
    g1=taylor(g(x),'ExpansionPoint',X0,'Order',2);
    g2=taylor(g(x),'ExpansionPoint',X0,'Order',3);
else
    g1=taylor(g(x),X0,2);
    g2=taylor(g(x),X0,3);
end
if annee_version>=2019
   g1=matlabFunction(g1); 
   g2=matlabFunction(g2); 
else
    eval(['g1=@(x)',vectorize(char(g1)),';']);
    eval(['g2=@(x)',vectorize(char(g2)),';']);
end


X=linspace(-pi,pi,1000);
Y1=feval(g1,X);
Y2=feval(g2,X);
subplot(2,2,1);
plot(X,f(X));
title('f');
subplot(2,2,2);
plot(X,Ff(X));
title('F');
subplot(2,2,3);
plot(X,g(X),X,Y1,X,Y2,X0,g(X0),'o');
legend('G','ordre 1','ordre 2','location','best');
title('G et ses deux DL');
subplot(2,2,4);
plot(X,Gf(X));
title('G''');