/*
 Mehrfachreflexion am Kreis x^2+y^2=1 mit punktfoermiger Lichtquelle in (1,0)
 mit der Lieberherr-Methode
 
 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;

// 5-fache Reflexion
for(long i=0;i<anz;i++)
    {

    pkt P(cos(5*w*i),sin(5*w*i));
    pkt Q(cos(6*w*i),sin(6*w*i));
    
    pixel_strecke(L,P,pix);
    pixel_strecke(P,Q,pix);
};

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

}