I have solved two equations with two unknowns and confined my answer in just ten answers. The unknowns are n and k, and I want to make histogram for them separately, but don't know how to separate the n and k?
My code is like this:
d = 2/100;
lambda = (1.67)*10^(-4);
alpha = Exp[-2*Pi*k*d/lambda];
delta = 4*Pi*n*d/lambda;
R = ((n - 1)^2 + k^2)/((n + 1)^2 + k^2);
Rp = 0.133;
Tp = 0.601;
pts = Table[
FindRoot[{R (1 + (alpha)^4 - 2*(alpha)^2*Cos[delta])/(1 -
2*(alpha)^2*R*Cos[delta] + (alpha)^4*R^2) ==
Rp, ((alpha)^2*(1 - R)^2)/(1 -
2*(alpha)^2*R*Cos[delta] + (alpha)^4*R^2) == Tp}, {n,
RandomReal[{1.5, 1.7}]}, {k, RandomReal[{0, 0.02}]}], {10}]
Grid[pts]
and what I got is
{n -> 1.664, k -> 0.0002223}, {n -> 1.609,
k -> 0.0002249}, {n -> 1.669, k -> 0.0002221}, {n -> 1.597,
k -> 0.0002255}, {n -> 1.568, k -> 0.0002268}, {n -> 1.592,
k -> 0.0002257}, {n -> 1.646, k -> 0.0002232}, {n -> 1.664,
k -> 0.0002223}, {n -> 1.568, k -> 0.0002268}, {n -> 1.663,
k -> 0.0002224}}



{n, k} /. pts– Edmund Mar 01 '16 at 21:00Histogram[n /. pts]– george2079 Mar 01 '16 at 21:47