% OCW Universidad de Cantabria
% Cálculo Simbólico y Numérico en Ecuaciones Diferenciales
% © Cursos 2001-2014,  M. Eugenia Pérez, UC
%
% y''-c(x)*y=d(x) y(0)=0,y(1)=0  
function [T,v] = ffinitge(FunF1, FunFcn, n)
h=1/(n+1);
x=0+h:h:1-h;
a=2*ones(n,1);
b=-ones((n-1),1);
Ann=diag(a)+diag(b,1)+diag(b,-1);
c=feval(FunF1,x);
An=Ann+h^2*diag(c);
g=feval(FunFcn, x);
f=-h^2*g(:);
u1=An\f;
t=[0,x,1]; 
T=t(:);
v=[0;u1;0];
plot(T,v)
end