function [dF ] = dfun( X )

% 牛顿拉普森迭代

syms A B F1 F2 F3 F4 F5 F6

A=X(1);

B=X(2);

F1=X(3);

F2=X(4);

F3=X(5);

F4=X(6);

F5=X(7);

F6=X(8);

F=fun(X);

dF=[-4.39E-07/(1-B^2)+1.40E-06*B/(1-B^2),(-4.39E-07*2*A*B+1.40E-06*A*(1+B^2))/(1-B^2)^2,-1,0,0,0,0,0

-5.72E-06/(1-B^2)+2.44E-06*B/(1-B^2),(-5.72E-06*2*A*B+2.44E-06*A*(1+B^2))/(1-B^2)^2,0,0,0,0,0,0

1.44E-06/(1-B^2)-2.40E-06*B/(1-B^2),(1.44E-06*2*A*B-2.40E-06*A*(1+B^2))/(1-B^2)^2,0,-1,0,0,0,0

1.82E-07/(1-B^2)-1.07E-06*B/(1-B^2),(1.82E-07*2*A*B-1.07E-06*A*(1+B^2))/(1-B^2)^2,0,0,0,0,0,0

-7.91E-07/(1-B^2)+9.21E-06*B/(1-B^2),(-7.91E-07*2*A*B+9.21E-06*A*(1+B^2))/(1-B^2)^2,0,0,0,0,0,0

2.81E-06/(1-B^2)-1.56E-06*B/(1-B^2),(2.81E-06*2*A*B-1.56E-06*A*(1+B^2))/(1-B^2)^2,0,0,0,0,0,0

1.38E-06/(1-B^2)-9.80E-06*B/(1-B^2),(1.38E-06*2*A*B-9.80E-06*A*(1+B^2))/(1-B^2)^2,0,0,0,0,0,0

-4.62E-06/(1-B^2)+1.05E-06*B/(1-B^2),(-4.62E-06*2*A*B+1.05E-06*A*(1+B^2))/(1-B^2)^2,0,0,0,0,0,0

-1.21E-06/(1-B^2)+3.41E-06*B/(1-B^2),(-1.21E-06*2*A*B+3.41E-06*A*(1+B^2))/(1-B^2)^2,0,0,-1,0,0,0

8.39E-08/(1-B^2)+1.94E-06*B/(1-B^2),(8.39E-08*2*A*B+1.94E-06*A*(1+B^2))/(1-B^2)^2,0,0,0,-1,0,0

-3.85E-07/(1-B^2)-1.82E-06*B/(1-B^2),(-3.85E-07*2*A*B-1.82E-06*A*(1+B^2))/(1-B^2)^2,0,0,0,0,-1,0

7.26E-06/(1-B^2)-2.80E-06*B/(1-B^2),(7.26E-06*2*A*B-2.80E-06*A*(1+B^2))/(1-B^2)^2,0,0,0,0,0,-1];

dF=dF';

end