with(linalg): v := vector([x,y,z]); v[1]+v[3]^2; u := vector(3); evalm(u); uPlusv := u + v; evalm(uPlusv); w := vector(7, i -> i^2-1); evalm( a*v ); innerprod(v,u); angle( [1,1] , [0,1] ); convert(",degrees); evalf( convert( angle( [-2,sqrt(2),1] , [2,-1,1] ) , degrees ) );