```
In [3]:
```from numpy.random import normal
from scipy.stats import norm

We have the following constants for the two lizard species.

```
In [10]:
```sexual_mean, sexual_standard_deviation = 1.1, 0.15
asexual_mean, asexual_standard_deviation = 1.2, 0.3

```
In [25]:
```pod_sexual = norm.cdf(0, loc=sexual_mean, scale=sexual_standard_deviation)
pod_asexual = norm.cdf(0, loc=asexual_mean, scale=asexual_standard_deviation)
print("The probability that the sexual lizards die off in a given generation is {}".format(pod_sexual))
print("The probability that the asexual lizards die off in a given generation is {}".format(pod_asexual))

```
```

```
In [26]:
```expected_generations_sexual = int(1/pod_sexual)
expected_generations_asexual = int(1/pod_asexual)
print("The sexual lizards are expected to survive {} generations".format(expected_generations_sexual))
print("The asexual lizards are expected to survive {} generations".format(expected_generations_asexual))

```
```