Re: [Pdl-general] rotate matrix
here is the big orbitals and all there combinations, I also would like to
show you that $n is important( I noticed you took it out)
in that it can give you all the colors...
$Best,
-Mark Baker
use PDL;
use PDL::Complex;
use PDL::Graphics::TriD;
use PDL::NiceSlice;
use PDL::Math; keeptwiddling3d;
$PDL::BIGPDL=1;
## remember this can crash the computer if you dont have enough RAM and
swap in linux or pagefile in windows
for (;;){
$c = (4) ; ## speed
for $phase(-360..360) { ## 360
$phases = $phase*0.001;
$frequency = ($c**-1);
$n= (6.28*$c*$frequency+($phases));
$photons = 4178; ### [Density] 2444 4178
$t=rvals(exp(zeros($photons)));
$cz=-1**$t*$c; # -1**$t*$c
$cy=-1**sin($t*$c); #-1**$c*sin($t*$c)
$cx=-1**$c*sin(rvals($t))*$c; #-1**$c*bessj0(rvals($t,{Centre=>[0]}))*$c;
$w=$cz-$cy-$cx;
$g=sin($w); #sin
$r=cos($cy+$c+$cz); #cos
$b=cos($w); #cos
$i=($cz-$cx-$cy);
$q=$c*($i*$n);
### altitude longitude latitude
$x = $b*sin($q); ### cos change the spin of the electron here
$y = $r*cos($q); ### sin
$z = $g*sin($q); ### cos
$rad=sum(sqrt(($x**2)+($y**2)+($z**2)));
$theta= sin($y*$x*$z)**-1;
$phi=cos($z*$rad)->rvals;
$x = $rad*cos($theta)*sin($phi); ### cos change the spin of the electron
here
$y = $rad*sin($theta)*sin($phi); ### sin
$z = $rad*cos($phi); ### cos
$xx = $b*sin($b); ### cos change the spin of the electron here
$yy = $r*($q); ### sin
$zz = $g*cos($g); ### cos
$dd++;
## torus
#$xxx = $xx*sin($yy)*cos($zz); ### cos change the spin of the electron
here
#$yyy = $xx*cos($yy)*cos($zz); ### sin
#$zzz = $xx*sin($zz); ### cos
# Draw a line
### d orbital i think
$ = $xx*sin($yy)*cos($zz)+$xx; ### cos change the spin of the
electron here
$ = $xx*cos($yy)*cos($zz)+$xx; ### sin
$ = $xx*sin($zz); ### cos
# Draw a line
p orbital
$xxx = $xx*cos($yy)*sin($zz); ### cos change the spin of the electron
here
$yyy = $xx*sin($yy)*sin($zz); ### sin
$zzz = $xx*cos($zz);
$vv=sequence(4178);
### wavefunction #exp(-sqrt(($x**2)+($y**2)+($z**2)))
#$xxx=rotate($xxx,$phase);
#$yyy=rotate($yyy,$phase);
#$zzz=rotate($zzz,$phase);
#$xxx=$xxx*$vv;
#$yyy=$yyy*$vv;
#$zzz=$zzz*$vv;
#$=$*$vv;
#$=$*$vv;
#$=$*$vv;
$matrixa = ones(100);
$matrixb = ones(100);
$matrixc = ones(100);
$colora = ones(100);
$colorb = ones(100);
$colorc = ones(100);
#
# changing the colors can show the second electron in orbit
$matrixa->range(5000,[4178]) .=(($))-.2 ;
$matrixb->range(5000,[4178]) .=((-$))-.2 ;
$matrixc->range(5000,[4178]) .=(($))-.2 ;
$colora->range(5000,[4178]) .= $r+$n ;
$colorb->range(5000,[4178]) .= $g+$n ;
$colorc->range(5000,[4178]) .= $b+$n ;
$matrixa->range(1,[4178]) .=((-$))-.2 ;
$matrixb->range(1,[4178]) .=(($))-.2 ;
$matrixc->range(1,[4178]) .=((-$))-.2 ;
$colora->range(1,[4178]) .= $r+$n ;
$colorb->range(1,[4178]) .= $g+$n ;
$colorc->range(1,[4178]) .= $b+$n ;
$matrixa->range(15000,[4178]) .=((-$))-.2 ;
$matrixb->range(15000,[4178]) .=(($))-.2 ;
$matrixc->range(15000,[4178]) .=(($))-.2 ;
$colora->range(15000,[4178]) .= $r-$n ;
$colorb->range(15000,[4178]) .= $g-$n ;
$colorc->range(15000,[4178]) .= $b+$n ;
$matrixa->range(2,[4178]) .=(($))-.2 ;
$matrixb->range(2,[4178]) .=((-$))-.2 ;
$matrixc->range(2,[4178]) .=((-$))-.2 ;
$colora->range(2,[4178]) .= $r-$n ;
$colorb->range(2,[4178]) .= $g-$n ;
$colorc->range(2,[4178]) .= $b+$n ;
$matrixa->range(25000,[4178]) .=((-$))-.2 ;
$matrixb->range(25000,[4178]) .=(($))-.2 ;
$matrixc->range(25000,[4178]) .=(($))-.2 ;
$colora->range(25000,[4178]) .= $r+$n ;
$colorb->range(25000,[4178]) .= $g-$n ;
$colorc->range(25000,[4178]) .= $b-$n ;
$matrixa->range(3,[4178]) .=(($))-.2 ;
$matrixb->range(3,[4178]) .=((-$))-.2 ;
$matrixc->range(3,[4178]) .=((-$))-.2 ;
$colora->range(3,[4178]) .= $r+$n ;
$colorb->range(3,[4178]) .= $g-$n ;
$colorc->range(3,[4178]) .= $b-$n ;
$matrixa->range(35000,[4178]) .=(($))-.2 ;
$matrixb->range(35000,[4178]) .=((-$))-.2 ;
$matrixc->range(35000,[4178]) .=((-$))-.2 ;
$colora->range(35000,[4178]) .= $r-$n ;
$colorb->range(35000,[4178]) .= $g+$n ;
$colorc->range(35000,[4178]) .= $b+$n ;
$matrixa->range(4,[4178]) .=((-$))-.2 ;
$matrixb->range(4,[4178]) .=(($))-.2 ;
$matrixc->range(4,[4178]) .=(($))-.2 ;
$colora->range(4,[4178]) .= $r-$n ;
$colorb->range(4,[4178]) .= $g+$n ;
$colorc->range(4,[4178]) .= $b+$n ;
$matrixa->range(45000,[4178]) .=(($))-.2 ;
$matrixb->range(45000,[4178]) .=(($))-.2 ;
$matrixc->range(45000,[4178]) .=(($))-.2 ;
$colora->range(45000,[4178]) .= $r+$n ;
$colorb->range(45000,[4178]) .= $g+$n ;
Re: [Pdl-general] rotate matrix
i found the torus matrix here
$xxx = $xx*sin($yy)*cos($zz);
$yyy = $xx*cos($yy)*cos($zz);
$zzz = $xx*sin($zz);
and if you add the $xx to $xxx and $yyy you get bigger orbs
like
$xxx = $xx*sin($yy)*cos($zz)+$xx;
$yyy = $xx*cos($yy)*cos($zz)+$xx;
$zzz = $xx*sin($zz);
would that be all the orbitals i need ???
Best,
-Mark Baker
On Sun, Oct 15, 2023 at 5:46 PM Mark NanoNebulas
wrote:
> the s orbital is x y z the p orbital is the d orbitals are
> xxx yyy zzz
> is that right and can they be used to describe all atoms ???
>
> -Mark Baker
>
> On Sun, Oct 15, 2023 at 5:26 PM Mark NanoNebulas
> wrote:
>
>> have you used the range function ??? here i have all orbitals that can
>> be without using hold3d();
>>
>> could you put the Rotxz into 4178 dimensions I dont understand how your
>> doing it in 3
>> here are all the orbitals i got 14 of them total in your opinion can i
>> use these to describe all the atoms
>>
>> check out the range function it lets you use a lot of piddles in one
>> matrix...
>>
>> Best,
>> -Mark
>> #
>>
>> use PDL;
>> use PDL::Complex;
>> use PDL::Graphics::TriD;
>> use PDL::NiceSlice;
>> use PDL::Math; keeptwiddling3d;
>> $PDL::BIGPDL=1;
>> ## remember this can crash the computer if you dont have enough RAM and
>> swap in linux or pagefile in windows
>>
>>
>>
>> for (;;){
>> $c = (4000) ; ## speed
>> for $phase(-360..360) { ## 360
>> $phases = $phase*0.001;
>> $frequency = ($c**-1);
>> $n= (6.28*$c*$frequency+($phases));
>>
>> $photons = 4178; ### [Density] 2444 4178
>>
>>
>> $t=rvals(exp(zeros($photons)));
>> $cz=-1**$t*$c; # -1**$t*$c
>> $cy=-1**sin($t*$c); #-1**$c*sin($t*$c)
>> $cx=-1**$c*sin(rvals($t))*$c; #-1**$c*bessj0(rvals($t,{Centre=>[0]}))*$c;
>>
>>
>> $w=$cz-$cy-$cx;
>> $g=sin($w); #sin
>> $r=cos($cy+$c+$cz); #cos
>> $b=cos($w); #cos
>> $i=($cz-$cx-$cy);
>> $q=$c*($i*$n);
>>
>> ### altitude longitude latitude
>> $x = $b*sin($q); ### cos change the spin of the electron here
>> $y = $r*cos($q); ### sin
>> $z = $g*sin($q); ### cos
>>
>>
>> $xx = $b*sin($b); ### cos change the spin of the electron here
>> $yy = $r*($q); ### sin
>> $zz = $g*cos($g); ### cos
>> $dd++;
>> $xxx = $xx*cos($yy)*sin($zz)+sin($xx*(1/sqrt 2)); ### cos change the
>> spin of the electron here
>> $yyy = $xx*sin($yy)*sin($zz)+sin($xx*(1/sqrt 2)); ### sin
>> $zzz = $xx*cos($zz); ### cos
>> # Draw a line
>>
>> $ = $xx*cos($yy)*sin($zz); ### cos change the spin of the electron
>> here
>> $ = $xx*sin($yy)*sin($zz); ### sin
>> $ = $xx*cos($zz);
>>
>>
>>
>> $vv=sequence(4178);
>> ### wavefunction #exp(-sqrt(($x**2)+($y**2)+($z**2)))
>>
>> #$xxx=rotate($xxx,$phase);
>> #$yyy=rotate($yyy,$phase);
>> #$zzz=rotate($zzz,$phase);
>>
>> #$xxx=$xxx*$vv;
>> #$yyy=$yyy*$vv;
>> #$zzz=$zzz*$vv;
>>
>> #$=$*$vv;
>> #$=$*$vv;
>> #$=$*$vv;
>>
>>
>>
>>
>> $matrixa = ones(100);
>> $matrixb = ones(100);
>> $matrixc = ones(100);
>> $colora = ones(100);
>> $colorb = ones(100);
>> $colorc = ones(100);
>>
>> $matrixa->range(5000,[4178]) .=0.8*(($xxx))-.2 ;
>> $matrixb->range(5000,[4178]) .=0.8*(($yyy))-.2 ;
>> $matrixc->range(5000,[4178]) .=0.8*(($zzz))-.2 ;
>> $colora->range(5000,[4178]) .= $r ;
>> $colorb->range(5000,[4178]) .= $g ;
>> $colorc->range(5000,[4178]) .= $b ;
>>
>> $matrixa->range(1,[4178]) .=0.8*((-$xxx))-.2 ;
>> $matrixb->range(1,[4178]) .=0.8*((-$yyy))-.2 ;
>> $matrixc->range(1,[4178]) .=0.8*((-$zzz))-.2 ;
>> $colora->range(1,[4178]) .= $r ;
>> $colorb->range(1,[4178]) .= $g ;
>> $colorc->range(1,[4178]) .= $b ;
>>
>>$matrixa->range(15000,[4178]) .=0.8*((-$xxx))-.2 ;
>> $matrixb->range(15000,[4178]) .=0.8*((-$yyy))-.2 ;
>> $matrixc->range(15000,[4178]) .=0.8*(($zzz))-.2 ;
>> $colora->range(15000,[4178]) .= $r ;
>> $colorb->range(15000,[4178]) .= $g ;
>> $colorc->range(15000,[4178]) .= $b ;
>>
>> $matrixa->range(2,[4178]) .=0.8*(($xxx))-.2 ;
>> $matrixb->range(2,[4178]) .=0.8*((-$yyy))-.2 ;
>> $matrixc->range(2,[4178]) .=0.8*(($zzz))-.2 ;
>> $colora->range(2,[4178]) .= $r ;
>> $colorb->range(2,[4178]) .= $g ;
>> $colorc->range(2,[4178]) .= $b ;
>>
>>$matrixa->range(25000,[4178]) .=0.8*((-$xxx))-.2 ;
>> $matrixb->range(25000,[4178]) .=0.8*(($yyy))-.2 ;
>> $matrixc->range(25000,[4178]) .=0.8*((-$zzz))-.2 ;
>> $colora->range(25000,[4178]) .= $r ;
>> $colorb->range(25000,[4178]) .= $g ;
>> $colorc->range(25000,[4178]) .= $b ;
>>
>>$matrixa->range(3,[4178]) .=0.8*(($xxx))-.2 ;
>> $matrixb->range(3,[4178]) .=0.8*((-$yyy))-.2 ;
>> $matrixc->range(3,[4178]) .=0.8*((-$zzz))-.2 ;
>> $colora->range(3,[4178]) .= $r ;
>> $colorb->range(3,[4178]) .= $g ;
>> $colorc->range(3,[4178]) .= $b ;
>>
>> $matrixa->range(35000,[4178]) .=0.8*(($xxx))-.2 ;
>> $matrixb->range(35000,[4178]) .=0.8*(($yyy))-.2 ;
>>
Re: [Pdl-general] rotate matrix
the s orbital is x y z the p orbital is the d orbitals are
xxx yyy zzz
is that right and can they be used to describe all atoms ???
-Mark Baker
On Sun, Oct 15, 2023 at 5:26 PM Mark NanoNebulas
wrote:
> have you used the range function ??? here i have all orbitals that can be
> without using hold3d();
>
> could you put the Rotxz into 4178 dimensions I dont understand how your
> doing it in 3
> here are all the orbitals i got 14 of them total in your opinion can i use
> these to describe all the atoms
>
> check out the range function it lets you use a lot of piddles in one
> matrix...
>
> Best,
> -Mark
> #
>
> use PDL;
> use PDL::Complex;
> use PDL::Graphics::TriD;
> use PDL::NiceSlice;
> use PDL::Math; keeptwiddling3d;
> $PDL::BIGPDL=1;
> ## remember this can crash the computer if you dont have enough RAM and
> swap in linux or pagefile in windows
>
>
>
> for (;;){
> $c = (4000) ; ## speed
> for $phase(-360..360) { ## 360
> $phases = $phase*0.001;
> $frequency = ($c**-1);
> $n= (6.28*$c*$frequency+($phases));
>
> $photons = 4178; ### [Density] 2444 4178
>
>
> $t=rvals(exp(zeros($photons)));
> $cz=-1**$t*$c; # -1**$t*$c
> $cy=-1**sin($t*$c); #-1**$c*sin($t*$c)
> $cx=-1**$c*sin(rvals($t))*$c; #-1**$c*bessj0(rvals($t,{Centre=>[0]}))*$c;
>
>
> $w=$cz-$cy-$cx;
> $g=sin($w); #sin
> $r=cos($cy+$c+$cz); #cos
> $b=cos($w); #cos
> $i=($cz-$cx-$cy);
> $q=$c*($i*$n);
>
> ### altitude longitude latitude
> $x = $b*sin($q); ### cos change the spin of the electron here
> $y = $r*cos($q); ### sin
> $z = $g*sin($q); ### cos
>
>
> $xx = $b*sin($b); ### cos change the spin of the electron here
> $yy = $r*($q); ### sin
> $zz = $g*cos($g); ### cos
> $dd++;
> $xxx = $xx*cos($yy)*sin($zz)+sin($xx*(1/sqrt 2)); ### cos change the
> spin of the electron here
> $yyy = $xx*sin($yy)*sin($zz)+sin($xx*(1/sqrt 2)); ### sin
> $zzz = $xx*cos($zz); ### cos
> # Draw a line
>
> $ = $xx*cos($yy)*sin($zz); ### cos change the spin of the electron
> here
> $ = $xx*sin($yy)*sin($zz); ### sin
> $ = $xx*cos($zz);
>
>
>
> $vv=sequence(4178);
> ### wavefunction #exp(-sqrt(($x**2)+($y**2)+($z**2)))
>
> #$xxx=rotate($xxx,$phase);
> #$yyy=rotate($yyy,$phase);
> #$zzz=rotate($zzz,$phase);
>
> #$xxx=$xxx*$vv;
> #$yyy=$yyy*$vv;
> #$zzz=$zzz*$vv;
>
> #$=$*$vv;
> #$=$*$vv;
> #$=$*$vv;
>
>
>
>
> $matrixa = ones(100);
> $matrixb = ones(100);
> $matrixc = ones(100);
> $colora = ones(100);
> $colorb = ones(100);
> $colorc = ones(100);
>
> $matrixa->range(5000,[4178]) .=0.8*(($xxx))-.2 ;
> $matrixb->range(5000,[4178]) .=0.8*(($yyy))-.2 ;
> $matrixc->range(5000,[4178]) .=0.8*(($zzz))-.2 ;
> $colora->range(5000,[4178]) .= $r ;
> $colorb->range(5000,[4178]) .= $g ;
> $colorc->range(5000,[4178]) .= $b ;
>
> $matrixa->range(1,[4178]) .=0.8*((-$xxx))-.2 ;
> $matrixb->range(1,[4178]) .=0.8*((-$yyy))-.2 ;
> $matrixc->range(1,[4178]) .=0.8*((-$zzz))-.2 ;
> $colora->range(1,[4178]) .= $r ;
> $colorb->range(1,[4178]) .= $g ;
> $colorc->range(1,[4178]) .= $b ;
>
>$matrixa->range(15000,[4178]) .=0.8*((-$xxx))-.2 ;
> $matrixb->range(15000,[4178]) .=0.8*((-$yyy))-.2 ;
> $matrixc->range(15000,[4178]) .=0.8*(($zzz))-.2 ;
> $colora->range(15000,[4178]) .= $r ;
> $colorb->range(15000,[4178]) .= $g ;
> $colorc->range(15000,[4178]) .= $b ;
>
> $matrixa->range(2,[4178]) .=0.8*(($xxx))-.2 ;
> $matrixb->range(2,[4178]) .=0.8*((-$yyy))-.2 ;
> $matrixc->range(2,[4178]) .=0.8*(($zzz))-.2 ;
> $colora->range(2,[4178]) .= $r ;
> $colorb->range(2,[4178]) .= $g ;
> $colorc->range(2,[4178]) .= $b ;
>
>$matrixa->range(25000,[4178]) .=0.8*((-$xxx))-.2 ;
> $matrixb->range(25000,[4178]) .=0.8*(($yyy))-.2 ;
> $matrixc->range(25000,[4178]) .=0.8*((-$zzz))-.2 ;
> $colora->range(25000,[4178]) .= $r ;
> $colorb->range(25000,[4178]) .= $g ;
> $colorc->range(25000,[4178]) .= $b ;
>
>$matrixa->range(3,[4178]) .=0.8*(($xxx))-.2 ;
> $matrixb->range(3,[4178]) .=0.8*((-$yyy))-.2 ;
> $matrixc->range(3,[4178]) .=0.8*((-$zzz))-.2 ;
> $colora->range(3,[4178]) .= $r ;
> $colorb->range(3,[4178]) .= $g ;
> $colorc->range(3,[4178]) .= $b ;
>
> $matrixa->range(35000,[4178]) .=0.8*(($xxx))-.2 ;
> $matrixb->range(35000,[4178]) .=0.8*(($yyy))-.2 ;
> $matrixc->range(35000,[4178]) .=0.8*((-$zzz))-.2 ;
> $colora->range(35000,[4178]) .= $r ;
> $colorb->range(35000,[4178]) .= $g ;
> $colorc->range(35000,[4178]) .= $b ;
>
> $matrixa->range(4,[4178]) .=0.8*((-$xxx))-.2 ;
> $matrixb->range(4,[4178]) .=0.8*(($yyy))-.2 ;
> $matrixc->range(4,[4178]) .=0.8*(($zzz))-.2 ;
> $colora->range(4,[4178]) .= $r ;
> $colorb->range(4,[4178]) .= $g ;
> $colorc->range(4,[4178]) .= $b ;
>
> $matrixa->range(45000,[4178]) .=(($))-.2 ;
> $matrixb->range(45000,[4178]) .=(($))-.2 ;
>
Re: [Pdl-general] rotate matrix
have you used the range function ??? here i have all orbitals that can be
without using hold3d();
could you put the Rotxz into 4178 dimensions I dont understand how your
doing it in 3
here are all the orbitals i got 14 of them total in your opinion can i use
these to describe all the atoms
check out the range function it lets you use a lot of piddles in one
matrix...
Best,
-Mark
#
use PDL;
use PDL::Complex;
use PDL::Graphics::TriD;
use PDL::NiceSlice;
use PDL::Math; keeptwiddling3d;
$PDL::BIGPDL=1;
## remember this can crash the computer if you dont have enough RAM and
swap in linux or pagefile in windows
for (;;){
$c = (4000) ; ## speed
for $phase(-360..360) { ## 360
$phases = $phase*0.001;
$frequency = ($c**-1);
$n= (6.28*$c*$frequency+($phases));
$photons = 4178; ### [Density] 2444 4178
$t=rvals(exp(zeros($photons)));
$cz=-1**$t*$c; # -1**$t*$c
$cy=-1**sin($t*$c); #-1**$c*sin($t*$c)
$cx=-1**$c*sin(rvals($t))*$c; #-1**$c*bessj0(rvals($t,{Centre=>[0]}))*$c;
$w=$cz-$cy-$cx;
$g=sin($w); #sin
$r=cos($cy+$c+$cz); #cos
$b=cos($w); #cos
$i=($cz-$cx-$cy);
$q=$c*($i*$n);
### altitude longitude latitude
$x = $b*sin($q); ### cos change the spin of the electron here
$y = $r*cos($q); ### sin
$z = $g*sin($q); ### cos
$xx = $b*sin($b); ### cos change the spin of the electron here
$yy = $r*($q); ### sin
$zz = $g*cos($g); ### cos
$dd++;
$xxx = $xx*cos($yy)*sin($zz)+sin($xx*(1/sqrt 2)); ### cos change the
spin of the electron here
$yyy = $xx*sin($yy)*sin($zz)+sin($xx*(1/sqrt 2)); ### sin
$zzz = $xx*cos($zz); ### cos
# Draw a line
$ = $xx*cos($yy)*sin($zz); ### cos change the spin of the electron
here
$ = $xx*sin($yy)*sin($zz); ### sin
$ = $xx*cos($zz);
$vv=sequence(4178);
### wavefunction #exp(-sqrt(($x**2)+($y**2)+($z**2)))
#$xxx=rotate($xxx,$phase);
#$yyy=rotate($yyy,$phase);
#$zzz=rotate($zzz,$phase);
#$xxx=$xxx*$vv;
#$yyy=$yyy*$vv;
#$zzz=$zzz*$vv;
#$=$*$vv;
#$=$*$vv;
#$=$*$vv;
$matrixa = ones(100);
$matrixb = ones(100);
$matrixc = ones(100);
$colora = ones(100);
$colorb = ones(100);
$colorc = ones(100);
$matrixa->range(5000,[4178]) .=0.8*(($xxx))-.2 ;
$matrixb->range(5000,[4178]) .=0.8*(($yyy))-.2 ;
$matrixc->range(5000,[4178]) .=0.8*(($zzz))-.2 ;
$colora->range(5000,[4178]) .= $r ;
$colorb->range(5000,[4178]) .= $g ;
$colorc->range(5000,[4178]) .= $b ;
$matrixa->range(1,[4178]) .=0.8*((-$xxx))-.2 ;
$matrixb->range(1,[4178]) .=0.8*((-$yyy))-.2 ;
$matrixc->range(1,[4178]) .=0.8*((-$zzz))-.2 ;
$colora->range(1,[4178]) .= $r ;
$colorb->range(1,[4178]) .= $g ;
$colorc->range(1,[4178]) .= $b ;
$matrixa->range(15000,[4178]) .=0.8*((-$xxx))-.2 ;
$matrixb->range(15000,[4178]) .=0.8*((-$yyy))-.2 ;
$matrixc->range(15000,[4178]) .=0.8*(($zzz))-.2 ;
$colora->range(15000,[4178]) .= $r ;
$colorb->range(15000,[4178]) .= $g ;
$colorc->range(15000,[4178]) .= $b ;
$matrixa->range(2,[4178]) .=0.8*(($xxx))-.2 ;
$matrixb->range(2,[4178]) .=0.8*((-$yyy))-.2 ;
$matrixc->range(2,[4178]) .=0.8*(($zzz))-.2 ;
$colora->range(2,[4178]) .= $r ;
$colorb->range(2,[4178]) .= $g ;
$colorc->range(2,[4178]) .= $b ;
$matrixa->range(25000,[4178]) .=0.8*((-$xxx))-.2 ;
$matrixb->range(25000,[4178]) .=0.8*(($yyy))-.2 ;
$matrixc->range(25000,[4178]) .=0.8*((-$zzz))-.2 ;
$colora->range(25000,[4178]) .= $r ;
$colorb->range(25000,[4178]) .= $g ;
$colorc->range(25000,[4178]) .= $b ;
$matrixa->range(3,[4178]) .=0.8*(($xxx))-.2 ;
$matrixb->range(3,[4178]) .=0.8*((-$yyy))-.2 ;
$matrixc->range(3,[4178]) .=0.8*((-$zzz))-.2 ;
$colora->range(3,[4178]) .= $r ;
$colorb->range(3,[4178]) .= $g ;
$colorc->range(3,[4178]) .= $b ;
$matrixa->range(35000,[4178]) .=0.8*(($xxx))-.2 ;
$matrixb->range(35000,[4178]) .=0.8*(($yyy))-.2 ;
$matrixc->range(35000,[4178]) .=0.8*((-$zzz))-.2 ;
$colora->range(35000,[4178]) .= $r ;
$colorb->range(35000,[4178]) .= $g ;
$colorc->range(35000,[4178]) .= $b ;
$matrixa->range(4,[4178]) .=0.8*((-$xxx))-.2 ;
$matrixb->range(4,[4178]) .=0.8*(($yyy))-.2 ;
$matrixc->range(4,[4178]) .=0.8*(($zzz))-.2 ;
$colora->range(4,[4178]) .= $r ;
$colorb->range(4,[4178]) .= $g ;
$colorc->range(4,[4178]) .= $b ;
$matrixa->range(45000,[4178]) .=(($))-.2 ;
$matrixb->range(45000,[4178]) .=(($))-.2 ;
$matrixc->range(45000,[4178]) .=(($))-.2 ;
$colora->range(45000,[4178]) .= $r ;
$colorb->range(45000,[4178]) .= $g ;
$colorc->range(45000,[4178]) .= $b ;
$matrixa->range(5,[4178]) .=((-$))-.2 ;
$matrixb->range(5,[4178]) .=((-$))-.2 ;
$matrixc->range(5,[4178]) .=((-$))-.2 ;
$colora->range(5,[4178]) .= $r ;
$colorb->range(5,[4178]) .= $g ;
$colorc->range(5,[4178]) .= $b ;
$matrixa->range(55000,[4178]) .=(($))-.2 ;
