<?php
	# ------- The graph values in the form of associative array
	function DensitySp($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) {
	
	$values=array("34" => 0.0,"35" => 0.0,"36" => 0.0,"37" => 0.0,"38" => 0.0,"39" => 0.0,"40" => $s1,"41" => $s2,"42" => $s3,"43" => $s4,"44" => $s5,"45" => $s6,"46" => $s7,"47" => $s8,"48" => $s9,"49" => $s10,"50" => $s11,"51" => $s12,"52" => $s13,"53" => $s14,"54" => $s15,"55" => $s16,"56" => $s17,"57" => $s18,"58" => $s19,"59" => $s20,"60" => $s21,"61" => $s22,"62" => $s23,"63" => $s24,"64" => $s25,"65" => $s26,"66" => $s27,"67" => $s28,"68" => $s29,"69" => $s30,"70" => $s31, "71" => 0.0,"72" => 0.0,"73" => 0.0,"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("34" => 0.0,"35" => 0.0,"36" => 0.0,"37" => 0.0,"38" => 0.0,"39" => 0.0,"40" => 0.0,"41" => 0.0,"42" => 0.0,"43" => 0.0,"44" => 0.0,"45" => 0.0,"46" => 0.0,"47" => 0.0,"48" => 0.0,"49" => 0.0,"50" =>0.0,"51" => 0.0,"52" => 0.0,"53" => 0.0,"54" => $s15,"55" => 0.0,"56" => 0.0010 ,"57" => 0.0600 ,"58" => 0.4498,"59" => 29.9916 
,"60" => 100.0000 
,"61" => 84.9180 
,"62" => 54.9541 
,"63" => 25.0035 
,"64" => 10.0000 
,"65" => 5.0003 
,"66" => 1.4997 
,"67" => 0.5000 
,"68" => 0.2999 
,"69" => 0.1500 
,"70" => 0.0500 
, "71" => 0.0100 
,"72" =>0.0010 
,"73" => 0.0,"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);

	$SpectraM=array("34" => 0.0,"35" => 0.0,"36" => 0.0,"37" => 0.0,"38" => 0.0,"39" => 0.0,"40" => 0.0,"41" => 0.0,"42" => 0.0,"43" => 0.0,"44" => 0.0,"45" => 0.0,"46" => 0.0,"47" => 0.0010 
,"48" => 1.0000 
,"49" => 5.0003 
,"50" =>27.9898 
,"51" => 68.0769 
,"52" => 92.0450 
,"53" => 100.0000 
,"54" => 87.9023 
,"55" => 66.0693 
,"56" => 41.9759 
 ,"57" => 21.9786 
 ,"58" => 8.9950 
,"59" => 2.5003
,"60" => 0.6998  
,"61" => 0.0899  
,"62" => 0.0100  
,"63" => 0.0010  
,"64" => 0.0 
,"65" => 0.0 
,"66" => 0.0 
,"67" => 0.0 
,"68" => 0.0 
,"69" => 0.0 
,"70" => 0.0 
, "71" => 0.0 
,"72" =>0.0 
,"73" => 0.0,"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);
	
	$SpectraY=array("34" => 0.0,"35" => 0.0100 
,"36" => 0.0200 
,"37" => 0.1000 
,"38" => 0.2999 
,"39" => 1.4997 
,"40" => 5.9979 
,"41" => 16.9824 
,"42" => 39.9945 
,"43" => 59.9791 
,"44" => 82.0352 
,"45" => 93.9723 
,"46" => 100.0000 
,"47" => 97.0510  
,"48" => 84.9180  
,"49" => 65.0130  
,"50" =>39.9945  
,"51" => 17.9887  
,"52" => 5.0003 
,"53" => 0.2000  
,"54" => 0.0400  
,"55" => 0.0010  
,"56" => 0.0000  
 ,"57" => 0.00 
 ,"58" => 0.00 
,"59" => 0.00
,"60" => 0.00  
,"61" => 0.00  
,"62" => 0.00  
,"63" => 0.00  
,"64" => 0.0 
,"65" => 0.0 
,"66" => 0.0 
,"67" => 0.0 
,"68" => 0.0 
,"69" => 0.0 
,"70" => 0.0 
, "71" => 0.0 
,"72" =>0.0 
,"73" => 0.0,"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);	

	$SpectraV=array("34" => 0.0,"35" => 0.0000 
,"36" => 0.0 
,"37" => 0.0 
,"38" => 0.0 
,"39" => 0.0 
,"40" => 0.0100  
,"41" => 0.0210  
,"42" => 0.0820  
,"43" => 0.2799  
,"44" => 0.6471  
,"45" => 1.2303  
,"46" => 2.2182  
,"47" => 3.8194   
,"48" => 6.5766   
,"49" => 10.9901   
,"50" =>18.8799   
,"51" => 32.5837   
,"52" => 50.3501  
,"53" => 66.8344   
,"54" => 80.3526   
,"55" => 90.5733   
,"56" => 97.4990   
 ,"57" => 100.0000  
 ,"58" => 97.4990  
,"59" => 90.3649 
,"60" => 79.7995   
,"61" => 67.1429   
,"62" => 53.8270   
,"63" => 39.1742   
,"64" => 27.1019  
,"65" => 17.2982  
,"66" => 10.3039  
,"67" => 5.6105  
,"68" => 3.0903  
,"69" => 1.5417  
,"70" => 0.7962  
, "71" => 0.4188  
,"72" =>0.2158  
,"73" => 0.1099 
,"74" => 0.0540 
,"75" => 0.0270 
,"76" => 0.0140 
,"77" => 0.0100 
,"78" => 0.0,"79" => 0.0,"80" => 0.0,"81" => 0.0,"82" => 0.0,"83" => 0.0);
 

//-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];
$sPSum=$sProduct+$sPSum;}
$sumV=array_sum($SpectraV);}
elseif ($V=="C"){
for ( $i = 34; $i <= 83; $i += 1) {
$sProduct=$values[$i]*$SpectraC[$i];
$sPSum=$sProduct+$sPSum;}
$sumV=array_sum($SpectraC);}
elseif ($V=="M"){
for ( $i = 34; $i <= 83; $i += 1) {
$sProduct=$values[$i]*$SpectraM[$i];
$sPSum=$sProduct+$sPSum;}
$sumV=array_sum($SpectraM);}
elseif ($V=="Y"){
for ( $i = 34; $i <= 83; $i += 1) {
$sProduct=$values[$i]*$SpectraY[$i];
$sPSum=$sProduct+$sPSum;}
$sumV=array_sum($SpectraY);}
elseif ($V=="K"){
for ( $i = 34; $i <= 83; $i += 1) {
$sProduct=$values[$i]*$SpectraV[$i];
$sPSum=$sProduct+$sPSum;}
$sumV=array_sum($SpectraV);}

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

?>