Построить график интенсивности ДОЭ вдоль оптической оси и выполнить анализ хроматизма

Выполнить анализ хроматизма линзы, вычислив зависимости распределения интенсивности вдоль оптической оси I(y) на длинах волн ?=400,500,600,700 нм при y=(0..1,5)f. Предложить методы коррекции хроматической аберрации дифракционных линз.
Исходные данные:
Ширина линзы L=6 мм
Фокусное расстояние f=55 мм
Рабочая длина волны ?=635 нм
Падение оптического излучения – нормальное
Распределение интенсивности падающего пучка – равномерное
Линза изготовлена на поверхности стеклянной подложки с показателем преломления n=1,45.
Программы для расчета рельефно-фазовой функции и распределения интенсивности в фокальной плоскости я написала. Здесь немного в тупик встала, а времени мало до сдачи. Надо писать в Scilab в соответствии с моей программой расчета рельефно-фазовой функции