solver1.m

One program to simulation light passing through a glass sphere. It is interesting and can be used i

function F=solver1(Pw,Pc,delC,L,beta)
F=[dot((Pw-Pc),delC)-norm(Pw-Pc)*norm(delC)*cos(beta);...
        dot((Pw-Pc),cross(L,delC));...
        norm(Pw-Pc)-1];