Tspan = linspace(0,300e-7,1000);

c_H=linspace(0.8,1.1,20);
k=1;

for i = c
   disp(i)
   [t u] = ode15s('diff', Tspan, [1 10 0 0 c_H]);
   pOH(k) = -log10(u(1000,4));
   k=k+1;
end

plot(c,14-pOH)  %ekvivalenspunkten c) iv)

%figure(2) %b)
%plot(t,u(:,1),t,u(:,3:5))
%title('[NH_3] = 1 M  [H^+] = 0.5 M ( c) i) )')
%title('Ammoniaks protolys c)')
%ylabel('Koncentration (M)')
%xlabel('Tid (s)')
%grid
%legend('NH_3','NH_4^+, H_2O','OH^-','H^+')

%figure(2) %b) i)
%plot(t,-log10(u(:,4)),t,14+log10(u(:,4)))
%title('pH och pOH för 0.1 M NH_3') 
%axis([0 1e-7 0 14])
%grid
%xlabel('Tid (s)')
%legend('pOH','pH')

%b) ii) löses genom en FOR-loop över ODE-lösaren
%A = linspace(0,0.05,30);
%pOH(i) = -log10(u(1000,4));
%plot(A,14-pOH)
%xlabel('[NH_3]  (M)')
%ylabel('pH')
%title('pH för olika lösningar av NH_3  ( b) ii) )')

%b ) iii)
%plot(t,u(:,3)./u(:,1)*100)
%title('Dissociationsgraden av 0.1M Ammoniak')
%ylabel('Procent (%)')
%xlabel('Tid (s)')
%grid

%b) iv)
%plot(t,u(:,3).*u(:,4)./u(:,1))
%title('K = 1.76*10^-^5')
%xlabel('Tid (s)')
%ylabel('[NH_4^+][OH^-]/[NH_3]')
%grid
