Supplement 2. The Matlab source code containing the optimization algorithm of pollen morphology used in this paper, and a short description of the program.

<h2>File List</h2><blockquote> <p><a href="pollen_optimization.m">pollen_optimization.m</a> -- (MD5: d24091e4790ac8d987991bf27ee743fb)</p> <p><a href="num_int_ell.m">num_int_ell.m</a> -- (MD5: 60c749bee5938bb138447ecabce826ed)</p> <p><a href="newnumint2.m">newnumint2.m</a> -- (MD5: c77139938b0734a8cbb0e8e59e1ba2df)</p> <p><a href="grule2dgen.m">grule2dgen.m</a> -- (MD5: 5beca3e5899d4fd3db3eced617ea4ba9)</p> <p><a href="grule.m">grule.m</a> -- (MD5: 690023ef7e1549a8f2820cc056ca5e9c)</p> <p><a href="gquad2dgen.m">gquad2dgen.m</a> -- (MD5: 1e951611314354c760d7c1130e677a4b)</p> <p><a href="gquad.m">gquad.m</a> -- (MD5: b9e62facef85dedf2041c5854c938a88)</p> <p><a href="dehydration.m">dehydration.m</a> -- (MD5: 433c4f74bd906abdc633521cbb2f5c2e)</p> </blockquote><h2>Description</h2><blockquote> <p>The code included in pollen_optimization.m allows one to calculate fitness landscape on the polar-equatorial axis plane according to formula 13 given in the model description. To run the program place file pollen_optimization.m with functions: num_int_ell.m, newnumint2.m, grule2dgen.m, grule.m, gquad2dgen.m, gquad.m, dehydration.m in one folder and run pollen_optimization.m under MATLAB with Statistic and Symbolic Math Toolbox. All variables used in the program were described in pollen_optimization.m.</p> </blockquote>