Add numerical experiment code for the paper #502

Job	Run time
Julia 1 - windows-latest - x64 - pull_request	10m 37s
Julia 1 - ubuntu-latest - x64 - pull_request	9m 50s
Julia 1.10 - ubuntu-latest - x86 - pull_request	8m 30s
Julia 1.10 - ubuntu-latest - x64 - pull_request	7m 39s
	36m 36s

Provide feedback