Solve doesn't seem to work

clc;
clear;
epi0=8.854e-12;
mew0=4*pi*(10^(-7));
z0=sqrt(mew0/epi0);
c0=3.8e08;
alpha=15;
b=15.5677e-03;
L=40e-03;

epir=3.4;
epi1=epir*epi0;
n1=sqrt(epi1/epi0);
mewr=1;
mew1=mewr*mew0;
cntr=1;
for fcntr= 0.1e09:0.10e09:12.5e09
            k0a=((2*pi*fcntr)/(3e08));
            k1a=((epi1/epi0)^(1/2))*(2*pi*fcntr/c0);
            lambda0=c0/fcntr;
            omega=2*pi*fcntr;
            freq(cntr)=fcntr;
            rsq(cntr)=((2*(b+L*tand(alpha))/lambda0)^2)*(pi^2)*(mewr*epir-1);
            Vfn = @(W) sqrt(rsq(cntr) - W.^2 );
            syms W
            Wa = solve((epir*(diff(besselj(1,Vfn(W)))/(Vfn(W)*besselj(1,Vfn(W))))-(diff(besselk(1,W))/(W*besselk(1,W))))*...
                (mewr*(diff(besselj(1,Vfn(W)))/(Vfn(W)*besselj(1,Vfn(W))))-(diff(besselk(1,W))/(W*besselk(1,W))))-((Vfn(W)^2+W^2)*...
                (Vfn(W)^2+mewr*epir*W^2))/(Vfn(W)^4*W^4), W );
            Va=Vfn(Wa);
            Vb(cntr)=double(Wa);
            Wb(cntr)=double(Va);
            beta(cntr)=sqrt(omega^2*mew1*epi1-(Vb(cntr)/(b+L*tand(alpha)))^2);
            beta2(cntr)=sqrt(omega^2*mew0*epi0+(Wb(cntr)/(b+L*tand(alpha)))^2);
            beta0(cntr)=beta(cntr)/k0a;
            beta20(cntr)=beta2(cntr)/k0a;
            cntr=cntr+1;
            fcntr
end
Getting error: "Warning: Explicit solution could not be found. 
> In solve at 179"
Also, tried : eqns=[(epir*(diff(besselj(1,V))/(V*k0a*besselj(1,V)))-(diff(besselk(1,W))/(W*k0a*besselk(1,W))))*...
                    (mewr*(diff(besselj(1,V))/(V*k0a*besselj(1,V)))-(diff(besselk(1,W))/(W*k0a*besselk(1,W))))==((V^2+W^2)*(V^2+mewr*epir*W^2))/(V^4*W^4), ...
                    ((2*(b+L*tand(alpha))/lambda0)^2)*(pi^2)*(mewr*epir-1)==(V^2+W^2)];            
                    [Va, Wa]=solve(eqns,[V W]);
Getting error: "Warning: 4 equations in 2 variables. 
> In C:\Program Files\MATLAB\R2013a\toolbox\symbolic\symbolic\symengine.p>symengine at 56
  In mupadengine.mupadengine>mupadengine.evalin at 97
  In mupadengine.mupadengine>mupadengine.feval at 150
  In solve at 170
  In propa_HE11_Chen_solve at 51 
Warning: Explicit solution could not be found. 
> In solve at 179"
Please help! Can't seem to use solve. The roots can be complex. Tried fzero,fsolve and Muller. But they all failed terribly.