Richard's Curve and Derivative Probability Plot

These routines create probability plots of the Richard's Curve and its derivative when given the minimum and maximum values of its 95% confidence interval parameters.
It assumes that the growth is contained within a calendar year and that the units of time are given in “decimal years”.

The code can be amended if the time units are different (see code for line number) but be aware that the scaling of the x-axis and the tick marks labels for the resulting figure still assume that the calculations were done in decimal years so those will need to be changed as well.

The first part of the routine introduces Richard's Curve, defines its parameters, prompts the user to input their values, and calculates a family of curves using the user input.

The second part of the routine will generate a matrix with a frequency analysis of of the matrix generated in the first part. It counts how many times does the value of the regression falls between two values at a given time of the year.

The third part of the routine plots the results. It uses the "Custom Colormap" routine by Martinez-Cagigal (2020) to create a colourblind-friendly colour map for the figure

References:
Víctor Martínez-Cagigal (2020). Custom Colormap (https://www.mathworks.com/matlabcentral/fileexchange/69470-custom-colormap), MATLAB Central File Exchange. Retrieved June 19, 2020.