It's a shame nobody ever answered this question years ago when asked. As an exercise for myself I will attempt it.
There's a useful answer with some python here: https://stackoverflow.com/a/43213692/1335793 That shows some varied parameterisations to shift or scale the curve.
In your examples, these two are equivalent:
$w = w.max / 1 + e ^-k (t - tm)$
$y = C / 1 + Ae ^-Bx$
In the first one, the sigmoid is applied to growth over time, so the x axis is assumed to be time. t-tm is the difference between the start time and some point later in time, which is really the same as an x coordinate. So t-tm is the same as x in the second formula.
What IS (B) actually representing here?
B is positive or negative growth rate. If you make it negative, the curve will be a mirrored, i.e. it will start high and get lower over time. If it is positive the curve will start low and grow higher over time. In the paper you linked http://www.cs.xu.edu/math/math120/01f/logistic.pdf the first page shows a positive growth curve, the second page a negative curve. There is a relationship between the growth rate and the inflection point. If you modify the value of the growth rate, it changes the steepness of the slope in the middle of the curve. The paper explains that if you know the inflection point you can then back-calculate the growth rate.
And how, why, and where did they assign it the value of 0.0266 ???
This is a bit frustrating to learn about, as it seems arbitrary and in fact it is to some extent. They started with the data points, and then recognized that they kind-of followed a logistic distribution, and then attempted to "fit" the right function to the data by modifying the parameters of the logistic function. This is usually done with computational methods like ordinary-least-squares or gradient-descent, which are akin to starting with a random guess, then changing it slightly until the curve is as close to each of the data points as is possible within the constraints of the chosen function. It could just as well be done by trial-and-error, mucking around with the parameters until it looks right. If that doesn't feel very elegant then you would be right, it is a practical brute-force approach to applied mathematics.
And lastly, how can we determine what parameters to use for a sigmoid
function?
The convention for functions that describe curves in two dimensions is to use some form of
$y = f(x)$
where the $y$ coordinate is calculated by applying some function to the $x$ coordinate. This assumes you know what all the $x$ coordinates are, but that is not a problem for idealized functions like the logistic (sigmoid distribution). In the growth-of-population example, you know that time is continuous and the unknown is "what will the population be at any point in time?". Using the known (time) value you can then try to predict the unknown (population) value.
Choosing the syntax to use doesn't really matter that much, although it does help to use conventions, or for specific examples some variables that have meaning in the context. In the generic formula
$y = C / 1 + Ae ^-Bx$
$y$ will be the population value, $C$ will be the highest possible population value, and $e$ is Euler's number.
The generic logistic model is symmetrical around 0 but can be shifted left-right (x axis offset), up-down (y-axis offset), scaled (amplitude) or "stretched", changing the slope of the growth. You need to add additional parameters to the model to control those aspects, which makes the model describe your observed data. The approach of starting with a generic model and customizing it to data needs to balance goodness of fit with complexity (not too many parameters). The process of model selection within the field of statistical modelling is dedicated to doing just that.
Here's a simple example: $f(x,y)=xy$. You can change the powers of $x$ and $y$ to 'bias' one of $x$ or $y$.
Here's another function that is "sigmoid-shaped" and is close to your required values, as per your comments:
$$f(x,y)=\left( \frac{1}{1+e^{-100x+50}}\right)\left( \frac{1}{1+e^{-100y+50}}\right)$$
This is approximately $0$ at $(0,0)$ and $1$ and $(1,1)$, and is sigmoid-shaped in the sense that when either $x=1$ or $y=1$ the resulting slice of the 3D-plot is a sigmoid, as seen here:
Theoretically, you could also make this exactly equal to $1$ at $(1,1)$ and $0$ at $(0,0)$ by multiplying by a constant and changing the values of $50$ in the exponents appropriately to ensure that the function went through the points $(0,0,0)$ and $(1,1,1)$.
Best Answer
We can start with a simple soft staircase function:
$$ f(x) = x - sin \space x $$
and then feed it into itself:
$$ y(x) = f(f(x)) $$
then again:
$$ y(x) = f(f(f(x))) $$
and again:
$$ y(x) = f^4(x) $$
As you can see, each iteration makes the "flat" part of the step longer, and the rise steeper.
The period and the height of each step is $ 2 \pi $, so multiply $x$ by $2 \pi / w$ and $y$ by $h / 2 \pi$ to reach your desired scale.
In reality, the curve is only truly flat (zero derivative) at the centre of each step — at every $ 2 \pi k $ — and only close to flat on either side of that point.
Configurability is limited: The softness of the step can only be specified in integer amounts (the number of times we reapply $f$ to itself), and it requires many/infinite applications to make the step really sharp.