Вот код
- with(plots):
- with(LinearAlgebra):
- with(VectorCalculus):
-
- p:=(5/4)/(1-9/25*(cos(t))^2)^(1/2):
- z:=3/2*sin(5*t):
- spacecurve([p,t,z],t=-Pi..Pi,color=yellow,coords=cylindrical);
- A:=Vector([p*cos(t),p*sin(t),z]):
- a:=diff(A,t):
- x:=p*cos(t):
- y:=p*sin(t):
- z:=z:
- E:=Vector([x/(x^2+y^2+z^2)^(5/2),y/(x^2+y^2+z^2)^(5/2),z/(x^2+y^2+z^2)^(5/2)]);
-
- f:=DotProduct(a,E);
- int(f,t=-Pi..Pi);
Не считает последний определенный интеграл, он просто его пишет символьно, без подсчета значения
что делать?
-- 02.06.2014, 00:22 --Код, который заработал
- with(plots):
- with(LinearAlgebra):
- with(VectorCalculus):
-
-
-
- q:=1;
- p:=(5/4)/(1-9/25*(cos(t))^2)^(1/2):;
- z:=3/2*sin(5*t):
- spacecurve([p,t,z],t=-Pi..Pi,color=yellow,coords=cylindrical);
- A:=Vector([p*cos(t),p*sin(t),z]):
- a:=diff(A,t):
- x:=p*cos(t):
- y:=p*sin(t):
- z:=z:
- E:=Vector([q*x/(x^2+y^2+z^2)^(3/2),q*y/(x^2+y^2+z^2)^(3/2),q*z/(x^2+y^2+z^2)^(3/2)]):
-
- f:=DotProduct(a,E):
- V:=Int(f,t=-Pi..Pi):
- evalf(V);
-- 02.06.2014, 00:24 --а почему в прошлом случае не считал?