SAS 11.11

In China, 11.11 is Single’s Day.
I found this day very interesting, although 1111 looks like a solitude date on the calendar. It is not!

It has a prime factorization 1111=101*11,  and both 11 and 101 are primes. Therefore I think 1111 is a cute number and a Prime relationship day!

Some people may be like Composite Numbers, always have friends(factors). Others may like prime numbers, only solitude accompanies them. 
However, I believe all the primes number can meet up with other primes one day, and make a Composite Number! 
For those did not meet with other primes, I think maybe they are the precious Twin Prime.

data vday;
do n=1 to 6283*16;
t=n/10000;
x=16*sin(t)*sin(t)*sin(t);
y=13*cos(t)-5*cos(2*t)-2*cos(3*t)-cos(4*t);
output;
end;
run;
goptions colors=(pink);

proc gplot data=vday ;
plot y*x;
run;

The SAS Code is originated from Joe Wang .