#include <iostream>
#include <stdio.h>
#include <math.h>
using namespace std;
int main(int argc, char *argv[])
{
double w,A,x,y,Z,sigma;
x = atof(argv[1]); // Between -1:1
y = atof(argv[2]); // Between -1:1
w = atof(argv[3]); // Between 1:100
A = atof(argv[4]); // Between 0:20
sigma = atof(argv[5]); // Between 0:1
Z = A * sin((pow(x,2) + pow(y,2))*w) * cos((pow(x,2) + pow(y,2))*w);
Z = Z * exp(-(pow(x,2) + pow(y,2))/sigma);
cout << Z << ", ";
return 0;
}
I2luY2x1ZGUgPGlvc3RyZWFtPgojaW5jbHVkZSA8c3RkaW8uaD4KI2luY2x1ZGUgPG1hdGguaD4KdXNpbmcgbmFtZXNwYWNlIHN0ZDsKCmludCBtYWluKGludCBhcmdjLCBjaGFyICAqYXJndltdKQp7CgogICAgZG91YmxlIHcsQSx4LHksWixzaWdtYTsKICAgIHggPSBhdG9mKGFyZ3ZbMV0pOyAvLyBCZXR3ZWVuIC0xOjEKICAgIHkgPSBhdG9mKGFyZ3ZbMl0pOyAvLyBCZXR3ZWVuIC0xOjEKICAgIHcgPSBhdG9mKGFyZ3ZbM10pOyAvLyBCZXR3ZWVuIDE6MTAwCiAgICBBID0gYXRvZihhcmd2WzRdKTsgLy8gQmV0d2VlbiAwOjIwCiAgICBzaWdtYSA9IGF0b2YoYXJndls1XSk7IC8vIEJldHdlZW4gMDoxCgogICAgWiA9IEEgKiBzaW4oKHBvdyh4LDIpICsgcG93KHksMikpKncpICogY29zKChwb3coeCwyKSArIHBvdyh5LDIpKSp3KTsKICAgIFogPSBaICogZXhwKC0ocG93KHgsMikgKyBwb3coeSwyKSkvc2lnbWEpOwogICAgY291dCA8PCBaIDw8ICIsICI7CgogICAgcmV0dXJuIDA7Cn0K