%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%   use ode45 solve the Riccati Equation
%   y' = sin(x)y+cos(x)y^2+x
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

close all
clear all
  
[x, y] = ode45(@Riccati_func, [-2,  5], 2);
 
plot(x, y, 'k-',    'linewidth', 3);
box on
grid on
xlabel('x', 'fontsize', 24);
ylabel('y', 'fontsize', 24);
title('Numerical solution of Riccati Equation');