function [A_ss,B_ss,C_ss,D_ss]=crea_DLTI_StateSpace_calore_1D(m,Dx,c)

n=m-1;

A=sparse(n+1,n+1);

B_ss=zeros(size(A,1),1);

% 1(Nf-N0): dT/dx(t,0)=f0(t), forzante nota, f0 inteso come flusso di calore

%           dT/dx (t,L)=0, condizioni di Neumnann omogenee



for i=2:n % nodi interni: discretizzazione del laplaciano

    A(i,i) = -2/Dx^2;

    if i > 1      A(i,i-1) = 1/Dx^2;   end

    if i < n+1    A(i,i+1) = 1/Dx^2;   end

end







A(1,1)=-2/Dx^2; A(1,2)=2/Dx^2;

A(n+1,n+1)=-2/Dx^2; A(n+1,n)=2/Dx^2;

forzante=0;

save forzante forzante;

B_ss(1)=2/Dx;



% creo le matrici del sistema state-space continuo:

A_ss=A;

C_ss=zeros(1,size(A,1));

C_ss(1)=1;

D_ss=[0];



