simulate function code:
function [V] = simulate(I , params) %parameters alphaa = params(1); k = params(2); ka = params(3); delb = params(4); x4 = params(5); x3 =