<?php
	# ------- The graph values in the form of associative array
	function DensitySp38($V,$s1,$s2,$s3,$s4,$s5,$s6,$s7,$s8,$s9,$s10,$s11,$s12,$s13,$s14,$s15,$s16,$s17,$s18,$s19,$s20,$s21,$s22,$s23,$s24,$s25,$s26,$s27,$s28,$s29,$s30,$s31,$s32,$s33,$s34,$s35,$s36) {
	
	$values=array("34" => 0.0,"35" => 0.0,"36" => 0.0,"37" => 0.0,"38" => $s1,"39" =>$s2,"40" => $s3,"41" => $s4,"42" => $s5,"43" => $s6,"44" => $s7,"45" => $s8,"46" => $s9,"47" => $s10,"48" => $s11,"49" => $s12,"50" => $s13,"51" => $s14,"52" => $s15,"53" => $s16,"54" => $s17,"55" => $s18,"56" => $s19,"57" => $s20,"58" => $s21,"59" => $s22,"60" => $s23,"61" => $s24,"62" => $s25,"63" => $s26,"64" => $s27,"65" => $s28,"66" => $s29,"67" => $s30,"68" => $s31,"69" => $s32,"70" => $s33, "71" => $s34,"72" => $s35,"73" => $s36,"74" => 0.0,"75" => 0.0,"76" => 0.0,"77" => 0.0,"78" => 0.0,"79" => 0.0,"80" => 0.0,"81" => 0.0,"82" => 0.0,"83" => 0.0);

	$SpectraC=array(0.00000000 ,0.00000000 ,0.00000000 ,0.00000000 ,0.00000000 ,0.00000000 ,0.00000000 ,0.00000000 ,0.00000000 ,0.00000000 ,0.00000000 ,0.00000000 ,0.00000000 ,0.00000000 ,0.00000000 ,0.00000000 ,0.00000000 ,0.00000000 ,0.00000000 ,0.00000000 ,0.00000000 ,0.00000000 ,0.00100000 ,0.05997911 ,0.44977985 ,29.99162519 ,100.00000000 ,84.91804750 ,54.95408739 ,25.00345362 ,10.00000000 ,5.00034535 ,1.49968484 ,0.50003453 ,0.29991625 ,0.14996848 ,0.05000345 ,0.01000000 ,0.00100000 ,0.00000000 ,0.00000000 ,0.00000000 ,0.00000000 ,0.00000000 ,0.00000000 ,0.00000000 ,0.00000000 ,0.00000000 ,0.00000000 ,0.00000000);

	$SpectraM=array(0.00000000 ,0.00000000 ,0.00000000 ,0.00000000 ,0.00000000 ,0.00000000 ,0.00000000 ,0.00000000 ,0.00000000 ,0.00000000 ,0.00000000 ,0.00000000 ,0.00000000 ,0.00100000 ,1.00000000 ,5.00034535 ,27.98981320 ,68.07693587 ,92.04495718 ,100.00000000 ,87.90225168 ,66.06934480 ,41.97589840 ,21.97859873 ,8.99497582 ,2.50034536 ,0.69984200 ,0.08994976 ,0.01000000 ,0.00100000 ,0.00000000 ,0.00000000 ,0.00000000 ,0.00000000 ,0.00000000 ,0.00000000 ,0.00000000 ,0.00000000 ,0.00000000 ,0.00000000 ,0.00000000 ,0.00000000 ,0.00000000 ,0.00000000 ,0.00000000 ,0.00000000 ,0.00000000 ,0.00000000 ,0.00000000 ,0.00000000);
	
	$SpectraY=array(0.00000000 ,0.01000000 ,0.01999862 ,0.10000000 ,0.29991625 ,1.49968484 ,5.99791076 ,16.98243652 ,39.99447498 ,59.97910763 ,82.03515443 ,93.97233106 ,100.00000000 ,97.05099672 ,84.91804750 ,65.01296903 ,39.99447498 ,17.98870915 ,5.00034535 ,0.19998619 ,0.03999447 ,0.00100000 ,0.00000000 ,0.00000000 ,0.00000000 ,0.00000000 ,0.00000000 ,0.00000000 ,0.00000000 ,0.00000000 ,0.00000000 ,0.00000000 ,0.00000000 ,0.00000000 ,0.00000000 ,0.00000000 ,0.00000000 ,0.00000000 ,0.00000000 ,0.00000000 ,0.00000000 ,0.00000000 ,0.00000000 ,0.00000000 ,0.00000000 ,0.00000000 ,0.00000000 ,0.00000000 ,0.00000000 ,0.00000000);	

	$SpectraV=array(0.00000000 ,0.00000000 ,0.00000000 ,0.00000000 ,0.00000000 ,0.00000000 ,0.01000000 ,0.02098940 ,0.08203515 ,0.27989813 ,0.64714262 ,1.23026877 ,2.21819642 ,3.81944271 ,6.57657837 ,10.99005839 ,18.87991349 ,32.58367010 ,50.35006088 ,66.83439176 ,80.35261222 ,90.57326009 ,97.49896377 ,100.00000000 ,97.49896377 ,90.36494737 ,79.79946873 ,67.14288529 ,53.82697825 ,39.17418771 ,27.10191632 ,17.29816359 ,10.30386120 ,5.61047976 ,3.09029543 ,1.54170045 ,0.79615935 ,0.41879357 ,0.21577444 ,0.10990058 ,0.05395106 ,0.02697739 ,0.01399587 ,0.01000000 ,0.00000000 ,0.00000000 ,0.00000000 ,0.00000000 ,0.00000000 ,0.00000000);
 

//-LOG10(SUMPRODUCT(AZ2:AZ51,AV2:AV51)/SUM(AV2:AV51))   =-LOG10(SUMPRODUCT(Sample,TT)/SUM(Den_StatusT_R))
$sPSum=0;
$sProduct=0;
if ($V=="V") {
for ( $i = 34; $i <= 83; $i += 1) {
$sProduct=$values[$i]*$SpectraV[$i-34];
$sPSum=$sProduct+$sPSum;}
$sumV=array_sum($SpectraV);}
elseif ($V=="C"){
for ( $i = 34; $i <= 83; $i += 1) {
$sProduct=$values[$i]*$SpectraC[$i-34];
$sPSum=$sProduct+$sPSum;}
$sumV=array_sum($SpectraC);}
elseif ($V=="M"){
for ( $i = 34; $i <= 83; $i += 1) {
$sProduct=$values[$i]*$SpectraM[$i-34];
$sPSum=$sProduct+$sPSum;}
$sumV=array_sum($SpectraM);}
elseif ($V=="Y"){
for ( $i = 34; $i <= 83; $i += 1) {
$sProduct=$values[$i]*$SpectraY[$i-34];
$sPSum=$sProduct+$sPSum;}
$sumV=array_sum($SpectraY);}
elseif ($V=="K"){
for ( $i = 34; $i <= 83; $i += 1) {
$sProduct=$values[$i]*$SpectraV[$i-34];
$sPSum=$sProduct+$sPSum;}
$sumV=array_sum($SpectraV);}

 //$sumduct=$values[40]; 
 
 $den=-log10($sPSum/$sumV);
 return round($den,4);
	}

?>