I have just recently started to use freefem3d to simulate mass diffusion inside a certain 3d shape. My question here is that i need to define a step function to set the initial conditions for the concentration. I would need a function that says: for x<= x1 conc(x,y,z)=3(for example), and for x>x1 conc(x,y,z)=0. And i can not figure out a way to do it.
In freefem++ is possible to define such a function by suitable if clauses inside a function statement since the return instruction gives back the result of the call to the function.
I am using the distributed binary for Macintosh.
Thank you very much!
Emilio
-- ############################################### Dr. Emilio Melero García Laboratorio de Estudios Cristalográficos, IACT Instituto de parasitología Lopez Neyra,
P.T. Ciencias de la Salud, E-18100 Armilla, Granada, España