P3=subs(4*cos(2*l*((m/2 - n/2)^2 + k^2)^(1/2))^2*cos(k*l)^2*sin(k*l)^2 + (k^2*sin(2*l*((m/2 - n/2)^2 + k^2)^(1/2))^2*cos(k*l)^4)/((m/2 - n/2)^2 + k^2) + (k^2*sin(2*l*((m/2 - n/2)^2 + k^2)^(1/2))^2*sin(k*l)^4)/((m/2 - n/2)^2 + k^2) - (2*k^2*sin(2*l*((m/2 - n/2)^2 + k^2)^(1/2))^2*cos(k*l)^2*sin(k*l)^2)/((m/2 - n/2)^2 + k^2) - (4*k*cos(2*l*((m/2 - n/2)^2 + k^2)^(1/2))*sin(2*l*((m/2 - n/2)^2 + k^2)^(1/2))*cos(k*l)*sin(k*l)^3)/((m/2 - n/2)^2 + k^2)^(1/2) + (4*k*cos(2*l*((m/2 - n/2)^2 + k^2)^(1/2))*sin(2*l*((m/2 - n/2)^2 + k^2)^(1/2))*cos(k*l)^3*sin(k*l))/((m/2 - n/2)^2 + k^2)^(1/2),{(m/2 - n/2),k,l},{1/lambda,1960,1800})
lambda=1.30:0.01:1.60;
ezplot(P3,lambda)

最后有K，l，lambda 三个变量，我的目的是想找到K，l的值使的P3的值随lambda变化不明显，几乎与lambda无关，lambda的值在1310nm到1600nm，我上面用subs函数是假设K，l为某一个定值，可是这样手动的改，找不到最优解，也没有说服力，有没有高手帮帮我该怎么办

可以让电脑自己找这些数值 然后给lambda数值 最后一个和第一个变化不大