Hi guys, Let me start by saying that I am quite new to Julia, and this is my first serious attempt at learning this great language :) In order to do that, I have set out to solve a pair of coupled, time-dependent PDEs for a physical model that I use in my PhD thesis. The details of the problem are here: scicomp.stackexchange.com/questions/11772/numerical-approach-to-coupled-nonlinear-pdes-with-time-dependence
What I have come up with so far is this: http://nbviewer.ipython.org/github/berceanu/notebooks/blob/master/RK_OPO_julia.ipynb I am addressing the community because my idea was that, once I get this in a nice, optimal working form, I could push it to the examples repo. So the notebook I linked to is a working version, but I am quite sure that I am not really using Julia to its advantage really. I am looking forward to hear your suggestions and comments on how I can improve/optimize the code, and also make it more "julia-like" (I come from a python/numpy/scipy background). Thanks a lot in advance!
