FM / jpgraph

/var/www/doncode.com.ua/www/obmen/include/jpgraph/

/var/www/doncode.com.ua/www/obmen/include/jpgraph

Name	Size	Mode	Actions
lang/	-	0775	rm
CHANGELOG	89389	0664	edit dl rm
flags.dat	982587	0664	edit dl rm
flags_thumb35x35.dat	214310	0664	edit dl rm
flags_thumb60x60.dat	375469	0664	edit dl rm
flags_thumb100x100.dat	687530	0664	edit dl rm
gd_image.inc.php	68914	0664	edit dl rm
imgdata_balls.inc.php	53731	0664	edit dl rm
imgdata_bevels.inc.php	4500	0664	edit dl rm
imgdata_diamonds.inc.php	8329	0664	edit dl rm
imgdata_pushpins.inc.php	27889	0664	edit dl rm
imgdata_squares.inc.php	6733	0664	edit dl rm
imgdata_stars.inc.php	6578	0664	edit dl rm
jpg-config.inc.php	9243	0664	edit dl rm
jpgraph.php	185282	0664	edit dl rm
jpgraph_antispam-digits.php	11963	0664	edit dl rm
jpgraph_antispam.php	39948	0664	edit dl rm
jpgraph_bar.php	30740	0664	edit dl rm
jpgraph_canvas.php	3071	0664	edit dl rm
jpgraph_canvtools.php	14000	0664	edit dl rm
jpgraph_date.php	16501	0664	edit dl rm
jpgraph_errhandler.inc.php	7518	0664	edit dl rm
jpgraph_error.php	4360	0664	edit dl rm
jpgraph_flags.php	12801	0664	edit dl rm
jpgraph_gantt.php	130736	0664	edit dl rm
jpgraph_gb2312.php	120006	0664	edit dl rm
jpgraph_gradient.php	12975	0664	edit dl rm
jpgraph_iconplot.php	5333	0664	edit dl rm
jpgraph_imgtrans.php	7319	0664	edit dl rm
jpgraph_led.php	10655	0664	edit dl rm
jpgraph_line.php	17806	0664	edit dl rm
jpgraph_log.php	8542	0664	edit dl rm
jpgraph_mgraph.php	11890	0664	edit dl rm
jpgraph_pie.php	40020	0664	edit dl rm
jpgraph_pie3d.php	25728	0664	edit dl rm
jpgraph_plotband.php	18461	0664	edit dl rm
jpgraph_plotmark.inc.php	14262	0664	edit dl rm
jpgraph_polar.php	23007	0664	edit dl rm
jpgraph_radar.php	19682	0664	edit dl rm
jpgraph_regstat.php	5668	0664	edit dl rm
jpgraph_scatter.php	6603	0664	edit dl rm
jpgraph_stock.php	5060	0664	edit dl rm
jpgraph_ttf.inc.php	10971	0664	edit dl rm
jpgraph_utils.inc.php	15633	0664	edit dl rm

Edit: /var/www/doncode.com.ua/www/obmen/include/jpgraph/jpgraph_regstat.php (5668B)

<?php 
/*=======================================================================
// File:	JPGRAPH_REGSTAT.PHP
// Description: Regression and statistical analysis helper classes
// Created: 	2002-12-01
// Ver:		$Id: jpgraph_regstat.php 782 2006-10-08 08:09:02Z ljp $
//
// Copyright (c) Aditus Consulting. All rights reserved.
//========================================================================
*/

//------------------------------------------------------------------------
// CLASS Spline
// Create a new data array from an existing data array but with more points.
// The new points are interpolated using a cubic spline algorithm
//------------------------------------------------------------------------
class Spline {
    // 3:rd degree polynom approximation

var $xdata,$ydata;   // Data vectors
    var $y2;		 // 2:nd derivate of ydata	
    var $n=0;

function Spline($xdata,$ydata) {
	$this->y2 = array();
	$this->xdata = $xdata;
	$this->ydata = $ydata;

$n = count($ydata);
	$this->n = $n;
	if( $this->n !== count($xdata) ) {
	    JpGraphError::RaiseL(19001);
//('Spline: Number of X and Y coordinates must be the same');
	}

// Natural spline 2:derivate == 0 at endpoints
	$this->y2[0]    = 0.0;
	$this->y2[$n-1] = 0.0;
	$delta[0] = 0.0;

// Calculate 2:nd derivate
	for($i=1; $i < $n-1; ++$i) {
	    $d = ($xdata[$i+1]-$xdata[$i-1]);
	    if( $d == 0  ) {
		JpGraphError::RaiseL(19002);
//('Invalid input data for spline. Two or more consecutive input X-values are equal. Each input X-value must differ since from a mathematical point of view it must be a one-to-one mapping, i.e. each X-value must correspond to exactly one Y-value.');
	    }
	    $s = ($xdata[$i]-$xdata[$i-1])/$d;
	    $p = $s*$this->y2[$i-1]+2.0;
	    $this->y2[$i] = ($s-1.0)/$p;
	    $delta[$i] = ($ydata[$i+1]-$ydata[$i])/($xdata[$i+1]-$xdata[$i]) - 
		         ($ydata[$i]-$ydata[$i-1])/($xdata[$i]-$xdata[$i-1]);
	    $delta[$i] = (6.0*$delta[$i]/($xdata[$i+1]-$xdata[$i-1])-$s*$delta[$i-1])/$p;
	}

// Backward substitution
	for( $j=$n-2; $j >= 0; --$j ) {
	    $this->y2[$j] = $this->y2[$j]*$this->y2[$j+1] + $delta[$j];
	}
    }

// Return the two new data vectors
    function Get($num=50) {
	$n = $this->n ;
	$step = ($this->xdata[$n-1]-$this->xdata[0]) / ($num-1);
	$xnew=array();
	$ynew=array();
	$xnew[0] = $this->xdata[0];
	$ynew[0] = $this->ydata[0];
	for( $j=1; $j < $num; ++$j ) {
	    $xnew[$j] = $xnew[0]+$j*$step;
	    $ynew[$j] = $this->Interpolate($xnew[$j]);
	}
	return array($xnew,$ynew);
    }

// Return a single interpolated Y-value from an x value
    function Interpolate($xpoint) {

$max = $this->n-1;
	$min = 0;

// Binary search to find interval
	while( $max-$min > 1 ) {
	    $k = ($max+$min) / 2;
	    if( $this->xdata[$k] > $xpoint ) 
		$max=$k;
	    else 
		$min=$k;
	}

// Each interval is interpolated by a 3:degree polynom function
	$h = $this->xdata[$max]-$this->xdata[$min];

if( $h == 0  ) {
	    JpGraphError::RaiseL(19002);
//('Invalid input data for spline. Two or more consecutive input X-values are equal. Each input X-value must differ since from a mathematical point of view it must be a one-to-one mapping, i.e. each X-value must correspond to exactly one Y-value.');
	}

$a = ($this->xdata[$max]-$xpoint)/$h;
	$b = ($xpoint-$this->xdata[$min])/$h;
	return $a*$this->ydata[$min]+$b*$this->ydata[$max]+
	     (($a*$a*$a-$a)*$this->y2[$min]+($b*$b*$b-$b)*$this->y2[$max])*($h*$h)/6.0;
    }
}

//------------------------------------------------------------------------
// CLASS Bezier
// Create a new data array from a number of control points
//------------------------------------------------------------------------
class Bezier {
/**
 * @author Thomas Despoix, openXtrem company
 * @license released under QPL
 * @abstract Bezier interoplated point generation,
 * computed from control points data sets, based on Paul Bourke algorithm :
 * http://astronomy.swin.edu.au/~pbourke/curves/bezier/
 */
    var $datax = array();
    var $datay = array();
    var $n=0;
 
    function Bezier($datax, $datay, $attraction_factor = 1) {
	// Adding control point multiple time will raise their attraction power over the curve   
	$this->n = count($datax);
	if( $this->n !== count($datay) ) {
	    JpGraphError::RaiseL(19003);
//('Bezier: Number of X and Y coordinates must be the same');
	}
	$idx=0;
	foreach($datax as $datumx) {
	    for ($i = 0; $i < $attraction_factor; $i++) {
		$this->datax[$idx++] = $datumx;
	    }
	}
   	$idx=0;
	foreach($datay as $datumy) {
	    for ($i = 0; $i < $attraction_factor; $i++) {
		$this->datay[$idx++] = $datumy;
	    }
	}
	$this->n *= $attraction_factor;
    }

function Get($steps) {
	$datax = array();
	$datay = array();
	for ($i = 0; $i < $steps; $i++) {
	    list($datumx, $datumy) = $this->GetPoint((double) $i / (double) $steps);       
	    $datax[] = $datumx;
	    $datay[] = $datumy;
	}
   
	$datax[] = end($this->datax);
	$datay[] = end($this->datay);
   
	return array($datax, $datay);
    }
 
    function GetPoint($mu) {
	$n = $this->n - 1;
	$k = 0;
	$kn = 0;
	$nn = 0;
	$nkn = 0;
	$blend = 0.0;
	$newx = 0.0;
	$newy = 0.0;

$muk = 1.0;
	$munk = (double) pow(1-$mu,(double) $n);

for ($k = 0; $k <= $n; $k++) {
	    $nn = $n;
	    $kn = $k;
	    $nkn = $n - $k;
	    $blend = $muk * $munk;
	    $muk *= $mu;
	    $munk /= (1-$mu);
	    while ($nn >= 1) {
		$blend *= $nn;
		$nn--;
		if ($kn > 1) {
		    $blend /= (double) $kn;
		    $kn--;
		}
		if ($nkn > 1) {
		    $blend /= (double) $nkn;
		    $nkn--;
		}
	    }
	    $newx += $this->datax[$k] * $blend;
	    $newy += $this->datay[$k] * $blend;
	}

return array($newx, $newy);
    }
}

// EOF
?>

cmd: