classic paper. that was the first paper i was asked to read as a part of my project. the amount of work contained in it is enormous! they didn't have computers to help them with numerics at the time. simply phenomenal!
ah.. that's an interesting question. i think since the finite difference scheme is stable, there would be an error but it will still largely satisfy the original pdes, the stability condition being D*ht/hx/hx being less that 0.25. this can vary with the boundary conditions but for this set of...
this is in continuation to the previous post. i attached the plots i obtained as well as the laplacian for the code i used. its taken from leppanen's code of course. i seem to have to gave obtained the reverse of what you did. i got stripes where you got spots ans spots where you got stripes.
hello jackmell. i wasnt able to log on till now as i had my semester end exams for the last couple of weeks. i would definitely like to know how the code in c works. i have never used c. the languages i know are fortran, basic and matlab. however, i followed up your links and obtained the code...
ahh.. that's too bad.. i was hoping to read some work at a masters level on this topic.. i don't know anyone in my university who works on this infact.. hm.. i am not sure where the problem is occurring but i have attached a code to this thread. the first post. i wrote it for the explicit...
hello filip. thank you so much for your response. in fact i have read this particular paper once before. though i got quite lost halfway through the paper! its quite mathematically involved! i haven't even reached a stage where i can analyse equations in such a comprehensive manner. i am still...
there was a typo in the equations i sent you. the second equation has the term Bu not bu. there is no b in the equations.
the equations are
du/dt = Du*(d2u/dx2 + d2u/dy2) + A - (B + 1)*u + v*u^2
dv/dt = Dv*(d2v/dx2 + d2v/dy2) + B*u - v*u^2
A = 3,
B = 9,
Du = 5;
Dv = 12;
The...
thank you so much jackmell for your response. but there's a msicommunication here. the code is not a mathematica code. its a MATLAB code. and are you sure i can generate turing patterns using the initial conditions like the one you have set? correct me if i am wrong but i thought the initial...
so my project is on reaction diffusion equations in 2d. i have been asked to reproduce known and published work to start off with and the system is the brusselator equations with diffusion terms. to put it simply my program doesn't work. now the parameter values for the equations are taken from...