Files
optical-trap-forces/Matlab/reflectivity.m
2023-03-20 22:50:56 +04:00

6 lines
252 B
Matlab

function R = reflectivity(th, psi, n1, n2)
theta_refl = th_r(th, n1, n2);
R = (tan(th - theta_refl).^2 ./ tan(th + theta_refl).^2) .* ...
cos(psi).^2 + ...
(sin(th - theta_refl).^2 ./ sin(th + theta_refl).^2) .* sin(psi).^2;
end