Найдите максимальную площадь (MatLab)

Рассмотрим код ниже (MatLab):

w = 0 : 0.0001 : 9.4978;
a = [1    11    46    95   109    74    24];
b = [-1 3 4 3 1];
mu = 1;
a0 = a(7) ;a1 = a(6) ;a2 = a(5); a3 = a(4) ; a4 = a(3) ; a5 = a(2); a6 = a(1);
b0 = b(5);b1 = b(4);b2 = b(3) ; b3 = b(2); b4 = b(1) ;
De = -a6*w.^6 + a4*w.^4 - a2*w.^2 + a0;
Do = a5*w.^4 - a3*w.^2 + a1;
Ne = b4*w.^4 - b2*w.^2 + b0;
No = -b3*w.^2 + b1;
T = 0.01;
e = real((1i*w).^mu);
f = imag((1i*w).^mu);
A = Ne.*cos(T*w) + w.*No.*sin(T*w);  
B = e.*(Ne.*cos(T*w) + w.*No.*sin(T*w)) - f.*(w.*No.*cos(T*w) - Ne.*sin(T*w));
C = w.*No.*cos(T*w) - Ne.*sin(T*w);
D = e.*(w.*No.*cos(T*w) - Ne.*sin(T*w)) + f.*(Ne.*cos(T*w) + w.*No.*sin(T*w));
Kp = (-De.*D + w.*Do.*B)./(f.*(Ne.^2 + w.^2.*No.^2));
Kd = (-w.*Do.*A + De.*C)./(f.*(Ne.^2 + w.^2.*No.^2));
figure
plot(Kp,Kd)
line([-24 -24],[-2.24 9.813])

Запустив код, мы получили этот рисунок:

Я хочу нарисовать касательные линии на указанной части кривой (красная часть, w принадлежит [0.6342,0.9985]):

после этого моя цель - найти максимальную площадь направленной внутрь полуплоскости, определяемой этой линией, и кривой между всеми возможными областями, полученными касательной линией (например):

другой пример с другой касательной в другой точке:

и мы можем заключить, что первая область больше второй. Этот подход должен работать для всех точек в красной части.

Как я могу сделать это с помощью MatLab?

Надеюсь, мой вопрос понятен. Любая идея будет оценена.