function [ior,muam,iorr, iori, Cp,density,thermalconductivity,thermaldiffusivity,bulkmodulus,thermalexpansion,youngmodulus,poissonratio]=...
    melanin(lambda);
%
% ior=melanin(lambda); 
% [ior,muam,iorr, iori, Cp,density,thermalconductivity,...
% thermaldiffusivity,bulkmodulus,thermalexpansion,youngmodulus,...
% poissonratio]=melanin(lambda); 
%
% melanin properties
%
% by Chuck DiMarzio
%    Northeastern University
%    February 2012
%
%  Mechanical properties added by Jason Kellicker, January 2013
%  Also rearranged the order of the returned parameters
%
% Melanin index of refraction (real part)
% from Cosmetics Applications of Laser & Light-Based Systems
% By Gurpreet S. Ahluwalia
% Page 90, Eqq. 3.11
%
% Imaginary part computed from mu_a taken from Jacques, McAuliffe, Photochem
% and Photobiol. 53, No. 6, 1991. 
%
% Alternative equation from kurtz1986 is
% imag(n)=0.861*(1-0.990*(lambda*1e6)).^2
% This has a very different magnitude and shape, and has the
% disadvantage of going to zero at 1/0.990 microns. 
%
% lambda in m
%
% returns ior as complex index
% Cp specific heat, J/(kg*K)
% density, kg/m3
% thermalconductivity= 0.209; % W/(m*K)
% thermaldiffusivity=0.14/1000^2; % m2 /s
% 
% ior Valid from 350 to 800 nm

% Real part 
lambdanm=lambda*1e9;
iorr=1.6840...  %Real Part
    -1.8723e4./lambdanm.^2+1.0964e10./lambdanm.^4-8.6484e14./lambdanm.^6;
mua  =  (6.6E11)*(lambdanm.^(-3.33));  %1/cm
%mua  =  (6.6E11)*(lambdanm.^(-3.33))*10;  %1/cm - JK - for 10x absorption FDTD simulation
muam = mua*10^2; % 1/m
% Above from http://omlc.ogi.edu/news/jan98/skinoptics.html 
iori=(mua*100).*(lambda)/2/pi; % *100 to convert to m^{-1}
ior=iorr+1i*iori;
% density=1.32; % g/cm3.
density=1.32/1000*100^3; % kg/m3.
thermalconductivity= 0.209; % W/(m*K)
%thermaldiffusivity=0.14; % mm2 /s
thermaldiffusivity=0.14/1000^2; % m2 /s
Cp=1131; % J/(kg*K)n

%%% Below added Fri Apr  6 14:45:35 2012 from faraggi2005 (JBO Nov/Dec
%2005, Vol 10.  No 6. Pg. 064029-4.
bulkmodulus=29.4e9; % Pascals
thermalexpansion=2.98e-5; % Per Kelvin

%%% JK - 8/5/12 - Moses 2006 - Young's modulus. Poisson ratio calculated
%%% from Young's Modulus and Bulk Modulus.

youngmodulus = 1.5e9; %Pascals
poissonratio = (3*bulkmodulus - youngmodulus)/(6*bulkmodulus); %Unitless
specificvolume = 1/density;
volcoeffexp = thermalexpansion*3;  %Assumption from linear coefficient.
