One common model is that the resistance force is proportional to the speed.
Under that model, an object falling, under gravity has acceleration -g+ kv (k is the proportionality constant, v the speed. Since that is a function of v, it give the linear differential equation mdv/dt= -mg+ kv. The general solution to that is v(t)= Ce-kt/m-mg/k. For very large t, that exponential (with negative exponent) goes to 0 and the "terminal velocity" is -mg/k.
Another common model is to set the resistance force proportional to the square of the speed. That means the net force is -g+ kv2 and v satisfies the differential equation mdv/dt= -g+ kv2. That's a non-linear differential equation but is separable and first order. We can integrate it by writing
dv/(kv2-g)/m= (-1/2√(g))(1/(√(k)v+√(g))dv/m+(1/2√(g))(√(k)v-&radic(g))dv/m= dt.
Integrating both sides, we get (1/2√(kg))ln((√(k)v-√(g))/(√(k)v+√(g))= mt+ C. For large t, the denominator on the left must go to 0: the terminal velocity is -√(g/k) which, you will notice, is independent of m. This model is typically used for very light objects falling through air or objects falling through water.