Advanced Search

Question: fsolve-not giving solutions

Question:fsolve-not giving solutions

arunousephthoma 5 Maple
complex fsolve
February 07 2013
1

hi all,

> e0 := -21+.7*I;
                          -21. + 0.7 I
> e2 := e0;
                          -21. + 0.7 I
> e1 := 3.9;
                              3.9
> e3 := 1;
                               1
> t1 := 10*10^(-9);
                               1    
                           ---------
                           100000000
> t2 := 50*10^(-9);
                               1    
                            --------
                            20000000
> lambda = 780*10^(-9);
                                       39   
                       lambda = --------
                                      50000000
> k0 := (2*3.14)/lambda;
                                        6.28 
                                        ------
                                         lambda
> g0 := sqrt(e0*k0^2-x^2);
> g1 := sqrt(e1*k0^2-x^2);

> g2 := g0;

> g3 := sqrt(e3*k0^2-x^2);

> theta1 := 0+I*g1*t1;
> theta2 := 0+I*g2*t2;

> m := (e2*g3-e3*g2)*(e1*g2+e2*g1)*(e0*g1-e1*g0)*e^(theta1+theta2)+(e2*g3-e3*g2)*(e1*g2-e2*g1)*(e0*g1+e1*g0)*e^(theta2-theta1)+(e2*g3+e3*g2)*(e1*g2-e2*g1)*(e0*g1-e1*g0)*e^(theta1-theta2)+(e2*g3+e3*g2)*(e1*g2+e2*g1)*(e0*g1+e1*g0)*e^(-theta1-theta2)=0;
 
gives
 
m := ((-21.+.7*I)*sqrt(39.4384/lambda^2-x^2)-sqrt((-828.2064+27.60688*I)/lambda^2-x^2))*(3.9*sqrt((-828.2064+27.60688*I)/lambda^2-x^2)+(-21.+.7*I)*sqrt(153.80976/lambda^2-x^2))*((-21.+.7*I)*sqrt(153.80976/lambda^2-x^2)-3.9*sqrt((-828.2064+27.60688*I)/lambda^2-x^2))*2.72^((1/100000000*I)*sqrt(153.80976/lambda^2-x^2)+(1/20000000*I)*sqrt((-828.2064+27.60688*I)/lambda^2-x^2))+((-21.+.7*I)*sqrt(39.4384/lambda^2-x^2)-sqrt((-828.2064+27.60688*I)/lambda^2-x^2))*(3.9*sqrt((-828.2064+27.60688*I)/lambda^2-x^2)+(21.-.7*I)*sqrt(153.80976/lambda^2-x^2))*(3.9*sqrt((-828.2064+27.60688*I)/lambda^2-x^2)+(-21.+.7*I)*sqrt(153.80976/lambda^2-x^2))*2.72^((1/20000000*I)*sqrt((-828.2064+27.60688*I)/lambda^2-x^2)-(1/100000000*I)*sqrt(153.80976/lambda^2-x^2))+((-21.+.7*I)*sqrt(39.4384/lambda^2-x^2)+sqrt((-828.2064+27.60688*I)/lambda^2-x^2))*(3.9*sqrt((-828.2064+27.60688*I)/lambda^2-x^2)+(21.-.7*I)*sqrt(153.80976/lambda^2-x^2))*((-21.+.7*I)*sqrt(153.80976/lambda^2-x^2)-3.9*sqrt((-828.2064+27.60688*I)/lambda^2-x^2))*2.72^((1/100000000*I)*sqrt(153.80976/lambda^2-x^2)-(1/20000000*I)*sqrt((-828.2064+27.60688*I)/lambda^2-x^2))+((-21.+.7*I)*sqrt(39.4384/lambda^2-x^2)+sqrt((-828.2064+27.60688*I)/lambda^2-x^2))*(3.9*sqrt((-828.2064+27.60688*I)/lambda^2-x^2)+(-21.+.7*I)*sqrt(153.80976/lambda^2-x^2))^2*2.72^(-(1/100000000*I)*sqrt(153.80976/lambda^2-x^2)-(1/20000000*I)*sqrt((-828.2064+27.60688*I)/lambda^2-x^2))
 
>fsolve(m = 0, x, x = -10^(-4) .. 1.10^4, complex,maxsols=20);
 
is not giving me solutions. i tried with different ranges too.suggestions please.
thanks
Please Wait...