MATLAB: Translate equation to matlab code

equationsnell law

hellow I would like to translate this equation to matlab code, I am new to matlab language and syntax and would very much appreciate help.

Best Answer

s2_hat= n1/n2*(N_hat*(-N_hat*s1_hat))-N_hat*sqrt(1-(n1/n2)^2*(N_hat*s1_hat)*(N_hat*s1_hat))

Related Solutions

MATLAB: Solve on vector equation

Try this (solution exists not for any s1 vector)

function main
nt = 1;
ni = 2;
s1 = [1 2 2];
s1 = s1/norm(s1);
    function y = F(Nx)
        p = cross(-Nx,s1);
        f = nt/ni*cross(Nx,p) - Nx*sqrt(1-(nt/ni)^2*dot(p,p));
        y = f' - [1;0;0];
    end
N = fsolve(@F,[1 1 1]);
p = cross(-N,s1);
t = [s1
    N*sqrt(1-(nt/ni)^2*dot(p,p))
    nt/ni*cross(N,p)];
v0 = zeros(3,1);
quiver3(v0,v0,v0,t(:,1),t(:,2),t(:,3),1)        % all vectors together
hold on
quiver3(t(2,1),t(2,2),t(2,3),1,0,0, 1,'r')      % vector [1 0 0]
% quiver3(0,0,0,N(1),N(2),N(3),'g')
hold off
text(s1(1),s1(2),s1(3),'s1 vector')
view(0,0)
axis equal
end

MATLAB: Solving System of Equations – Complex

Try this:


kp1 = 2;
kp2 = 3;
p   = 1;
syms xco2 xco xo2 xo
eqns = [(xo2^(1/2)* p^(1/2)*xco)/xco2 == kp1, 
        (xo^2*p)/xo2 == kp2, 
         xco2 + xco + xo2 + xo == 1,
        (xco2 + xco)/(xco + xo + 2*xo2) == 1/2] ; 
sol = solve(eqns,[xco2, xco, xo2, xo])

Best Answer

Related Solutions

MATLAB: Solve on vector equation

MATLAB: Solving System of Equations – Complex

Related Question