Question: optimization

How can i solve this complicated problem????restart:
R:=3:
kappa:=2.4048/R:
a:=-mu^2/kappa^2:
sigma1:=Int(t^2*BesselJ(0,t)*BesselJ(1,t)*exp(a*t^2),t=0..2.4048):
sigma2:=Int(t^3*BesselJ(0,t)^2*exp(a*t^2),t=0..2.4048):
sigma3:=int(t*BesselJ(0,t)^2*exp(a*t^2),t=0..2.4048):
sigma4:=Int(t*BesselJ(0,t)^2*BesselJ(0,(q*t)/kappa)*exp(a*t^2),t=0..2.4048):
s[q]:=(sigma4/sigma3)^2:
epsilon:=lambda^2/2+kappa^2+mu^2*(2-(2*sigma1-(mu/kappa)2*sigma2)/sigma3):
h:=(s[q]*exp(q^2/(2*lambda^2))*erfc(q/(sqrt(2)*lambda))):

x:=Int(h,y=0..infinity):
 E:=epsilon-(x):
with (Optimization);
  Minimize(E,lambda=0..10000,mu=0..10000);
       

 

Please Wait...