Then my question could simply become, how to approximate a double integral using Adaptive Simspon's Quadrature rule, which is a math question. All I really need to know is how to do only a single example, then the rest becomes easier. A pseduo-algorithm, like the ones given above, use the methods necessary to compute the math problems, but they will work for all general cases, instead of a particular one.