Help me outperform Matlab in numerical solution of semidiscretized PDE system as much as possible

[ChrisRackauckas]

@RobertS

@RobertS

• I tried to modellingtoolkitize my DAE function, but that failed with the following error: LoadError: ArgumentError: invalid index: α₈ of type Num . What did I do wrong? And by the way, what benefit can I expect when using the modellingtoolkitize d system?

This would be the same issue as the sparsity one.