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