What I understand is that in your case u(2) is not an independent variable.
You might want to try: d [u(1)]/dt = f( d[u(1)]/dx ) d [constant*u(1)]/dt = f( d[u(1)]/dx )
Find the treasures in MATLAB Central and discover how the community can help you!
Start Hunting!