Freefem++

[naria]

Salut

Je besoin d'aide s'il vous plaît !!
Je veux programmer cette formulation mixte de problème d'obstacle mais je n'arrive pas.

Code : Tout sélectionner

problem pb([p1,p2],[q1,q2])=int2d(th)(p1*q1+p2*q2)
                  +int2d(th)((lamda*(div(q1,q2))+psi*(div(q1,q2))))
                 -int1d(th)(0.7*(q1*N.x+q2*N.y))+on(1,p1=0.7,p2=0.7);  \\ problème variationnelle on freefem

on a le problème est : trouver (p,$\lambda$) $p\in H(div,\Omega)$ et $\lambda$ est dans $P1$
$<p,q>+<\lambda,div(q)>=-<\psi,u-\psi>+<g,q.n> \quad \forall q \in H(div,\Omega)$
$p=gard(u) \quad \lambda =u-\psi \quad \Omega=/{(x,y) |x^2+y^2 <1\} \quad \psi=-x^2-y^2+0.3 \quad g=0.7$

Merci d'avance.