elfin

Esta es una vista previa del archivo. Inicie sesión para ver el archivo original

function solucion=elfin(c,f,a,b,n)
 h=(b-a)/n;
 %Vector de los x gorro
 xx=zeros(n,1);
 for i=1:n
 xx(i)=a+(i-0.5)*h; 
 end
 %Matriz de sistema
 M=diag((2/h+2*c*h/3)*ones(n-1,1))+diag((-1/h+c*h/6)*ones(n-2,1),1)+diag((-1/h+c*h/6)*ones(n-2,1),-1);
 
 %Vector del lado derecho
 F=zeros(n-1,1);
 for i=1:(n-1)
 F(i)=0.5*h*(feval(f,xx(i))+feval(f,xx(i+1)));
 end
 solucion=[0;M\F;0];
end