Oh yes I was being dim
of course you just rearrage the formula to solve for x
%-- rearrange the formula
syms param1 param2 param3 param4 x y
eqn = (y == param1+(param2-param1)./(1+10.^((param3-x)*param4)));
solx = solve(eqn, x)
%---- vreat function to solve for x
fsigm_for_x = @(param,y) param(3) - log(-(param(2) - y)/(param(1) - y))/(param(4)*log(10));
param=[0 1 5 1]; % "min", "max", "x50", "slope"
x = fsigm_for_x(param, 0.5)
