% tema 5, test 2, 521230, S2-2011

% llamando a esquema pc de orden 2 para resolver problema en tema 5
h = 0.1;
% valor inicial dado
y0 = [-1 0];
% 2do valor inicial necesitado para esquema pc, se calcula evaluando la
% solucion exacta y su derivada en h
yh = [exp(h)/4+exp(-3*h)/12-4/3 exp(h)/4-exp(-3*h)/4];

% llamado a esquema pc
[x,y] = pc_orden2('func1_tema5',[0 4],y0,yh,h);

% graficar aproximacion a solucion y solucion exacta
plot(x,y(:,1));
hold on
plot(x,exp(x)/4+exp(-3*x)/12-4/3,'r');

% graficar aproximacion a derivada, haciendo zoom en el grafico se observa
% que ella es 1 cerca de 1.38
figure(2);
plot(x,y(:,2));
% la derivada es 1 cerca de 1.38