-
Notifications
You must be signed in to change notification settings - Fork 0
/
main.py
57 lines (43 loc) · 1.29 KB
/
main.py
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
# -*- decoding: utf-8 -*-
from MinimosCuadrados import *
from Lector import *
from GaussJordan import *
from newton import *
import math
lec = Lector("archivo.txt")
m = lec.lee()
print "Hola, este programa obtiene una función polinomial a partir de una función tabular"
newt = newton(m)
a = newt.intervalo
if a != -1:
print "Los incrementos no son constantes, favor de ingresar el orden deseado"
orden = int(raw_input("Favor de introducir el orden"))
if orden > len(m) - 1:
orden = len(m) - 1
minimos = MinimosCuadrados(orden)
print "Las ecuaciones: "
var = minimos.formatearEcuaciones()
matriz = minimos.transfomaEcuaciones(var)
for x in xrange(0,len(var)):
print var[x].printPolinomio()
gauss = GaussJordan(matriz)
m = gauss.GaussJordan(matriz, orden)
n = minimos.obtenerResultados(m)
p = minimos.arregloPolinomio(n)
print "Resultado: "
p.printPolinomio()
else:
newt.DiferenciasFinitas()
orden = newt.k + 1
minimos = MinimosCuadrados(orden)
print "Las ecuaciones: "
var = minimos.formatearEcuaciones()
matriz = minimos.transfomaEcuaciones(var)
for x in xrange(0,len(var)):
print var[x].printPolinomio()
gauss = GaussJordan(matriz)
m = gauss.GaussJordan(matriz, orden)
n = minimos.obtenerResultados(m)
p = minimos.arregloPolinomio(n)
print "Resultado: "
p.printPolinomio()