-- First use Tools| Generate a reference table to make
-- a table named "ref1" with x values from 0 to 3.0, spaced 0.01 apart.

drop table if exists levytemp;
create table levytemp (
    x double,
    y double ,
   y2 double, 
   y3 double,
   y4 double,
   y5 double
);

insert into levytemp(x) 
  select item from ref1;

set @c = 0.5;
update levytemp set y = sqrt(@c / (PI() * 2.0)) * (exp(-@c/ (2.0 * x)))/(pow(x, 1.5 )) ;

set @c = 1.0;
update levytemp set y2 = sqrt(@c / (PI() * 2.0)) * (exp(-@c/ (2.0 * x)))/(pow(x, 1.5 )) ;

set @c = 2.0;
update levytemp set y3 = sqrt(@c / (PI() * 2.0)) * (exp(-@c/ (2.0 * x)))/(pow(x, 1.5 )) ;

set @c = 4.0;
update levytemp set y4 = sqrt(@c / (PI() * 2.0)) * (exp(-@c/ (2.0 * x)))/(pow(x, 1.5 )) ;

set @c = 8.0;
update levytemp set y5 = sqrt(@c / (PI() * 2.0)) * (exp(-@c/ (2.0 * x)))/(pow(x, 1.5 )) ;

plot xychart
    x, red  line, dark green  line, 
        blue  line
  with
    title "Levy Distributions"
    scale y 0 1.0 0.25
    format y decimal "#.##"

    scale y2 0 1.0 0.25
    no labels y2
    ticks y2 left

    scale x 0 3.0 1.0
    scale x2 0 3.0 1.0
    no labels x2
    ticks x2 bottom
    light gray gridlines
    legend bottom

--  Push the frame out to the edges for iPhone etc.
--    FRAME 60 41 380 174

select 
  x,
  y as "c = 0.5",
  y2 as "c = 1.0",
  y3 as "c = 2.0"
 from levytemp;