I am just getting started with Mathematica, and I have what should be a fairly simple question about substitutions, but I cannot get it to work.
I would like to find the Euler-Lagrange equations for the functional of the function phi [x, y], and then make a replacement for the function phi [x, y]
If I enter the following:
VariationalD[tau*phi[x, y]^2 - 2*phi[x, y]^4 + phi[x, y]^6 + Dot[D[phi[x, y], {{x, y}}], D[phi[x, y, {{x, y}}]]], phi[x, y], {x, y}]
I get
Plus[Times[2,tau,phi[x,y]],Times[-8,Power[phi[x,y],3]],Times[6,Power[phi[x,y],5]],Times[-2,Plus[Derivative[0,2][phi][x,y],Derivative[2,0][phi][x,y]]]]
Now, if I try % /. phi[x,y] -> phi0[x,y] + psi[x,y], it makes a replacement for all polynomial members, but not for derived members.
How to force a substitution in these functions?