Help with plotting reflection Coefficient magnitude vs frequency for a single stub shunt (short circuit)

I am trying to use matlab to replicate this graph from the book.

The load impedance is Zl = 60-j80 ohms

The load is assumed to be matched at 2 GHz

Zo = 50 ohms

Range is from 1 to 3 GHz

The normalized impedance is zl = 1.2 - j1.6

d1 = 0.110 gamma

d2 = 0.260 gamma

y1 = 1.00 + j1.47

y2 = 1.00 - j1.47

[Yuehan Zhang]

hi there, just wanna know if you know how this one can be plotted? I am struggling with this question as well.