Can you help me please how to do this kind of things using freefem++. I am facing problem regarding this.
I wan to such type of computations. I need some demo examples so that i can understand how i can create this. Please help if anyone know.
Go through the example of Laplace equation and try to replace the Laplace equation with the weak formulation of your PDE. Then you can export the solution to paraview vtk file and use paraview for generating figures like that.
I recommend checking out the examples and Videos by TuxRider for FreeFEM. They would be helpful and teach you all the necessary stuff to use FreeFEM
load "iovtk" //place this in the start of the code
savevtk("output.vtu", Th, u1,u2, p, dataname="U V P", order=[1,1,1]);
//place this after the problem is solved and you have solution.
This will output a paraview readable file which you can load in Paraview and generated images like that. Paraview is recommended for better images otherwise you can simply plot the solution and save it.
If you are not aware about Paraview. It is a post-processing tool. Tux-riders channel has a series of tutorial on that too.
The code remains same whether you use FreeFEM on windows or linux. It doesn’t change.
In windows also you can use the same command to run the code through CMD.