th1 = dtor(70);

n1 = 1; n2 = 1.6;

th2 = asin(1/n2*sin(th1));