DECLARE FUNCTION b (n) DECLARE FUNCTION a (n) CONST pi = 3.141592654# CLS : SCREEN 12 LOCATE 1, 1: PRINT "Universidad Nacional de Educacion a distancia (UNED)." LOCATE 3, 1: PRINT "Departamento de ciencias y tecnicas fisico - quimicas." LOCATE 5, 1: PRINT "Asignatura Metodos teoricos de la quimica - fisica." LOCATE 7, 1: PRINT "Alumno: Francisco Cabot Pol n§ matricula 09-96-00766" LOCATE 11, 1: PRINT "Generamos la secuencia del trabajo numeros aleatorios uniformemente distribuidos en el intervalo(0, 1) mediante el algoritmo de las congruencias y con un periodo de 1024 elementos" LOCATE 15, 1: PRINT "Si el periodo es de 1024 (2^10) entonces; 2(B-2)=2^10 luego B=12 y por tanto M=B^12 y como vamos a tomar una lamda=5 tendremos que buscar numeros comgruentes de lamda elevado a n respecto el modulo 4096" LOCATE 19, 1: PRINT "Y una vez generada vamos a aplicar esta secuencia a la integracion del sen en medio ciclo teniendo en cuenta que la constante de normalizacion sera de pi/2 nos tendra que dar la integracion 1+- una desviacion que despues compararemos con la desviacion si aplicamos RND" DO WHILE INKEY$ = "" LOCATE 27, 1: PRINT "presione cualquier tecla y continue" LOOP CLS LINE (0, 150)-(500, 150) c = 0 FOR n = 0 TO 1023 STEP 1 IF b(n) > 1 THEN LOCATE 1, 1: PRINT n + 1, b(n) END IF y = (SIN(pi * b(n))) PSET (b(n) * 500, 150 - 100 * y) LINE (b(n) * 500, 150)-(b(n) * 500, 150 - 100 * y) c = c + y NEXT n LOCATE 10, 1: PRINT pi * c / 2048 LOCATE 11, 1: PRINT "final primera secuencia de numeros aleatorios ahora veamos con RND" DO WHILE INKEY$ = "" LOCATE 12, 1: PRINT "presione cualquier tecla y continue" LOOP LINE (0, 250)-(500, 250) c = 0 FOR n = 0 TO 1023 STEP 1 x = RND y = (SIN(pi * x)) PSET (x * 500, 250 - 100 * y) LINE (x * 500, 250)-(x * 500, 250 - 100 * y) c = c + y NEXT n LOCATE 17, 1: PRINT pi * c / 2048, "secuencia RND.............Luego la primera secuencia presenta menor desviacion que la RND" DO WHILE INKEY$ = "" LOCATE 18, 1: PRINT "presione cualquier tecla y continue" LOOP c = 0 FOR x = 0 TO 1 STEP .001 y = (SIN(pi * x)) PSET (x * 500, 350 - 100 * y) LINE (x * 500, 350)-(x * 500, 350 - 100 * y) c = c + y NEXT x LOCATE 24, 1: PRINT pi * c / 2048, "integrando a traves del incremento regular (0.001) de la variable x" FUNCTION a (n) a = LOG(5 ^ n / 4096) - INT(LOG(5 ^ n / 4096)) END FUNCTION FUNCTION b (n) b = EXP(a(n)) - INT(EXP(a(n))) END FUNCTION