/*
 Herzkurve als Kaustik mit der Lieberherr-Methode
 Kreis: x^2+y^2=1; Lichtquelle (1,0)
 
 Autor
 Daniel Baumgartner
 Kantonsschule Baden
 
 www.mathematicus.ch
 eMail <vorname>_<nachname>@bluewin.ch
*/ 

#include <iostream>
#include <cmath>
using namespace std;

#include "gsm.h"

pkt L(1,0);//Lichtquelle
 
main()
{
// Koordinatenbereich xmin,xmax,ymin,ymax
pixel_koord(-1.2,1.2,-1.2,1.2);

pix_X=1024;//Anzahl Pixel in x-Richtung

// Die Matrix festlegen
int m=(int) round(deltay/deltax*pix_X);
intmat pix(pix_X,m);

    
long anz=1024;//Anzahl Lichtstrahlen
double w=2*pi/anz;

for(long i=0;i<anz;i++)
    {
    
    pkt P(cos(w*i),sin(w*i));
    pkt Q(cos(2*w*i),sin(2*w*i));
    
    pixel_strecke(L,P,pix);
    pixel_strecke(P,Q,pix);
};

/* Die Graustufenmatrix erzeugen
 0=Mathematica-Syntax
 1=Octave-/Matlab-Syntax
*/
pixelbild(pix,"kardioide00",0);

}