import numpy as np import numpy.linalg as lin import ktbc K,T,B,C = ktbc.ktbc(4) print lin.inv(K)
import numpy as np import scipy.linalg as lin import matplotlib.pyplot as plt import ktbc K,T,B,C = ktbc.ktbc(3); print T h = 1./4. discrete = lin.solve( (1./h)**2 * T, [1.,1.,1.] ) discrete = np.insert(discrete, 0, discrete[0]) discrete = np.append(discrete, 0.) K,T,B,C = ktbc.ktbc(4); print T discrete_2 = lin.solve( (1./h**2)*T, [1./2.,1.,1.,1.] ) # add little diff for plotting # grafik ust uste binmesin diye azicik fark ekledik discrete_2 = discrete_2 + 0.01 discrete_2 = np.append(discrete_2, 0.) def u(x): return (1./2.)*(1. - x**2) p1 = plt.plot([u(0.0), u(0.25), u(0.5), u(0.75), u(1.)]) p2 = plt.plot(discrete) p3 = plt.plot(discrete_2) plt.legend([p1,p2,p3], ["analytical solution (analitik cozum)",
import numpy as np import scipy.linalg as lin import ktbc K,T,B,C = ktbc.ktbc(5) u,v=lin.eig(K) print u print 2-np.sqrt(3), 2-1, 2-0, 2+1, 2+np.sqrt(3) print 2*np.ones((5,1)).T - 2*np.cos((np.arange(5)+1) * np.pi/6)
import numpy as np import scipy.linalg as lin import matplotlib.pyplot as plt import ktbc K, T, B, C = ktbc.ktbc(3) print T h = 1. / 4. discrete = lin.solve((1. / h)**2 * T, [1., 1., 1.]) discrete = np.insert(discrete, 0, discrete[0]) discrete = np.append(discrete, 0.) K, T, B, C = ktbc.ktbc(4) print T discrete_2 = lin.solve((1. / h**2) * T, [1. / 2., 1., 1., 1.]) # add little diff for plotting # grafik ust uste binmesin diye azicik fark ekledik discrete_2 = discrete_2 + 0.01 discrete_2 = np.append(discrete_2, 0.) def u(x): return (1. / 2.) * (1. - x**2) p1 = plt.plot([u(0.0), u(0.25), u(0.5), u(0.75), u(1.)])