function [I_t,S_t] = run_SIS_hypergraph_disease(beta,delta,I_t,S_t,I,T,dt,K,fun)

n_steps = floor(T/dt);
for i=2:n_steps+1
    [I_t(:,i),S_t(:,i)] = SIS_hypergraph_disease_iterate(beta,delta,I_t(:,i-1),S_t(:,i-1),I,dt,K,fun);
end

end