% CALCOLO DEL DETERMINANTE DELLA PRIMA MATRICE DI DISPERSIONE

A(1,1)=exp(-h(1)*d2);

A(1,2)=-cos(h(2)*d2);

A(1,3)=-sin(h(2)*d2);

A(2,1)=h(1)*exp(-h(1)*d2);

A(2,2)=-h(2)*sin(h(2)*d2);

A(2,3)=h(2)*cos(h(2)*d2);

A(3,2)=cos(h(2)*g/2);

A(3,3)=sin(h(2)*g/2);

A(3,4)=-exp(h(3)*g/2);

A(3,5)=-exp(-h(3)*g/2);

A(4,2)=h(2)*sin(h(2)*g/2);

A(4,3)=-h(2)*cos(h(2)*g/2);

A(4,4)=h(3)*exp(h(3)*g/2);

A(4,5)=-h(3)*exp(-h(3)*g/2);

A(5,4)=exp(-h(3)*g/2);

A(5,5)=exp(h(3)*g/2);

A(5,6)=-cos(h(4)*g/2);

A(5,7)=sin(h(4)*g/2);

A(6,4)=-h(3)*exp(-h(3)*g/2);

A(6,5)=h(3)*exp(h(3)*g/2);

A(6,6)=h(4)*sin(h(4)*g/2);

A(6,7)=h(4)*cos(h(4)*g/2);

A(7,6)=cos(h(4)*d4);

A(7,7)=-sin(h(4)*d4);

A(7,8)=-exp(-h(5)*d4);

A(8,6)=-h(4)*sin(h(4)*d4);

A(8,7)=-h(4)*cos(h(4)*d4);

A(8,8)=h(5)*exp(-h(5)*d4);

detA=det(A);

function detA=mdisp2(h,d2,d4,g);