% tema 4, test 2, 521230 S2, 2011

% funcion con parte derecha de PVI en func1_tema4

[t,y] = ode45('func1_tema4',[0 pi],[1 1 -1]);
plot(t,y(:,1));
hold on
plot(t,y(:,2),'r');
plot(t,y(:,3),'k');
legend('y(x)','dy/dx','d(dy/dx)/dx');

% buscando valores pedidos
[t,y] = ode45('func1_tema4',[0 pi/6 pi/3 pi],[1 1 -1]);
%y(pi/6)
y(2,1)
%y'(pi/3)
y(3,2)
%y''(pi)
y(4,3)