And it's definitely used by real engineers. It's for making switching power supplies, so it can handle fairly complex circuits. I wouldn't venture into RF/high speed circuits with it though, since the included parasitics aren't sufficient.
And, custom functions are supported, so, I believe, you can put whatever maths in a component (I've never ventured this far): https://ltwiki.org/index.php?title=B_sources_(complete_refer...
Most EDA packages will include some sort of simulator, often SPICE-based. The fancier ones for RF work include or integrate with field solvers.
https://computationstructures.org/exercises/sandboxes/jade.h...
As far as the linked article "Fifty Nifty Variations of Two-Transistor Circuits: A tribute to the versatility of MOSFETs", I will mention that CircuitLab does NOT currently provide a MOSFET model that includes a separate body terminal.
Instead, the body (also sometimes called the "back gate") and source terminals are assumed to always be internally connected in the CircuitLab MOSFET models, resulting in a three-terminal device.
Almost all real-world discrete MOSFETs you can buy are also three-terminal (gate/drain/source), not four (gate/drain/source/body). Some discussion here https://electronics.stackexchange.com/questions/137161/why-a... and here https://electronics.stackexchange.com/questions/185109/mosfe... but not especially satisfying in my opinion.
As a result, Figures 2, 16, 27, 29, 46, 49, and 50 can NOT currently be realized within CircuitLab.