function [v] = r(x)
v = exp(-x^2);