function [T,v] = ffinitge(FunFcn, FunF1, n)

%%ffinitge calcula la aproximación mediante diferencias finitas
%% de y''-c(x)*y=d(x) y(0)=0,y(1)=0, en los puntos del intervalo x_i=i*h, h=1/(n+1), i=1,2,,,,n
%%lee las funciones c(x) y d(x), y  el numero n. 
%%%c y d  estan en los ficheros  fmif y fmifi (modificaciones de fmif.m)
%%%Nos devuelve la solución aproximada en los puntos x_i de la particion de [0,1] y la grafica


h=1/(n+1);
x=0+h:h:1-h;
%%defino la matriz A  del sistema
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