use PDL;

		my ($mean,$std_dev,$median,$min,$max,$adev,$rms) = stats($source);
		
		my $skewness = sclr(sum(   
		
						 ($source - $mean)**3
		
							) / (
	
						 (nelem($source)-1) * $std_dev**3 
						
					   ));
		
		my $exs_kurtosis = sclr(sum(   
		
						 ($source - $mean)**4
		
							) / (
	
						 (nelem($source)-1) * $std_dev**4 
						
					   ) -3);
		
		my $JB	= 
					( nelem($source) / 6 )
					*
					(
						$skewness**2
						+ ($exs_kurtosis**2 / 4)
					)
				  ;