@tomleslie 

If it is posible, please help me  more...

Thank you very much.

Good luck

Thank you tomleslie 1674

I tried to correct what was wrong in this mw file(bellow).

pdeProblem.R.mw

1). Undefined parameter is assumed to be codtant. like bellow:(GAMMA(omega[0]) is a function which have a difined value)

omega[0] := 12;
C := 3*10^8;
n[2] := 1.2;
n[1] := 1.56;
n[c] := 1.55;
a := 2*10^(-6);
alpha := .2;
beta[1] := 1;
beta[2] := .2;
beta[3] := 2.2;
k[0] := omega[0]/C;
T[R] := .2;
v := k[0]*a*(n[1]^2-n[c]^2)^(1/2);
w := (.65+1.619/v^(3/2)+2.879/v^6)*a;
F(x,y):=exp(-(x^(2)+y^(2))/(w^(2))):
GAMMA(omega[0]) := evalf(omega[0]*n[2]*(int(int(abs(F(x, y))^4, x = -infinity .. infinity), y = -infinity .. infinity))/(C*(int(int(abs(F(x, y))^2, x = -infinity .. infinity), y = -infinity .. infinity))));

2. boundary conditions should modified to:

3. A(z,t) is an amplitude of EM wave.  it is not complex but the PDE equation is complex.