21-Oct-2013


1> sqrt(21 + 20i)
ans =  5 + 2i


2> (-1)^(1/3)
ans =  0.50000 + 0.86603i


3> sqrt(3)/2
ans =  0.86603


4> p = [1 0 0 -1]
p =

   1   0   0  -1



5> roots(p)
ans =

  -0.50000 + 0.86603i
  -0.50000 - 0.86603i
   1.00000 + 0.00000i



6> p1 = [1 0 0 -1]
p1 =

   1   0   0  -1



7> p2 = [1 -3 3 -1]
p2 =

   1  -3   3  -1



8> roots(p2)
ans =

   1.00000 + 0.00001i
   1.00000 - 0.00001i
   0.99999 + 0.00000i



9> p = conv(p1,p2)
p =

   1  -3   3  -2   3  -3   1



10> roots(p)
ans =

  -0.50000 + 0.86603i
  -0.50000 - 0.86603i
   1.00010 + 0.00010i
   1.00010 - 0.00010i
   0.99990 + 0.00010i
   0.99990 - 0.00010i



11> [q,r] = deconv(p1,p2)
q =  1
r =

   0   3  -3   0



12> conv(q,p2) + r
ans =

   1   0   0  -1



13> p1
p1 =

   1   0   0  -1



14> polyder(p1)
ans =

   3   0   0



15> polyint(ans)
ans =

   1   0   0   0



16> x = linspace(-3,3)$
parse error:

  syntax error

>>> x = linspace(-3,3)$
                      ^



16> x = linspace(-3,3);


17> y = polyval(p1,x);


18> plot(x,y)


19> grid on


20> triangolo = [1+i 4+2i 3+4i 1+i]
triangolo =

   1 + 1i   4 + 2i   3 + 4i   1 + 1i



21> plot(triangolo)


22> axis('equal')


23> hold on


24> z = cos(pi/3) + i*sin(pi/3)
z =  0.50000 + 0.86603i


25> plot(z*triangolo, 'r')


26> % girare il triangolo per 123 gradi


26> w = cos(123*pi/180) + i*sin(123*pi/180)
w = -0.54464 + 0.83867i


27> plot(w*triangolo, 'g')


28> quit