Az nekem nem jó, ha konkrét értéket adok, mert nekem csak egy teoretikus számítás kell. Egy képlet, amit késõbb fel tudok használni.
Leírom a teljes programot, remélem így többet tudsz majd mondani róla.
#
> restar; with(LinearAlgebra);
> r[1] := `<,>`([r*cos(theta)+u*cos(gamma)*sin(theta), r*sin(theta)+u*cos(gamma)*cos(theta), p*(theta-z*`ϕ`)-u*sin(gamma), 1]);
> n[1] := `<,>`([sin(gamma)*sin(theta), -sin(gamma)*cos(theta), cos(gamma)]);
> e := f*`ϕ`/(2*Pi);
> a := sqrt(R^2-e^2)-R+r;
> M[21] := Matrix(4, 4, {(1, 1) = 1, (1, 2) = 0, (1, 3) = 0, (1, 4) = a, (2, 1) = 0, (2, 2) = -sin(gamma), (2, 3) = -cos(gamma), (2, 4) = 0, (3, 1) = 0, (3, 2) = cos(gamma), (3, 3) = -sin(gamma), (3, 4) = (1/2)*f*`ϕ`/Pi, (4, 1) = 0, (4, 2) = 0, (4, 3) = 0, (4, 4) = 1});
> r[2] := evalm(`&*`(M[21], r[1]));
> M[32] := Matrix(4, 4, {(1, 1) = cos(`ϕ`), (1, 2) = sin(`ϕ`), (1, 3) = 0, (1, 4) = 0, (2, 1) = -sin(`ϕ`), (2, 2) = cos(`ϕ`), (2, 3) = 0, (2, 4) = 0, (3, 1) = 0, (3, 2) = 0, (3, 3) = 1, (3, 4) = 0, (4, 1) = 0, (4, 2) = 0, (4, 3) = 0, (4, 4) = 1});
> r[3] := evalm(`&*`(M[32], r[2]));
> L[21] := Minor(M[21], 4, 4, output = matrix);
> n[2] := Multiply(L[21], n[1]);
> w := `<,>`([0, 0, omega]);
> R[2] := `<,>`([r[2][1], r[2][2], r[2][3]]);
> v[2] := CrossProduct(w, R[2]);
> v[h] := `<,>`([0, p*z*omega*cos(gamma), p*z*omega*sin(gamma)]);
> v[f] := `<,>`([0, 0, f*omega/(2*Pi)]);
> v[r] := `<,>`([-v[f][3]*e/sqrt(R^2-e^2), 0, 0]);
> v[rel] := v[h]+v[f]+v[r]-v[2];
> F := DotProduct(v[rel]^%T, n[2], conjugate = false);
> solve(F, u);
Error, (in simplify/table) too many levels of recursion
>
Egy-két dolgot elég furán ír ki, de remélem érthetõ lesz.
Midnen jót!