Файл: BMonSphere.jpg

Это изображение было загружено в формате JPEG , хотя оно не является фотографическим . Эту информацию можно было бы более эффективно и точно хранить в формате PNG или SVG . Если возможно, загрузите версию этого изображения в формате PNG или SVG без артефактов сжатия , полученную из источника, отличного от JPEG (или с удаленными существующими артефактами ). После этого пометьте версию JPEG {{Superseded|NewImage.ext}}и удалите этот тег. Этот тег не следует применять к фотографиям или отсканированным изображениям. Для получения дополнительной информации см. {{BadJPEG}} .

 Это изображение было создано с помощью Blender .

	Этот файл находится под лицензией Creative Commons Attribution-Share Alike 2.5 Generic .
	Ты свободен: делиться - копировать, распространять и передавать произведение ремикс - адаптировать произведение При следующих условиях: атрибуция - вы должны указать соответствующий источник, предоставить ссылку на лицензию и указать, были ли внесены изменения. Вы можете сделать это любым разумным способом, но не любым способом, который предполагает, что лицензиар одобряет вас или ваше использование. делиться одинаково - если вы ремикшируете, трансформируете или опираетесь на материал, вы должны распространять свои материалы по той же или совместимой лицензии, что и оригинал. https://creativecommons.org/licenses/by-sa/2.5CC BY-SA 2.5 Creative Commons Attribution-Share Alike 2.5 правдаправда

код

Возможно, вы берете исходный код со страницы «редактирования» без викиформинга.

GNU R

Это создает пути и сохраняет их в текстовые файлы, которые может прочитать блендер. Существуют также пути для БМ на торе.

# вычислить броуновское движение на сфере; вывод - это список# состоящий из:# Z ... BM на сфере# Y ... тангенциальный BM, см. Price & Williams# b ... независимый 1D BM (см. Прайс и Вильямс)# B ... создание 3D BM# n ... количество временных шагов дискретизации# T ... вышеупомянутые процессы даны на однородной сетке размера# n на [0, T]эйлер = функция (x0, T, n) { # инициализировать объекты dt = T / (n-1); дБ = матрица (rep (0,3 * (n-1)), ncol = 3, nrow = n-1); дБ [, 1] = rnorm (n-1, 0, sqrt (dt)); дБ [, 2] = rnorm (n-1, 0, sqrt (dt)); дБ [, 3] = rnorm (n-1, 0, sqrt (dt)); Z = матрица (rep (0,3 * n), ncol = 3, nrow = n); dZ = матрица (rep (0,3 * (n-1)), ncol = 3, nrow = n-1); Y = матрица (rep (0,3 * n), ncol = 3, nrow = n); B = матрица (rep (0,3 * n), ncol = 3, nrow = n); b = rep (0, n); Z [1,] = x0; # делаем вычисления for (k in 2: n) { B [k,] = B [k-1,] + дБ [k-1,]; dZ [k-1,] = cross (Z [k-1,], дБ [k-1,]) - Z [k-1,] * dt; Z [k,] = Z [k-1,] + dZ [k-1,]; Y [k,] = Y [k-1,] - крест (Z [k-1,], dZ [k-1,]); b [k] = b [k-1] + точка (Z [k-1,], дБ [k-1,]); } return (список (Z = Z, Y = Y, b = b, B = B, n = n, T = T));}# записываем вывод euler в csv-файлыeuler.write = function (bms, files = c ("Z.csv", "Y.csv", "b.csv", "B.csv"), steps = bms $ n) { bigsteps = round (seq (1, bms $ n, length = шаги)) write.table (bms $ Z [bigsteps,], file = files [1], col.names = F, row.names = F, sep = ",", dec = "."); write.table (bms $ Y [bigsteps,], file = files [2], col.names = F, row.names = F, sep = ",", dec = "."); write.table (bms $ b [bigsteps], file = files [3], col.names = F, row.names = F, sep = ",", dec = "."); write.table (bms $ B [bigsteps,], file = files [4], col.names = F, row.names = F, sep = ",", dec = ".");}# вычисляем броуновское движение на трехмерном торе с внешним# радиус R и внутренний радиус reulerTorus = function (x0, r, R, t, n) { # инициализировать объекты dt = t / (n-1); дБ = матрица (rep (0,3 * (n-1)), ncol = 3, nrow = n-1); дБ [, 1] = rnorm (n-1, 0, sqrt (dt)); дБ [, 2] = rnorm (n-1, 0, sqrt (dt)); дБ [, 3] = rnorm (n-1, 0, sqrt (dt)); Z = матрица (rep (0,3 * n), ncol = 3, nrow = n); B = матрица (rep (0,3 * n), ncol = 3, nrow = n); dZ = матрица (rep (0,3 * (n-1)), ncol = 3, nrow = n-1); Z [1,] = x0; nT = rep (0,3); # делаем вычисления for (k in 2: n) { B [k,] = B [k-1,] + дБ [k-1,]; nT = nTorus (Z [k-1,], r, R); dZ [k-1,] = cross (nT, дБ [k-1,]) + HTorus (Z [k-1,], r, R) * nT * dt; Z [k,] = Z [k-1,] + dZ [k-1,]; } return (list (Z = Z, B = B, n = n, t = t));}# записываем вывод euler в csv-файлыtorus.write = function (bmt, files = c ("tZ.csv", "tB.csv"), steps = bmt $ n) { bigsteps = round (seq (1, bmt $ n, length = шаги)) write.table (bmt $ Z [bigsteps,], file = files [1], col.names = F, row.names = F, sep = ",", dec = "."); write.table (bmt $ B [bigsteps,], file = files [2], col.names = F, row.names = F, sep = ",", dec = ".");}# "определяющая" функция тораfTorus = function (x, r, R) { return ((x [1] ^ 2 + x [2] ^ 2 + x [3] ^ 2 + R ^ 2-r ^ 2) ^ 2 - 4 * R ^ 2 * (x [1] ^ 2 + x [2] ^ 2));}# вектор нормали трехмерного тора с внешним радиусом R и внутренним радиусом rnTorus = функция (x, r, R) { c1 = x [1] * (x [1] ^ 2 + x [2] ^ 2 + x [3] ^ 2-R ^ 2-r ^ 2) / (3 * x [1] ^ 4 * x [ 2] ^ 2 + 3 * x [3] ^ 4 * x [2] ^ 2 + 3 * x [3] ^ 4 * x [1] ^ 2 + 6 * x [3] ^ 2 * x [1] ^ 2 * x [2] ^ 2 + 3 * x [1] ^ 2 * x [2] ^ 4 + 3 * x [3] ^ 2 * x [2] ^ 4 -2 * x [3] ^ 2 * R ^ 2 * r ^ 2-4 * x [1] ^ 2 * x [2] ^ 2 * R ^ 2 + x [1] ^ 6 + x [2] ^ 6 + x [3] ^ 6 + 3 * x [3] ^ 2 * x [1] ^ 4 -4 * x [1] ^ 2 * x [2] ^ 2 * r ^ 2-4 * x [1] ^ 2 * x [3] ^ 2 * r ^ 2 + 2 * R ^ 2 * x [1 ] ^ 2 * г ^ 2 -4 * x [2] ^ 2 * x [3] ^ 2 * r ^ 2 + 2 * R ^ 2 * x [2] ^ 2 * r ^ 2-2 * x [1] ^ 4 * R ^ 2 -2 * x [1] ^ 4 * r ^ 2 + R ^ 4 * x [1] ^ 2 + x [1] ^ 2 * r ^ 4-2 * x [2] ^ 4 * R ^ 2-2 * x [2] ^ 4 * r ^ 2 + R ^ 4 * x [2] ^ 2 + x [2] ^ 2 * r ^ 4 + x [3] ^ 2 * R ^ 4 + x [3] ^ 2 * r ^ 4-2 * x [3] ^ 4 * r ^ 2 + 2 * x [3] ^ 4 * R ^ 2) ^ (1/2); c2 = x [2] * (x [1] ^ 2 + x [2] ^ 2 + x [3] ^ 2-R ^ 2-r ^ 2) / (3 * x [1] ^ 4 * x [ 2] ^ 2 + 3 * x [3] ^ 4 * x [2] ^ 2 + 3 * x [3] ^ 4 * x [1] ^ 2 + 6 * x [3] ^ 2 * x [1] ^ 2 * x [2] ^ 2 + 3 * x [1] ^ 2 * x [2] ^ 4 + 3 * x [3] ^ 2 * x [2] ^ 4 -2 * x [3] ^ 2 * R ^ 2 * r ^ 2-4 * x [1] ^ 2 * x [2] ^ 2 * R ^ 2 + x [1] ^ 6 + x [2] ^ 6 + x [3] ^ 6 + 3 * x [3] ^ 2 * x [1] ^ 4-4 * x [1] ^ 2 * x [2] ^ 2 * r ^ 2-4 * x [1] ^ 2 * x [3] ^ 2 * r ^ 2 + 2 * R ^ 2 * x [1] ^ 2 * r ^ 2 -4 * x [2] ^ 2 * x [3] ^ 2 * r ^ 2 + 2 * R ^ 2 * x [2] ^ 2 * r ^ 2-2 * x [1] ^ 4 * R ^ 2 -2 * x [1] ^ 4 * r ^ 2 + R ^ 4 * x [1] ^ 2 + x [1] ^ 2 * r ^ 4-2 * x [2] ^ 4 * R ^ 2-2 * x [2] ^ 4 * r ^ 2 + R ^ 4 * x [2] ^ 2 + x [2] ^ 2 * r ^ 4 + x [3] ^ 2 * R ^ 4 + x [3] ^ 2 * r ^ 4-2 * x [3] ^ 4 * r ^ 2 + 2 * x [3] ^ 4 * R ^ 2) ^ (1/2); c3 = (x [1] ^ 2 + x [2] ^ 2 + x [3] ^ 2 + R ^ 2-r ^ 2) * x [3] / (3 * x [1] ^ 4 * x [ 2] ^ 2 + 3 * x [3] ^ 4 * x [2] ^ 2 + 3 * x [3] ^ 4 * x [1] ^ 2 + 6 * x [3] ^ 2 * x [1] ^ 2 * x [2] ^ 2 + 3 * x [1] ^ 2 * x [2] ^ 4 + 3 * x [3] ^ 2 * x [2] ^ 4 -2 * x [3] ^ 2 * R ^ 2 * r ^ 2 -4 * x [1] ^ 2 * x [2] ^ 2 * R ^ 2 + x [1] ^ 6 + x [2] ^ 6 + x [3] ^ 6 + 3 * x [3] ^ 2 * x [1] ^ 4 -4 * x [1] ^ 2 * x [2] ^ 2 * r ^ 2 -4 * x [1] ^ 2 * x [3] ^ 2 * r ^ 2 + 2 * R ^ 2 * x [1] ^ 2 * r ^ 2 -4 * x [2] ^ 2 * x [3] ^ 2 * r ^ 2 + 2 * R ^ 2 * x [2] ^ 2 * r ^ 2-2 * x [1] ^ 4 * R ^ 2 -2 * x [1] ^ 4 * r ^ 2 + R ^ 4 * x [1] ^ 2 + x [1] ^ 2 * r ^ 4-2 * x [2] ^ 4 * R ^ 2 -2 * x [2] ^ 4 * r ^ 2 + R ^ 4 * x [2] ^ 2 + x [2] ^ 2 * r ^ 4 + x [3] ^ 2 * R ^ 4 + x [3] ^ 2 * r ^ 4-2 * x [3] ^ 4 * r ^ 2 + 2 * x [3] ^ 4 * R ^ 2) ^ (1/2); return (c (c1, c2, c3));}# средняя кривизна трехмерного тора с внешним радиусом R и внутренним радиусом rHTorus = функция (x, r, R) { return (- (3 * x [1] ^ 4 * r ^ 4 + 4 * x [2] ^ 6 * x [3] ^ 2 + 4 * x [1] ^ 6 * x [2] ^ 2-3) * x [2] ^ 4 * x [3] ^ 2 * R ^ 2 -2 * x [1] ^ 6 * R ^ 2 + 4 * x [1] ^ 2 * x [3] ^ 6 + x [3] ^ 6 * R ^ 2 + 4 * x [2] ^ 4 * R ^ 2 * r ^ 2-x [1] ^ 2 * r ^ 6 -x [2] ^ 2 * r ^ 6 + x [2] ^ 4 * R ^ 4 + 4 * x [2] ^ 2 * x [3] ^ 2 * R ^ 4 + 6 * x [2] ^ 2 * x [3] ^ 2 * r ^ 4 -2 * x [1] ^ 2 * R ^ 2 * r ^ 4-x [1] ^ 2 * R ^ 4 * r ^ 2-9 * x [1] ^ 4 * x [2] ^ 2 * r ^ 2 -9 * x [1] ^ 4 * x [3] ^ 2 * r ^ 2 + 4 * x [1] ^ 4 * R ^ 2 * r ^ 2 + 12 * x [1] ^ 2 * x [3 ] ^ 4 * x [2] ^ 2 -3 * x [2] ^ 6 * r ^ 2 + 4 * x [1] ^ 6 * x [3] ^ 2 + 3 * x [3] ^ 4 * r ^ 4-x [3] ^ 4 * R ^ 4 -9 * x [2] ^ 4 * x [3] ^ 2 * r ^ 2 + 2 * x [2] ^ 2 * x [3] ^ 2 * R ^ 2 * r ^ 2 + 4 * x [1 ] ^ 2 * x [2] ^ 6 -6 * x [1] ^ 2 * x [3] ^ 2 * x [2] ^ 2 * R ^ 2-x [3] ^ 2 * r ^ 6 + 6 * x [2] ^ 4 * x [ 3] ^ 4 + x [3] ^ 8 + x [1] ^ 8 + x [2] ^ 8-3 * x [1] ^ 6 * r ^ 2 + 6 * x [1] ^ 4 * x [3] ^ 4 + 12 * x [1] ^ 2 * x [3] ^ 2 * x [2] ^ 4 -6 * x [1] ^ 2 * x [2] ^ 4 * R ^ 2-2 * x [3] ^ 4 * R ^ 2 * r ^ 2-2 * x [2] ^ 2 * R ^ 2 * r ^ 4-x [2] ^ 2 * R ^ 4 * r ^ 2 -9 * x [2] ^ 2 * x [3] ^ 4 * r ^ 2 + x [3] ^ 2 * R ^ 2 * r ^ 4 + x [3] ^ 2 * R ^ 4 * r ^ 2 -9 * x [1] ^ 2 * x [2] ^ 4 * r ^ 2 + 2 * x [1] ^ 2 * R ^ 4 * x [2] ^ 2 + 6 * x [1] ^ 2 * x [2] ^ 2 * r ^ 4-3 * x [1] ^ 4 * х [3] ^ 2 * R ^ 2 -6 * x [1] ^ 4 * x [2] ^ 2 * R ^ 2 + 4 * x [1] ^ 2 * x [3] ^ 2 * R ^ 4 + 6 * x [1] ^ 2 * х [3] ^ 2 * г ^ 4 -9 * x [1] ^ 2 * x [3] ^ 4 * r ^ 2 + 8 * x [1] ^ 2 * R ^ 2 * x [2] ^ 2 * r ^ 2 + 2 * x [1 ] ^ 2 * x [3] ^ 2 * R ^ 2 * r ^ 2 + x [1] ^ 4 * R ^ 4-3 * x [3] ^ 6 * r ^ 2-2 * x [2] ^ 6 * R ^ 2 + 6 * x [1] ^ 4 * x [2 ] ^ 4-x [3] ^ 2 * R ^ 6 -18 * x [1] ^ 2 * x [2] ^ 2 * x [3] ^ 2 * r ^ 2 + 4 * x [2] ^ 2 * x [3] ^ 6 + 12 * x [1] ^ 4 * x [3] ^ 2 * x [2] ^ 2 + 3 * x [2] ^ 4 * r ^ 4) / (3 * x [1] ^ 4 * x [2] ^ 2 + 3 * x [3] ^ 4 * x [2] ^ 2 + 3 * x [3] ^ 4 * x [1] ^ 2 + 6 * x [3] ^ 2 * x [1] ^ 2 * x [2] ^ 2 + 3 * x [1] ^ 2 * x [2] ^ 4 + 3 * x [3] ^ 2 * x [2] ^ 4 -2 * x [3] ^ 2 * R ^ 2 * r ^ 2-4 * x [1] ^ 2 * x [2] ^ 2 * R ^ 2 + x [1] ^ 6 + x [2] ^ 6 + x [3] ^ 6 + 3 * x [3] ^ 2 * x [1] ^ 4-4 * x [1] ^ 2 * x [2] ^ 2 * r ^ 2 -4 * x [1] ^ 2 * x [3] ^ 2 * r ^ 2 + 2 * R ^ 2 * x [1] ^ 2 * r ^ 2 -4 * x [2] ^ 2 * x [3] ^ 2 * r ^ 2 + 2 * R ^ 2 * x [2] ^ 2 * r ^ 2-2 * x [1] ^ 4 * R ^ 2 -2 * x [1] ^ 4 * r ^ 2 + R ^ 4 * x [1] ^ 2 + x [1] ^ 2 * r ^ 4-2 * x [2] ^ 4 * R ^ 2 -2 * x [2] ^ 4 * r ^ 2 + R ^ 4 * x [2] ^ 2 + x [2] ^ 2 * r ^ 4 + x [3] ^ 2 * R ^ 4 + x [3] ^ 2 * r ^ 4-2 * x [3] ^ 4 * r ^ 2 + 2 * x [3] ^ 4 * R ^ 2) ^ (3/2));}# вычислить векторное произведение двух трехмерных векторов# х и у. Нет проверки аргументов по соображениям производительностиcross = function (x, y) { res = rep (0,3); res [1] = x [2] * y [3] - x [3] * y [2]; res [2] = -x [1] * y [3] + x [3] * y [1]; res [3] = x [1] * y [2] - x [2] * y [1]; return (res);}# вычисляем внутреннее произведение двух векторов dim 3# возвращает число, а не матрицу 1x1!точка = функция (х, у) { возврат (сумма (х * у));}# вычислить векторное произведение двух трехмерных векторов# х и у. Нет проверки аргументов по соображениям производительностиcross = function (x, y) { res = rep (0,3); res [1] = x [2] * y [3] - x [3] * y [2]; res [2] = -x [1] * y [3] + x [3] * y [1]; res [3] = x [1] * y [2] - x [2] * y [1]; return (res);}################ main-teilset.seed (280180)et = eulerTorus (c (3,0,0), 3,5,19,10000)torus.write (et, steps = 9000)## bms = euler (c (1,0,0), 4,70000)# euler.write (bms, steps = 10000)

blender3d

Код блендера (python) для создания изображения, которое выглядит почти как это. Поиграйте ...

## импортировать данные из текстового файла matlab и нарисовать BM на S ^ 2## (c) 2007 Кристан Байер и Томас Штайнериз Blender импортировать Curve, Object, Scene, Window, BezTriple, Mesh, Material, Camera,Мириз математического импорта *## import der BM auf der Kugel aus einem csv-filedef importcurve (inpath = "Z.csv"): infile = open (inpath, 'r') lines = infile.readlines () vec = [] для i в строках: li = i.split (',') vec.append ([float (li [0]), float (li [1]), float (li [2] .strip ())]) infile.close () возврат (vec)## function um aus einem vektor (mit den x, y, z Koordinaten) eine Kurve zu machendef vec2Cur (curPts, name = "BMonSphere"): bztr = [] bztr.append (BezTriple.New (curPts [0])) bztr [0] .handleTypes = (BezTriple.HandleTypes.VECT, BezTriple.HandleTypes.VECT) cur = Curve.New (имя) ## TODO wenn es das Objekt schon gibt, dann nicht neu erzeugen cur.appendNurb (bztr [0]) для i в диапазоне (1, len (curPts)): bztr.append (BezTriple.New (curPts [i])) bztr [i] .handleTypes = (BezTriple.HandleTypes.VECT, BezTriple.HandleTypes.VECT) cur [0] .append (bztr [i]) возврат (cur)#erzeugt einen kreis, der später die BM umgibt (в переводе с из-Эбене)def circle (r, name = "tubus"): bzcir = [] bzcir.append (BezTriple.New (0., - r, -4. / 3. * r, 0., - r, 0., 0., - r, 4. / 3. * r)) bzcir [0] .handleTypes = (BezTriple.HandleTypes.FREE, BezTriple.HandleTypes.FREE) cur = Curve.New (имя) ## TODO wenn es das Objekt schon gibt, dann nicht neu erzeugen cur.appendNurb (bzcir [0]) #jetzt alle weietren pkte bzcir.append (BezTriple.New (0., r, 4. / 3. * r, 0., r, 0., 0., r, -4. / 3. * r)) bzcir [1] .handleTypes = (BezTriple.HandleTypes.FREE, BezTriple.HandleTypes.FREE) cur [0] .append (bzcir [1]) bzcir.append (BezTriple.New (0., - r, -4. / 3. * r, 0., - r, 0., 0., - r, 4. / 3. * r)) bzcir [2] .handleTypes = (BezTriple.HandleTypes.FREE, BezTriple.HandleTypes.FREE) cur [0] .append (bzcir [2]) возврат (cur)#erzeuge mit skript eine (glas) kugel (UVSphere)def sphGlass (r = 1.0, name = "Glaskugel", n = 40, smooth = 0): glass = Mesh.New (имя) ## TODO wenn es das Objekt schon gibt, dann nicht neu erzeugen для i в диапазоне (0, n): для j в диапазоне (0, n): x = sin (j * pi * 2.0 / (n-1)) * cos (-pi / 2.0 + i * pi / (n-1)) * 1.0 * r y = cos (j * pi * 2,0 / (n-1)) * (cos (-pi / 2,0 + i * pi / (n-1))) * 1,0 * r z = sin (-pi / 2.0 + i * pi / (n-1)) * 1.0 * r glass.verts.extend (x, y, z) для i в диапазоне (0, n-1):  для j в диапазоне (0, n-1): glass.faces.extend ([i * n + j, i * n + j + 1, (i + 1) * n + j + 1, (i + 1) * n + j]) glass.faces [i * (n-1) + j] .smooth = 1 возврат (стекло)def torus (r = 0,3, R = 1,4):  krGro = круг (r = R, name = "grTorusKreis") #jetzt das material änderndef verglasen (сетка): matGlass = Material.New ("glas") ## TODO wenn es das Objekt schon gibt, dann nichtNeu Erzeugen # matGlass.setSpecShader (0.6) matGlass.setHardness (30) # для спецификации: 30 matGlass.setRayMirr (0,15) matGlass.setFresnelMirr (4.9) matGlass.setFresnelMirrFac (1.8) matGlass.setIOR (1.52) matGlass.setFresnelTrans (3.9) matGlass.setSpecTransp (2.7) # glass.materials.setSpecTransp (1.0) matGlass.rgbCol = [0,66, 0,81, 0,85] matGlass.mode | = Material.Modes.ZTRANSP matGlass.mode | = Material.Modes.RAYTRANSP # matGlass.mode | = Material.Modes.RAYMIRROR mesh.materials = [matGlass] возврат (сетка)def maleBM (сетка): matDraht = Material.New ("roterDraht") ## TODO wenn es das Objekt schon gibt, dannNicht Neu Erzeugen matDraht.rgbCol = [1.0, 0.1, 0.1] mesh.materials = [matDraht] возврат (сетка)#eine solide Mesh-Ebene (Квадер)# auf der höhe ebh, dicke d, seitenlänge (quadratisch) 2 * grdef ebene (ebh = -2,5, d = 0,1, gr = 6,0, name = "Schattenebene"): quader = Mesh.New (имя) ## TODO wenn es das Objekt schon gibt, dann nicht neu erzeugen #obere ebene quader.verts.extend (gr, gr, ebh) quader.verts.extend (-gr, gr, ebh) quader.verts.extend (-gr, -gr, ebh) quader.verts.extend (gr, -gr, ebh) #untere ebene quader.verts.extend (gr, gr, ebh-d) quader.verts.extend (-gr, gr, ebh-d) quader.verts.extend (-gr, -gr, ebh-d) quader.verts.extend (gr, -gr, ebh-d) quader.faces.extend ([0,1,2,3]) quader.faces.extend ([0,4,5,1]) quader.faces.extend ([1,5,6,2]) quader.faces.extend ([2,6,7,3]) quader.faces.extend ([3,7,4,0]) quader.faces.extend ([4,7,6,5]) #die ebene einfärben matEb = Material.New ("ebenen_material") ## TODO wenn es das Objekt schon gibt, dannNicht Neu Erzeugen matEb.rgbCol = [0,53, 0,51, 0,31] matEb.mode | = Material.Modes.TRANSPSHADOW matEb.mode | = Material.Modes.ZTRANSP quader.materials = [matEb] возврат (квадроцикл)####################### основная часть ##### wechsel в режиме редактирования deneditmode = Window.EditMode ()если режим редактирования: Window.EditMode (0)dataBMS = importcurve ("C: / Dokumente und Einstellungen / thire / Desktop / bmsphere / Z.csv")# dataBMS = importcurve ("H: \ MyDocs \ Sphere \ Z.csv")BMScur = vec2Cur (dataBMS, "BMname")# dataStereo = importcurve ("H: \ MyDocs \ сфера \ stZ.csv")# stereoCur = vec2Cur (dataStereo, "SterName")cir = круг (r = 0,01)стекло = sphGlass ()стекло = verglasen (стекло)ebe = ebene ()#jetzt alles hinzufügenscn = Scene.GetCurrent ()obBMScur = scn.objects.new (BMScur, «BMonSphere»)obcir = scn.objects.new (cir, «круглый»)obgla = scn.objects.new (стекло, «Гласкугель»)obebe = scn.objects.new (ebe, "Эбене")# obStereo = scn.objects.new (stereoCur, "StereoCurObj")BMScur.setBevOb (obcir)BMScur.update ()BMScur = мужскойBM (BMScur)# stereoCur.setBevOb (obcir)# stereoCur.update ()cam = Object.Get ("Камера") # cam.setLocation (-5., 5.5, 2.9) # cam.setEuler (62.0, -1., 222.6)#alternativ, besser ??cam.setLocation (-3,3, 8,4, 1,7) cam.setEuler (74,0,200)world = World.GetCurrent ()world.setZen ([0.81,0.82,0.61])world.setHor ([0.77,0.85,0.66])if editmode: Window.EditMode (1) # необязательно, zurück n den letzten modus #ergebnis von# set.seed (24112000)# sbm = euler (c (0,0, -1), T = 1,5, n = 5000)# euler.write (sbm)