function disegna(ndom,imodo,gammas,coef,freprop,LL)
x=[];
y=[];

for k=1:ndom
    c=coef(4*(k-1)+1:4*k,1)';
    xx=[0:0.01:1].*LL(k);
    yy=c(1)*cos(gammas*xx)+c(2)*sin(gammas*xx)+c(3)*cosh(gammas*xx)+c(4)*sinh(gammas*xx);
    if k>1
        xx=xx+LL(k-1);
    end
    x=[x xx];
    y=[y yy];
end
figure(imodo+1)
plot(x,y);grid;title(['Modo n.: ' int2str(imodo) ' frequenza: ' num2str(freprop) ' Hz '])