// Calculate the threshold wavelength
clc;
Th_wavelength=1.24*10^-6/1.8'
disp(Th_wavelength,'Threshold wavelength (m)')