function vale=soluzione(x)
global u0 coef;
vale=u0;
for i=1:length(coef)
  vale=vale+coef(i)*x.^i;
end
