// Cell Distribution Parameters taken from (Lee 2019) and (Masri 2021)
// Polynomial in the form of:
// D = cell_c1*exp(cell_lambda1*d) + cell_c2*exp(cell_lambda2*d) + cell_c3*exp(cell_lambda3*d)
// Where D is cell density in (cells/mm^2), c_n and lambda_n are dimensionless coefficients and d is eccentricity from the fovea in (mm)
// Assuming constant radial cell densities

rod_c1 = 585500
rod_lambda1 = -0.1388
rod_c2 = -598900
rod_lambda2 = -0.2998
rod_c3 = 0
rod_lambda3 = 0

cone_c1 = 367300
cone_lambda1 = -7.828
cone_c2 = -200000
cone_lambda2 = -2000
cone_c3 = 20340
cone_lambda3 = -0.2164

horizontal_c1 = 20510
horizontal_lambda1 = -0.3850
horizontal_c2 = -42460
horizontal_lambda2 = -4.196
horizontal_c3 = 2672
horizontal_lambda3 = -0.01097

rod_bipolar_c1 = 291700
rod_bipolar_lambda1 = -0.2529
rod_bipolar_c2 = -293200
rod_bipolar_lambda2 = -0.2803
rod_bipolar_c3 = 0
rod_bipolar_lambda3 = 0

ON_bipolar_c1 = 178800
ON_bipolar_lambda1 = -0.4755
ON_bipolar_c2 = -183600
ON_bipolar_lambda2 = -0.5618
ON_bipolar_c3 = 8894
ON_bipolar_lambda3 = -0.01832

OFF_bipolar_c1 = 420800
OFF_bipolar_lambda1 = -1.225
OFF_bipolar_c2 = -455100
OFF_bipolar_lambda2 = -1.387
OFF_bipolar_c3 = 12300
OFF_bipolar_lambda3 = -0.06986

AII_amacrine_c1 = 129700
AII_amacrine_lambda1 = -1.114
AII_amacrine_c2 = -139200
AII_amacrine_lambda2 = -1.217
AII_amacrine_c3 = 2553
AII_amacrine_lambda3 = -0.07214

// Total equation for full cell types
pr_c1 = 585500
pr_lambda1 = -0.1388
pr_c2 = -598900
pr_lambda2 = -0.2998
pr_c3 = 0
pr_lambda3 = 0
pr_c4 = 367300
pr_lambda4 = -7.828
pr_c5 = -200000
pr_lambda5 = -2000
pr_c6 = 20340
pr_lambda6 = -0.2164

bipolar_c1 = 291700
bipolar_lambda1 = -0.2529
bipolar_c2 = -293200
bipolar_lambda2 = -0.2803
bipolar_c3 = 0
bipolar_lambda3 = 0
bipolar_c4 = 178800
bipolar_lambda4 = -0.4755
bipolar_c5 = -183600
bipolar_lambda5 = -0.5618
bipolar_c6 = 8894
bipolar_lambda6 = -0.01832
bipolar_c7 = 420800
bipolar_lambda7 = -1.225
bipolar_c8 = -455100
bipolar_lambda8 = -1.387
bipolar_c9 = 12300
bipolar_lambda9 = -0.06986