%photodiode responsivity

function [y]=photodiode(lambda);

wl=[400:50:800];
PD=[0.08, 0.17, 0.23, 0.3, 0.35, 0.4, 0.44, 0.48, 0.52];
y=interp1(wl,PD,lambda,'cubic')