metal clad optical waveguides

hi, i would like to ask for help for matlab codes for my project which is metal clad optical waveguide polarizer.I need the coding to plot the 2-D graph of log α(1/cm) vs ω(μm).The unit ω represents core layer thickness.i anyone has the coding or know where to get the related coding for it,let me know please.thank you