f=H/c0
plot(p(c),f)
xlabel("pc")
ylabel("Fraction ayant réagi")