lab 2 / vpasolve9

syms x
r=vpasolve(2*x.^3+2==8*sin(x),x)
r2=vpasolve(2*x.^3+2==8*sin(x),x,pi/2)
x=-4:0.1:4;
y=8*sin(x)
p=2*x.^3+2;
figure;
plot(x,y,'k-.',x,p,'r:');
grid on;
title('Two graphs on one axis')
xlabel('X_axis')
ylabel('Y_axis')
legend({'y(x)','p(x)'},'Location','best')

% SECOND FUNCTION USING VPASOLVE exp(-cos(x))+sin(x)-0.5*x.^2
syms x
r3=vpasolve(0.5*x.^2-sin(x)==exp(-cos(x)),x)
r4=vpasolve(0.5*x.^2-sin(x)==exp(-cos(x)),x,pi/2)
x=-4:0.1:4;
y2=exp(-cos(x));
p2=0.5*x.^2-sin(x);
figure;
plot(x,y2,'k-.',x,p2,'r:');
grid on;
title('Two graphs on one axis')
xlabel('X_axis')
ylabel('Y_axis')
legend({'y2(x)','p2(x)'},'Location','best')