Ezt a részt azért hagytam benne, hogy ellenőrízni lehessen a QWP helyes működését.

+-45 foknál cirkulárisan polarizált lesz, alatta elliptikus. Az előjelcsere megfordítja a forgási irányt, vagyis helyes a szimuláció.

        double dist1=0.0;
        for(int p=0;p<1000;p++)//quarterwaveplate
        {
            vec2d amp_dp,amp_ds;
            add_quarterwaveplate(&amp_ds,0.0/radian , dist1*M_PI/100.0,-45.0/radian,0.5, 0.5,0.0);

            line(200,200,200+200.0*amp_ds.x,200+200.0*amp_ds.y,0x00ff00);
            dist1+=5.0;
            
            XFlush(dpy);
            getchar();
        }