% preg1, tema 8, test 2, 521230 S2, 2011

% funcion con parte derecha de PVI en func1_tema8
% tamannos de paso que se utilizaran para resolver el PVI con parte derecha
% en func1_tema8
h = [0.5 0.4 0.3 0.2 0.1];
for i = 1:length(h)
    % llamando a rk
    [t,y] = rkpreg1_tema8('func1_tema8',[0 2*pi],[1 0],h(i));
    % graficando error en ambas componentes
    figure(i)
    plot(t,abs(y(:,1)-cos(t)),'.-');
    hold on
    plot(t,abs(y(:,2)-sin(t)),'r.-');
end