File:FatTrifolium.svg

/*

Batch file for Maxima CAS
save as a o.mac
run maxima :
       
 maxima
 
and then : 

batch("a.mac");

*/
kill(all);
remvalue(all);
ratprint:false; /*  It doesn't change the computing, just the warnings. */
display2d:false;
numer: true;





/* 

folium curve
https://www.mathcurve.com/courbes2d.gb/foliumsimple/foliumsimple.shtml

r = a* cos^3(t)

with rotation 
https://math.stackexchange.com/questions/154743/how-to-rotate-a-polar-equation 

*/

FoliumCurve_Polar(a,rotation):= 
	/* polar(radius, ang, minang, maxang) */
	polar(a*(sin(t + rotation))^3,t,0,2*%pi)$
	
	

FoliumCurves_Polar(aMin, aMax, aStep):= block(

	[a,l, ll],
	ll:[],
	for a:aMin step aStep thru aMax do (
		l:FoliumCurve_Polar(a, 0),
		ll: cons(l,ll),
		l:FoliumCurve_Polar(a, 2*%pi/3),
		ll: cons(l,ll) ,
		l:FoliumCurve_Polar(a, 4*%pi/3),
		ll: cons(l,ll)  
		 
	),
	
	return (ll)


)$




ll:FoliumCurves_Polar(1, 20, 0.5)$


load(draw);

/* strings */
path:"~/c/julia/parabolic/hawaian/3/"$ /*  if empty then file is in a home dir , path should end with "/" */ 
FileName:  "FatTrifolium"$


draw2d(
	terminal      = svg,
	file_name = sconcat(path, FileName,""),
	title = "Fat trifolium curves ",
	dimensions = [1000, 1000], 
	nticks=300,
	ll)$