Welcome to Fractal Forums

procedure MyFormula(var x, y, z, w: Double; PIteration3D: TPIteration3D);
begin
  ...
end;

Parameter double scale=2; { the =2 is just to say it's the default value}
Parameter double offset=2;
Constant double sqrt_1_2=sqrt(0.5);
Constant double sqrt_1_3=sqrt(1./3.);
Constant double sqrt_2_3=sqrt(2./3.);
procedure BenesiPineTree(var x,y,z,w :double; PIteration3D : TPIteration3D);
var
   start_x, start_y, start_z : double;
   ... etc.
begin
   start_x = x; start_y = y; start_z = z;
   ...etc.
end

procedure BenesiPineTree(var x,y,z,w :double; PIteration3D : TPIteration3D);
var
{added by preprocessor}
   scale, offset : double;
   sqrt_1_2, sqrt_1_3, sqrt_2_3 : double;
{/added by preprocessor}
   start_x, start_y, start_z : double;
   ... etc.
begin
{added by preprocessor}
   sqrt_1_2 = PDouble(Integer(PIteration3D^.PVar) + 0)^;
   sqrt_1_3 = PDouble(Integer(PIteration3D^.PVar) + 8)^;
   sqrt_1_2 = PDouble(Integer(PIteration3D^.PVar) + 16)^;

   offset = PDouble(Integer(PIteration3D^.PVar) - 16)^;
   scale = PDouble(Integer(PIteration3D^.PVar) - 24)^;
{/added by preprocessor}

   start_x = x; start_y = y; start_z = z;
   ...etc.
end

.Double Offset 1 = 0 {InProcName = Offset_1} {Comment = Use it to change the shape! is it clear?}

procedure BenesiPineTree(var x, y, z, w: Double; PIteration3D: TPIteration3D);
var
  // begin preprocessor
  sqrt_1_2: Double;
  sqrt_1_3: Double;
  sqrt_2_3: Double;
  Offset: Double;
  Scale: Double;
  // end preprocessor
  temp_x, temp_y, temp_z, temp: Double;  
  new_x, new_y, new_z: Double;
  x_temp, y_temp, z_temp: Double;
begin
  // begin preprocessor
  sqrt_1_2 := PDouble(Integer(PIteration3D^.PVar) + 0)^;
  sqrt_1_3 := PDouble(Integer(PIteration3D^.PVar) + 8)^;
  sqrt_2_3 := PDouble(Integer(PIteration3D^.PVar) + 16)^;
  Offset := PDouble(Integer(PIteration3D^.PVar) - 16)^;
  Scale := PDouble(Integer(PIteration3D^.PVar) - 24)^;
  // end preprocessor
  temp_x := x*sqrt_2_3-z*sqrt_1_3;
  new_z := x*sqrt_1_3 + z*sqrt_2_3; 
  new_x := (temp_x-y)*sqrt_1_2; 
  new_y := (temp_x+y)*sqrt_1_2;

  new_x := abs(new_x);
  new_y := abs(new_y);
  new_z := abs(new_z); 

  temp_x := (new_x+new_y)*sqrt_1_2;
  new_y := (-new_x+new_y)*sqrt_1_2; 
  new_x := temp_x*sqrt_2_3+new_z*sqrt_1_3;
  new_z := -temp_x*sqrt_1_3+new_z*sqrt_2_3; 

  x:=scale*new_x-offset; //scale =2 offset=2 
  y:=scale*new_y;
  z:=scale*new_z;

  x_temp:=x*x;
  y_temp:=y*y; 
  z_temp:=z*z;
  temp:=(x+x)/sqrt(y_temp+z_temp); 

  new_x:= x_temp - y_temp - z_temp;
  new_z:=temp*(y_temp-z_temp); //y and z variables are switched to add symmetry 
  new_y:=temp*2.0*y*z;

  x := new_x + PIteration3D^.J1; 
  y := new_y; 
  z := new_z;
end;

[OPTIONS]
.DEScale = 1
.SIPow = 2
.Version = 9
.Double Offset = 2
.Double Scale = 2
[CONSTANTS]
.Double sqrt_1_2 = 0.707106781186548
.Double sqrt_1_3 = 0.577350269189626
.Double sqrt_2_3 = 0.816496580927726
[SOURCE]
procedure BenesiPineTree(var x, y, z, w: Double; PIteration3D: TPIteration3D);
var
  temp_x, temp_y, temp_z, temp: Double;  
  new_x, new_y, new_z: Double;
  x_temp, y_temp, z_temp: Double;
begin
  temp_x := x*sqrt_2_3-z*sqrt_1_3;
  new_z := x*sqrt_1_3 + z*sqrt_2_3; 
  new_x := (temp_x-y)*sqrt_1_2; 
  new_y := (temp_x+y)*sqrt_1_2;

  new_x := abs(new_x);
  new_y := abs(new_y);
  new_z := abs(new_z); 

  temp_x := (new_x+new_y)*sqrt_1_2;
  new_y := (-new_x+new_y)*sqrt_1_2; 
  new_x := temp_x*sqrt_2_3+new_z*sqrt_1_3;
  new_z := -temp_x*sqrt_1_3+new_z*sqrt_2_3; 

  x:=scale*new_x-offset; //scale =2 offset=2 
  y:=scale*new_y;
  z:=scale*new_z;

  x_temp:=x*x;
  y_temp:=y*y; 
  z_temp:=z*z;
  temp:=(x+x)/sqrt(y_temp+z_temp); 

  new_x:= x_temp - y_temp - z_temp;
  new_z:=temp*(y_temp-z_temp); //y and z variables are switched to add symmetry 
  new_y:=temp*2.0*y*z;

  x := new_x + PIteration3D^.J1; 
  y := new_y; 
  z := new_z;
end;

[END]
Description:

procedure BenesiPineTree(var x, y, z, w: Double; PIteration3D: TPIteration3D);
begin
if a=1 then x:= int(x*Scale);
if b=1 then y:= int(y*Scale);
if c=1 then z:= int(z*Scale);
if a=2 then x:= frac(x*Scale);
if b=2 then y:= frac(y*Scale);
if c=2 then z:= frac(z*Scale);
if a=3 then x:= exp(x*Scale);
if b=3 then y:= exp(y*Scale);
if c=3 then z:= exp(z*Scale);
if a=4 then x:= cos(x*Scale);
if b=4 then y:= Cos(y*Scale);
if c=4 then z:= Cos(z*Scale);
if a=5 then x:= sin(x*Scale);
if b=5 then y:= sin(y*Scale);
if c=5 then z:= sin(z*Scale);
if a=6 then x:= ln(x*Scale);
if b=6 then y:= ln(y*Scale);
if c=6 then z:= ln(z*Scale);
if a=7 then x:= arctan(x*Scale);
if b=7 then y:= arctan(y*Scale);
if c=7 then z:= arctan(z*Scale);
if a=8 then x:= sqrt(x*Scale);
if b=8 then y:= sqrt(y*Scale);
if c=8 then z:= sqrt(z*Scale);
if a=9 then x:= arccos(x*Scale);
if b=9 then y:= arccos(y*Scale);
if c=9 then z:= arccos(z*Scale);
if a=10 then x:= arcsin(x*Scale);
if b=10 then y:= arcsin(y*Scale);
if c=10 then z:= arcsin(z*Scale);
if a=11 then x:= tan(x*Scale);
if b=11 then y:= tan(y*Scale);
if c=11 then z:= tan(z*Scale);
if a=12 then x:= cotan(x*Scale);
if b=12 then y:= cotan(y*Scale);
if c=12 then z:= cotan(z*Scale);
if a=13 then x:= secant(x*Scale);
if b=13 then y:= secant(y*Scale);
if c=13 then z:= secant(z*Scale);
if a=14 then x:= cosecant(x*Scale);
if b=14 then y:= cosecant(y*Scale);
if c=14 then z:= cosecant(z*Scale);
//if a=15 then x:= hypot(x*Scale);
//if b=15 then y:= hypot(y*Scale);
//if c=15 then z:= hypot(z*Scale);

if a=15 then x:= cot(x*Scale);
if b=15 then y:= cot(y*Scale);
if c=15 then z:= cot(z*Scale);
if a=16 then x:= sec(x*Scale);
if b=16 then y:= sec(y*Scale);
if c=16 then z:= sec(z*Scale);
if a=17 then x:= csc(x*Scale);
if b=17 then y:= csc(y*Scale);
if c=17 then z:= csc(z*Scale);
if a=18 then x:= cosh(x*Scale);
if b=18 then y:= cosh(y*Scale);
if c=18 then z:= cosh(z*Scale);
if a=19 then x:= sinh(x*Scale);
if b=19 then y:= sinh(y*Scale);
if c=19 then z:= sinh(z*Scale);
if a=20 then x:= tanh(x*Scale);
if b=20 then y:= tanh(y*Scale);
if c=20 then z:= tanh(z*Scale);
if a=21 then x:= coth(x*Scale);
if b=21 then y:= coth(y*Scale);
if c=21 then z:= coth(z*Scale);
if a=22 then x:= sech(x*Scale);
if b=22 then y:= sech(y*Scale);
if c=22 then z:= sech(z*Scale);
if a=23 then x:= csch(x*Scale);
if b=23 then y:= csch(y*Scale);
if c=23 then z:= csch(z*Scale);
if a=24 then x:= arccot(x*Scale);
if b=24 then y:= arccot(y*Scale);
if c=24 then z:= arccot(z*Scale);
if a=25 then x:= arcsec(x*Scale);
if b=25 then y:= arcsec(y*Scale);
if c=25 then z:= arcsec(z*Scale);
if a=26 then x:= arccsc(x*Scale);
if b=26 then y:= arccsc(y*Scale);
if c=26 then z:= arccsc(z*Scale);
if a=27 then x:= arccosh(x*Scale);
if b=27 then y:= arccosh(y*Scale);
if c=27 then z:= arccosh(z*Scale);
if a=28 then x:= arcsinh(x*Scale);
if b=28 then y:= arcsinh(y*Scale);
if c=28 then z:= arcsinh(z*Scale);
if a=29 then x:= arctanh(x*Scale);
if b=29 then y:= arctanh(y*Scale);
if c=29 then z:= arctanh(z*Scale);
if a=30 then x:= arccoth(x*Scale);
if b=30 then y:= arccoth(y*Scale);
if c=30 then z:= arccoth(z*Scale);
if a=31 then x:= arcsech(x*Scale);
if b=31 then y:= arcsech(y*Scale);
if c=31 then z:= arcsech(z*Scale);
if a=32 then x:= arccsch(x*Scale);
if b=32 then y:= arccsch(y*Scale);
if c=32 then z:= arccsch(z*Scale);

//if d=1 then let x:=(x*scale)
//if e=1 then let y:=(y*scale)
//if f=1 then let z:=(z*scale)

end;