function [F] = y(tt,y,hx,n,c);
%% Calcula o valor da funçao F=A*Y (F e a derivada em Y)
%% tk e yk sao os vectores para os quais se quer calcular a soluçao


%% Geraçao da matriz dos coeficientes
d0=-2*ones(n,1) ;                  % vector cujos elementos sao -2
d1=ones(n-1,1) ;                   % vector de 1
A=diag(d0)+diag(d1,1)+diag(d1,-1) ;% colocaçao dos vectores sobre as diagonais
A=(c/(hx)^2)*A;

%% Calculo do valor da funçao 
F=A*y;