L = 100;
H = 50;
D = sqrt(L*L-H*H);
k = 10;


dH = -30:0.1:130;

Fk = ... 
F  = ...

plot(dH,F); grid on; 