#ifndef OPTICS_H
#define OPTICS_H

double get_index(double p, double wavelength, double T);

#endif