[x,y] = meshgrid(-2:0.2:2, -2:0.2:2);

u = ones(size(x));

v = exp(-y).*cos(x);

figure


quiver(x,y,u,v)
title('$\frac{dy}{dx} = e^{-y} \cos(x)$', 'interpreter', 'latex');

hold on

t = [-2:0.01:2];
plot(t, log(sin(t)+exp(1)), 'Linewidth', 3)

plot(t, log(sin(t)+1.1), 'Linewidth', 3)
