x = newton('vattenfunk',[1 1 1 1 1]'*1e-10,1e-16);

%HCO3 =[10e-3 1.03e-3 0.35e-3 0.145e-3 0.087e-3 0.0645e-3 0.057e-3 0.055e-3 0.055e-3 0.054e-3 0.054e-3 0.054e-3];
%CO3 =[5.96e-6 5.96e-6 5.96e-6 5.96e-6 5.96e-6 5.96e-6 5.96e-6 5.96e-6 5.959e-6 5.96e-6 5.96e-6 5.95e-6];
OH = [ 2 3 3.5 4 4.5 5 5.5 6 6.5 7 7.5 8];
pH = [8.9002 7.8946 7.3832 6.8853 6.5602 6.4131 6.3586 6.3414 6.3361 6.3344 6.3339 6.3337];
xlabel('-log_1_0[OH^-]')
ylabel('pH')
title('Buffertverkan')
grid
