Knowing the distribution of stellar rotational velocities is essential for the understanding stellar evolution. Because we measure the projected rotational speed vsini, we need to solve an ill-posed problem given by a Fredholm integral of the first kind to recover the true rotational velocity distribution. After discretization of the Fredholm integral, we apply the Tikhonov regularization method to obtain directly the probability distribution function for stellar rotational velocities. We propose a simple and straightforward procedure to determine the Tikhonov parameter. We applied Monte Carlo simulations to prove that Tikhonov method is a consistent estimator and asymptotically unbiased. This method is applied to a sample of cluster stars. We obtain confidences intervals using bootstrap method. Our results are in good agreement with the one obtained using the Lucy method, in recovering the probability density distribution of rotational velocities. Furthermore, Lucy estimation lies inside our confidence interval. Tikhonov regularization is a very robust method that deconvolve the rotational velocity probability density function from a sample of vsini data straightforward without needing any convergence criteria.

To better understand the evolution of stars in binary systems as well as to constrain the formation of binary stars, it is important to know the binary mass-ratio distribution. However, in most cases, i.e. for single-lined spectroscopic binaries, the mass ratio cannot be measured directly but only derived as the convolution of a function that depends on the mass ratio and the unknown inclination angle of the orbit on the plane of the sky. We extend our previous method to deconvolve this inverse problem (Cure et al. 2014), i.e., we obtain as an integral the cumulative distribution function (CDF) for the mass ratio distribution. After a suitable transformation of variables it turns out that this problem is the same as the one for rotational velocities $v \sin i$, allowing a close analytic formulation for the CDF. We then apply our method to two real datasets: a sample of Am stars binary systems, and a sample of massive spectroscopic binaries in the Cyg OB2 Association.} {We are able to reproduce the previous results of Boffin (2010) for the sample of Am stars, while we show that the mass ratio distribution of massive stars shows an excess of small mass ratio systems, contrarily to what was claimed by Kobulnicky et al. (2014). Our method proves very robust and deconvolves the distribution from a sample in just a single step.

Rotational speed is an important physical parameter of stars and knowing the distribution of stellar rotational velocities is essential for the understanding stellar evolution. However, it cannot be measured directly but the convolution of the rotational speed and the sine of the inclination angle, $v \sin i$. We developed a method to deconvolve this inverse problem and obtain the cumulative distribution function (CDF) for stellar rotational velocities extending the work of Chandrasekhar & M\"unch (1950). This method is applied a) to theoretical synthetic data recovering the original velocity distribution with very small error; b) to a sample of about 12.000 field main--sequence stars, corroborating that the velocity distribution function is non--Maxwellian, but is better described by distributions based on the concept of maximum entropy, such as Tsallis or Kaniadakis distribution functions. This is a very robust and novel method that deconvolve the rotational velocity cumulative distribution function from a sample of $v \sin i$ data in just one single step without needing any convergence criteria.

We calculated the influence of the limb-darkened finite disk correction factor in the theory of radiation-driven winds from massive stars. We solved the 1-D m-CAK hydrodynamical equation of rotating radiation-driven winds for all three known solutions, i.e., fast, \Omega-slow and \delta-slow. We found that for the fast solution, the mass loss rate is increased by a factor \sim 10%, while the terminal velocity is reduced about 10%, when compared with the solution using a finite disk correction factor from a uniformly bright star. For the other two slow solutions the changes are almost negligible. Although, we found that the limb darkening has no effects on the wind momentum luminosity relationship, it would affect the calculation of synthetic line profiles and the derivation of accurate wind parameters.

We present a general method for solving the non--linear differential equation of monotonically increasing steady--state radiation driven winds. We graphically identify all the singular points before transforming the momentum equation to a system of differential equations with all the gradients explicitly give. This permits a topological classification of all singular points and to calculate the maximum and minimum mass--loss of the wind. We use our method to analyse for the first time the topology of the non--rotating frozen in ionisation m--CAK wind, with the inclusion of the finite disk correction factor and find up to 4 singular points, three of the x--type and one attractor--type. The only singular point (and solution passing through) that satisfies the boundary condition at the stellar surface is the standard m--CAK singular point.

The effects of rapid rotation and bi--stability upon the density contrast between the equatorial and polar directions of a B[e] supergiant are investigated. Based on a new slow solution for different high rotational radiation--driven winds and the fact that bi--stability allows a change in the line--force parameters ($\alpha$, $k$, and $\delta$), the equatorial densities are about $10^2$--$10^3$ times higher than the polar ones. These values are in qualitative agreement with the observations. This calculation also permits to obtain the aperture angle of the disk.

The effects of rapid rotation and bi-stability upon the density contrast between the equatorial and polar directions of a B[e] supergiant are re-investigated. Based upon a new slow solution for different high rotational radiation driven winds (Cur\'e 2004) and the fact that bi--stability allows a change in the line--force parameters ($\alpha$, $k$, and $\delta$), the equatorial densities are about $10^2$--$10^4$ times higher than the polar ones. These values are in qualitative agreement with the observations.

The topological analysis from Bjorkman (1995) for the standard model that describes the winds from hot stars by Castor, Abbott & Klein (1975) has been extended to include the effect of stellar rotation and changes in the ionization of the wind. The differential equation for the momentum of the wind is non--linear and transcendental for the velocity gradient. Due to this non--linearity the number of solutions that this equation possess is not known. After a change of variables and the introduction of a new physically meaningless independent variable, we manage to replace the non--linear momentum differential equation by a system of differential equations where all the derivatives are {\it{explicitely}} given. We then use this system of equations to study the topology of the rotating--CAK model. For the particular case when the wind is frozen in ionization ($\delta=0$) only one physical solution is found, the standard CAK solution, with a X--type singular point. For the more general case ($\delta \neq 0$), besides the standard CAK singular point, we find a second singular point which is focal--type (or attractor). We find also, that the wind does not adopt the maximal mass--loss rate but almost the minimal.