% definindo o sistema de EDOs
a=0.25;
b=0.1;
F = @(s,z) Mathieu(s,z,a,b);

% condicoes iniciais
x0 = 1;         % valor inicial de x
y0 = 0;         % valor inicial de x'

% parametros de resolucao
t0 = 0;         % tempo inicial
tf = 100;        % tempo final
h = 0.01;       % passo inicial de integracao
Tol = 1e-6;     % erro maximo tolerado

% aplicando o metodo RKF45
  
z0=[x0,y0];     % condição inicial

[t,z] = RKF45(t0,tf,h,z0,F,Tol);

% representacao grafica do resultado

figure
plot(t,z(:,1))

