Вместо генерации файла вы можете использовать специальное имя файла "++", см. help special
.Я думаю, что последний пример на демонстрационной странице gnuplot имеет ваш вариант использования.Упрощено с незначительными изменениями:
xx(u, v) = cos(v) * cos(u)
yy(u, v) = cos(v) * sin(u)
zz(u, v) = sin(v)
f(theta, phi) = sqrt(1 - (sin(theta)*sin(phi))**2) + 5 * sqrt(1 - (sin(theta)*cos(phi))**2)
set parametric
set isosamples 121, 61
set samples 121, 61
set urange [-pi:pi]
set vrange [-pi/2:pi/2]
set border 4095
set view equal xyz
set xyplane 0
splot "++" using (xx($1,$2)):(yy($1,$2)):(zz($1,$2)):(f($1,$2)) with pm3d notitle
Это результат:
Пожалуйста, дважды проверьте, соответствуют ли определения сферических координатматч.