// Sample and plot N-by-N points on a parametric surface
// Update and get geometry information
model.component(“comp1”).geom(“geom1”).run(“fin”);
GeomSequence geom1 = model.component(“comp1”).geom(“geom1”);
GeomFeature ps1 = geom1.feature(“ps1”);
double minMaxS[] = geom1.faceParamRange(faceNum);
double minS1 = minMaxS[0];
double maxS1 = minMaxS[1];
double minS2 = minMaxS[2];
double maxS2 = minMaxS[3];
double s1List[] = new double[N];
double s2List[] = new double[N];
for (int k = 0; k < N; k++) {
s1List[k] = (double) (N-1-k)/(double) (N-1)*minS1+k/(double) (N-1)*maxS1;
s2List[k] = (double) (N-1-k)/(double) (N-1)*minS2+k/(double) (N-1)*maxS2;
double s12List[][] = new double[N*N][2];
for (int i = 0; i < N; i++) {
for (int j = 0; j < N; j++) {
s12List[i+N*j][0] = s1List[i];
s12List[i+N*j][1] = s2List[j];
double[][] XY = geom1.faceX(faceNum, s12List);
String pgTag = model.result().uniquetag(“pg”);
ResultFeature pg = model.result().create(pgTag, 3);
ResultFeature plot = pg.create(“pt1”, “PointData”);
double[][] p = new double[3][N*N];
double[] color = new double[N*N];
for (int i = 0; i < N*N; i++) {
.set(“coloring”, “colortable”)
.set(“sphereradiusscale”, 1);