Calcular producto escalar y ángulo que forman dos vectores en Pascal

Hoy me he currado un ejercicio de EDI (Estructura de datos y de la Información). Se trata de un recordatorio de Pascal. Tenemos que hacer un programa que calcule el producto escalar y el ángulo que forman dos vectores de dimensión N.

Así parece sencillo, y en realidad no es muy complicado, pero creía que iba a ser un programa bastante cortito y no ha sido tan corto. Al final he tenido que hacer unos cuantos subprogramas: pasar de radianes a grados, producto escalar, módulo, arcCos, y subprogramas para introducir los vectores y mostrarlos por pantalla.

Pulsa en leer más para ver el código completo.

Aquí tenéis el código completo del programa:



<br />PROGRAM Ejercicio1TP;<br /><br />{Autor: Diego Navarro Mellén}<br />{2º I.T.I. Sistemas}<br />{Objetivo: Dados dos vectores de igual dimensión, calcula el<br />producto escalar de ambos y el ángulo que forman}<br /><br />CONST<br />DIMENSION=2; {dimensión del vector}<br />TYPE<br />{declaramos el tipo vector}<br />tVector=array [1..DIMENSION] of integer;<br /><br />VAR {Variables del Programa Principal}<br />vector1,vector2: tVector;<br /><br /><br /><br />FUNCTION RadToDeg(angulo:real):real;<br />{Convierte radianes en grados}<br />BEGIN<br />RadToDeg:=(angulo*180)/PI;<br />END; {RadToDeg}<br /><br />PROCEDURE RellenaVector(VAR v:tVector);<br />VAR<br />aux,i:integer;<br />BEGIN<br />FOR i:=1 TO DIMENSION DO BEGIN<br />write('Introduce la cordenada ',i,' del vector: ');<br />readln(v[i]);<br />END;<br />writeln('El vector ha sido leído correctamente');<br />writeln;<br />END; {RellenaVector}<br /><br />PROCEDURE EscribeVector(v:tVector);<br />VAR<br />i:integer;<br />BEGIN<br />write('[');<br />FOR i:=1 TO (DIMENSION-1) DO BEGIN<br />write(v[i],',');<br />END;<br />write(v[DIMENSION],']');<br />END; {EscribeVector}<br /><br />FUNCTION Potencia(n:real;exp:integer):real;<br />VAR<br />i:integer;<br />aux:real;<br />BEGIN<br />aux:=1;<br />FOR i:=1 TO exp DO BEGIN<br />aux:=aux*n;<br />END;<br />Potencia:=aux;<br />END; {Potencia}<br /><br />FUNCTION ProductEscalar(u,v:tVector):integer;<br />VAR<br />i:integer;<br />aux: integer;<br />BEGIN<br />aux:=0;<br />FOR i:=1 TO DIMENSION DO BEGIN<br />aux:=aux+u[i]*v[i];<br />END; {FOR}<br />ProductEscalar:=aux;<br />END; {ProductEscalar}<br /><br />FUNCTION Modulo(u:tVector):real;<br />VAR<br />i:integer;<br />aux:real;<br />BEGIN<br />aux:=0;<br />FOR i:=1 TO DIMENSION DO BEGIN<br />aux:=aux+(Potencia(u[i],2));<br />END;<br />Modulo:=sqrt(aux);<br /><br />END; {Modulo}<br /><br />FUNCTION ArcCos(coseno:real):real;<br />VAR<br />seno,tangente:real;<br />BEGIN<br />{De la relación sen^2(x)+cos^2(x)=1}<br />seno:=sqrt(1-Potencia(coseno,2));<br />{Calculamos la artan(x) a partir de la tan(x);<br />Al final obtendremos el valor del ángulo}<br />IF (coseno=0) THEN BEGIN<br />writeln('Error: División por 0')<br />END<br />ELSE BEGIN<br />tangente:=seno/coseno;<br />ArcCos:=Arctan(tangente);<br />END;<br /><br />END; {ArcCos}<br /><br />FUNCTION SacaAngulo(u,v:tVector):real;<br />VAR<br />coseno:real;<br />pEscalar:integer;<br />BEGIN<br />pEscalar:=ProductEscalar(u,v);<br />IF (pEscalar<>0) THEN BEGIN<br />coseno:=pEscalar/(Modulo(u)*Modulo(v));<br />SacaAngulo:=RadToDeg(ArcCos(coseno))<br />END<br />ELSE {Si u*v=0}<br />SacaAngulo:=90; {Forman 90º}<br />END; {SacaAngulo}<br /><br />BEGIN {Programa Principal}<br /><br />{Recogemos los datos de los vectores}<br />RellenaVector(vector1);<br />RellenaVector(vector2);<br /><br />{Mostramos los resultados por pantalla}<br />write('El producto escalar de ');<br />EscribeVector(vector1);<br />write(' y ');<br />EscribeVector(vector2);<br />writeln(' es ',ProductEscalar(vector1,vector2));<br />write('Forman un ángulo de ',SacaAngulo(vector1,vector2):0:2);<br />writeln('º');<br /><br />readln;<br />END. {Programa Principal}<br />



Que pena que no me ha respetado la identación del código, espero poder arreglarlo. Si encontráis algún bug o te ha parecido interesante, escríbeme un comentario.