It is shown that we can design numerical methods and create software that would be universally applicable for all kinds of interfaces and for any number of layers. In addition, accurate results are presented for various layer configurations and interfaces, which can serve as benchmarks for developers of new methods.