Exponenta.ru

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

Книга http://www.aladjev.narod.ru/Mathematica.htm более доступна и, скорее всего, содержит более развитые средства для обеспечения алгебраических подстановок. Используя том факт, что процедуры Mathematica относительно несложно конвертируются в Maple, подобные подходы можно использовать и в последнем...

Книга http://www.aladjev.narod.ru/Mathematica.htm более доступна и, скорее всего, содержит более развитые средства для обеспечения алгебраических подстановок.

Ledi88, прошу извинить за краткось по причине срочного отъезда утром на длительный срок. В прикрепленном файле найдете документ, где показана попытка упростить ваш интеграл путем замены переменных. Уверен, такой подход вам известен. В исходном виде интеграл не берется ни Maple, ни Mathematica, но и...

Рад бы помочь, да мне не известна функция f.

Вщ-первых, неверно закодировано Pi, во-вторых, интеграл int(f(b+i*t1)*exp((b+i*t1)*t),t1=-T..T); не берется ни в Maple 8-11, ни в Mathematica 7. Не известна функция f. Естественно, и предел вы не вычислите. Поэтому и результат возвращается невычисленным. Либо проверьте задачу, приводящую к такой кон...

Процедура приведена ниже. Правда, увидел уже решение, но решил привести и свое, использующее лишь базовые средства. Я уже и не рад, что ввязался в эту задачку, ниже станет яснее почему. > Test4 := proc(m::posint) local a, b, c, d; a, b := nextprime(10^(m - 1)), 0; do c := length(a); if m < c then re...

Более быстрая процедура на основе предыдущей. > Test3 := proc(m::posint) local a, b, c; a, b := nextprime(10^(m - 1)), 0; do c := length(a); if m < c then return b elif nops({op(convert("" || a, list))}) = c then b := b + 1 end if; a := nextprime(a) end do end proc: > Test3(6); 10239

Прошу извинить, но пкрвые 2 процедуры решают несколько иную задачу - находят все простые числа требуемого свойства с длиной не более 10. Ниже то, что вам нужно > Test2 := proc(m::posint) local a, b, c; a, b := 10^(m - 1), 0; do a := a + 1; c := length(a); if m < c then return b elif type(a, prime) a...

Несколько более быстрая процедура > Test1 := proc(m::posint) local a, b, c, d; a, d := 0, 0; do a := a + 1; b := ithprime(a); c := length(b); if m < c then return d elif nops({op(convert("" || b, list))}) = c then d := d + 1 end if end do end proc; > Test1(6); 13399

> Test := proc(m::posint) local a, b, c; a, b := 0, 0; do a := a + 1; c := length(a); if m < c then return b elif type(a, prime) and nops({op(convert(cat("", a), list))}) = c then b := b + 1 end if end do end proc; > Test(5); 3160

Вероятно, вы не писали сложных процедур коли таковыми считаете приведенные. Простые и кодируются очень просто, например: > P:= (x,y) -> `if`(x<y, sin(x)+cos(y)+sin(x)*cos(y), sin(x)-cos(y)-sin(x)*cos(y)); > P1:= (x) -> `if`(x<2, sin(x)+cos(x), sin(x)-cos(x)); > plot3d(P(x,y), x=0..Pi, y=0..Pi); > pl...

Неверное утверждение о единственно возможных способах вывода графиков процедур. Один я уже привел. Имеются и другие. Наиболее типичные см. в примерах > P:= proc(x,y) sin(x)+cos(y)+sin(x)*cos(y) end proc; > plot3d(P(x,y), x=0..Pi, y=0..Pi); > P1:= proc(x) sin(x)+cos(x) end proc; > plot(P1(x), x=0..Pi...

Кстати, определения процедур также следует кодировать верно, например, Dipole := (n1,n2,l1,l2,R)-> evalf(int(Func(n1, l1, R, x)*x^3*Func(n2, l2, R, x), x = 0 .. R))^2; Разберитесь также с реализованной в процедурах логикой. Вообще говоря, такиго типа алгоритмы легко и эффективнее реализуются единой ...

> h:=0.01: plots[listplot]([seq(StrDipole(1, 0, 2, 1, p), p = [(1+k*h)$k=0..4/h])]);