%Script que que gira un ángulo t un cuadrado en el plano
%Ultano Kindelán, Álgebra ETSIME, UPM. 09/2022

%Coordendas x,y,z de los cuatro vértices del cuadrado (el (0,0,0) repetido
%para cerrar el cuadrado).
xc=[0,0,0,0,0];yc=[0,1,1,0,0];zc=[0,0,1,1,0];
%Matriz 3X5 con las coordenadas de los cuatro vértices más el primero repetido.
p1=[xc;yc;zc]
%Dibujo del cuadrado inicial
plot3(p1(1,:),p1(2,:),p1(3,:));

%Defino la matriz de rotación alrededor del eje x un ángulo t en sentido positivo. 
syms t;
Rx = [1 0 0; 0 cos(t) -sin(t); 0 sin(t) cos(t)]

%Giro el cuadrado un ángulo t=45º en sentido positivo (antihorario).
Rxp1 = Rx*p1;
p1x45=subs(Rxp1, t, pi/4);

%Dibujo el cuadrado girado
hold on;
view([90,0])
plot3(p1x45(1,:),p1x45(2,:),p1x45(3,:),'r');
title('Giro de 45º mediante uan aplicación lineal');