#!/usr/bin/gnuplot
# Author: Ika
# Date: 11/20/2014

# last changes:
# 2018-01-28 10h00: (by JoKalliauer) instead of 50 plots with one line each, over each other, changed to 1 plot with 50 lines
# 2018-01-28 10h07: (by JoKalliauer) changed 'Times' to 'Liberation Serif'

set term svg font 'Liberation Serif, 24'
set output 'envelope_astroid.svg'

# This program gives the envelope of the family of lines connecting points (s,0), (0,t) with s^2 + t^2 = 1.
g = 10;
v = 10;
# n is the number of the trajectories
n = 51;

set multiplot
set samples 1000

set xrange [-1:1]
set yrange [-1:1]
set xtics 1
set ytics 1
set xzeroaxis lt -1 lc rgb "black" lw 1
set yzeroaxis lt -1 lc rgb "black" lw 1

unset label

set key off

#t = (1.0*i)/n;
#s = sqrt(1-t*t)
# The following draws the lines. 0 will be used to cut the lines.

plot for [i=1-n:n-1] x>0 ? (x<sqrt(1-((1.0*i)/n)*((1.0*i)/n)) ? ((1.0*i)/n) - x*((1.0*i)/n)/sqrt(1-((1.0*i)/n)*((1.0*i)/n)) : 0) : (x>-sqrt(1-((1.0*i)/n)*((1.0*i)/n)) ? ((1.0*i)/n) + x*((1.0*i)/n)/sqrt(1-((1.0*i)/n)*((1.0*i)/n)) : 0) lc rgb "red" lw 1;

unset multiplot

set term x11
set output