Clc
Clc
syms x y
subplot(3,2,1)
y1= dsolve ('Dy=sin(x) -2*y/x' , 'x' );
simplify (y1)
pretty (y1)
fplot(@(x)1/x^2 - (cos(x)*(x^2 - 2) - 2*x*sin(x))/x^2,[-pi
pi],'b','linewidth',1)
xlabel ('x(secondes)','FontName','Times','Fontsize',14);
ylabel ('y1','FontName','Times','Fontsize', 14)
title ( 'Fonction y1','FontName','Times','Fontsize', 14)
subplot(3,2,2);
y2=dsolve('D3y+D2y+Dy69=x^2+sin(x)-1','x');
simplify(y2)
pretty(y2)
fplot(@(x)x^2 - (2^(1/2)*cos(x + pi/4))/2 - 2*x +exp(-
x)+1,[-2*pi 2*pi],'b','linewidth',1)
xlabel ('x(secondes)','FontName','Times','Fontsize',14);
ylabel ('y2','FontName','Times','Fontsize', 14)
title ( 'Fonction y2','FontName','Times','Fontsize', 14)
subplot(3,2,3);
y3=dsolve('Dy+2*y/(x-1)=x/ (x-1)', 'y(0)=1','x');
simplify(y3)
pretty(y3)
fplot(@(x)(2*x^3 - 3*x^2 + 6)/(6*(x - 1)^2),[-3
3],'b','linewidth',1)
xlabel ('x(secondes)','FontName','Times','Fontsize',14);
ylabel ('y3','FontName','Times','Fontsize', 14)
title ( 'Fonction y3','FontName','Times','Fontsize', 14)