%-- 15.6.10 19:01 --%
maple()
residue
help residue
%-- 16.6.10 18:41 --%
residue('(3z+4)/(z*(z+2^2))
residue('(3z+4)/(z*(z+2^2))')
help residue
residue('(3z+4)','(z*(z+2^2))')
residue('(3z+4)','(z*(z+2)^2)')
%-- 25.9.10 12:15 --%
syms p q r
A=[0 -r q; r 0 -p; -q p 0]
I=eye(3)
A'
A*I*-A
%-- 26.9.10 20:33 --%
syms x
2*x^4+3*x^2+5*x+1
f=ans
g=x^2-5
f*x
f*g
pformat
syms a,b,c,d
syms a b c d
A=[0 1 0 0; 0 0 1 0; 0 0 0 1; a b c d]
det(A)
I=eye(4)
ds
s
syms s
det(s*I-A)
f*x^2
f*g
ans
A
B=[0 1 0 0; 0 0 1 0; 0 0 0 1; -18 9 11 -1]
eig(B)
B=[0 1 0 0; 0 0 1 0; 0 0 0 1; -24 -22 7 4]
eig(B)
f/g
polyx
polyvalm
help polyvalm
polyvalm(f*g)
polyvalm(f*g,x)
polyfit
help polyfit
tf
tf([2 0 -7 5 -14 -24 -2], [1 0 -5])
deconv([2 0 -7 5 -14 -24 -2], [1 0 -5])
[f,r]=deconv([2 0 -7 5 -14 -24 -2], [1 0 -5])
conv([2 0 3 5 1], [ 1 0 -5])
%-- 3.10.10 16:58 --%
R = [0.8047   -0.5059   -0.3106
0.3106    0.8047   -0.5059
0.5059    0.3106    0.8047];
% less approximate rotation
[U,D,V] = svd(R);
R = U*V';
% translation
o = [1;1;1];
R
plot3([0 1], [0 1] [0 1)
plot3([0 1], [0 1], [0 1)
plot3([0 1], [0 1], [0 1])
plot3([0 0], [0 0], [0 1])
hold
plot3([0 0], [0 1], [0 0])
plot3([0 1], [0 0], [0 0])
R^-1
(R^-1)*[1;1;1]
ob=ans
plot3([ob(1) ob(1)+0.8047], [ob(2) ob(2)-0.5059], [ob(3) ob(3)-0.3106])
plot3([ob(1) ob(1)+0.3106], [ob(2) ob(2)+0.8047], [ob(3) ob(3)-0.5059])
plot3([ob(1) ob(1)+0.5059], [ob(2) ob(2)+0.3106], [ob(3) ob(3)+0.8047])
help plot3
plot3([ob(1) 1], [ob(2) 2], [ob(3) 3],'r')
plot3([0 1], [0 2], [0 3],'g')
R*[1;2;3]-[1;1;1]
x2=ans
ob
(R^-1)*x2-ob
plot3([ob(1) ans(1)], [ob(2) ans(2)], [ob(3) ans(3)],'r')
plot3([ob(1) ans(1)], [ob(2) ans(2)], [ob(3) ans(3)],'g')
%-- 3.10.10 18:42 --%