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

clear all;
clf;

ti=0;  tf=1;               % intervalo de tempo         
ui=0;  uf=1;               % condições de fronteira
h=0.1;                     % espaçamento da malha
np=(tf-ti)/h+1;            % numero de pontos totais
n=np-2;                    % numero de pontos interiores

rhs=@(t,y) [10*y.^3+3*y+t.^2];              % u''=F(t,u,u´)
[t, u]=ode_finit_diff(rhs,ti,tf,ui,uf,n);   % met. dif. finitas

% grafico 
plot(t,u,'r*-'); hold on;
grid on;
xlabel('t'); ylabel('u');
title('Método das Diferenças Finitas');


