clear all;
clc;
fprintf('Choose a test function \n');
fprintf(' 1. sqrt((z+1)(1-z))cos(x/2+y/3+z/5)\n');
fprintf(' 2. cos(30xyz) if z<1/2 and sqrt(x^2+y^2+z^2)^3 if z>=1/2 \n');
fprintf(' 3. exp(-sin(2x^2+3y^2+5z^2))\n');
fprintf(' 4. exp(-sqrt(x^2+y^2+z^2))/(1+x^2+y^2+z^2)*cos(1+x^2)sin(1-y^2)exp(z)\n');
s=input('Test Function #');
fprintf('Choose a manifold \n');
fprintf(' 1. Cone\n');
fprintf(' 2. Cylinder\n');
fprintf(' 3. Torus\n');
fprintf(' 4. Sphere\n');
t=input('Manifold #');
fprintf('Choose how many points ');
N=input('#');
QMCManifold2(N,t,s)