function VAchar
I=[0:10:80];
U=[0,1,3:3:12,14,15,15];
plot(I,U,'ko');
grid
k3=polyfit(I,U,2);
disp('koeficienty')
disp(k3)
xj=[-0.5 : 0.1 : 80.5];
hp3=polyval(k3,xj);
hold on
plot(xj,hp3,'r--')
legend('body','pol. 2.st','Location','SouthEast')
hold off
end

