% Matematica Aplicada 2008-2009 Carlos Balsa
% exercicio 1 (1) da ficha nº 7
% resolução pelo metodo das tentativas

clear all;
clf;

ti=0;  tf=1;              % intervalo de tempo         
ui=0;  uf=1;              % condições de fronteira

%h=0.2;                    % passo de integração
%np=(tf-ti)/h+1;           % numero de pontos

pvi=@(t,z) [z(2) 10*z(1)^3+3*z(1)+t]';    % u''=F(t,u,u´)
[t, u]=ode_shoot(pvi,ti,tf,ui,uf);        % metodo das tentativas

% grafico 
plot(t,u(:,1),'r*-'); hold on;
plot(t,u(:,2),'b+-'); hold on;
grid on; 
xlabel('t'); 
ylabel('u');
legend('u_1','u_2');


