[jakubstary]

Hi,

I need to solve the first order macroscopic pedestrian flow model in 2D in FreeFem:

rho_t + div(rho*V(rho)*mu) = 0,

where V(rho) is some speed-density relation and mu is constant for my case, but I can't figure out, how to write the code of my problem in FreeFem.

I think that it will be:

problem pedestrian(rho, rhoh) =
int2d(Th)( rho*rhoh/dt )
+int2d(Th)( ....... )
+on (some boundary condition);

Please, can you help me with it? I'm a beginner in FreeFem.
Thanks.