\cat geom-ex.r | pr -t -n:2 - 1:N = 40; 2:x = linspace(-1,1,N); 3:y = x; 4:z = zeros(N, N); 5: 6:for (i in 1:N) 7:{ 8: for (j in 1:N) 9: { 10: r = sqrt(x[i]^2 + y[j]^2); 11: z[i;j]=.4*(exp(-r*r)*cos(2*pi*r)); 12: } 13:} 14: 15:GM = "|/usr/local/bin/geomview -c -"; 16:fprintf(GM,... 17: "(geometry rmesh { : rmesh })\n"); 18:fprintf(GM,... 19: "(read geometry { define rmesh\n"); 20:fprintf(GM,"MESH\n"); 21:fprintf(GM,"%i %i\n", N, N); 22: 23:for (i in 1:x.n) 24:{ 25: for (j in 1:y.n) 26: { 27: fprintf(GM,"%f %f %f\n",... 28: x[i],y[j],z[i;j]); 29: } 30:} 31: 32:fprintf(GM, "})\n");