As it can be seen in the image, the upper bound of the standard deviation is considering "time" as the number of datapoints (200) instead of the actual time ([0, 100] with intervals of 0.5).
Replication
Replicated using the Michaelis Menten example with a timespan of TimeSpan.linspace(t=100, num_points=201)
Proposed Fix
Could be solved by changing line 821 of gillespy2/core/results.py to
plt.plot(average_result['time'], lowerBound, average_result['time'], upperBound, color='grey', linestyle='dashed')
Error
As it can be seen in the image, the upper bound of the standard deviation is considering "time" as the number of datapoints (200) instead of the actual time ([0, 100] with intervals of 0.5).
Replication
Replicated using the Michaelis Menten example with a timespan of
TimeSpan.linspace(t=100, num_points=201)
Proposed Fix
Could be solved by changing line 821 of gillespy2/core/results.py to plt.plot(average_result['time'], lowerBound, average_result['time'], upperBound, color='grey', linestyle='dashed')