k=0.4; alpha=0.5; tau=50; h=30;
% Define the values within the interval for gamma
minValue = 2*h/pi;
maxValue = 100;
gamma =linspace(minValue, maxValue, 1000);
%Define the left hand side of the above equation
func = sqrt(1-(k./gamma).^2)-((k*alpha*tau)./(gamma).^2)-sin(h./gamma);
% plot
plot(1./gamma, func)
