from oppgave2 import generateA
import numpy as np
from scipy.linalg import solve

# Gitte bjelkeparametre.
length = 2.0
width = 0.3
thickness = 0.03
density = 480.0  # kg/m^3
# Andre konstanter.
g = 9.8  # gravity, m/s^2
E = 1.3 * pow(10, 10)  # Materialkonstanten N/m^2
I = (width * pow(thickness, 3)) / 12.0 #Arealmomentet.
f = - density * width * thickness * g


# Losning paa Ay=b med storrelse (n)
def solve_3(n):
    # Deler bjelkens lengde saa alle segmentene er like.
    h = length / n
    # Genererer baandmatrise (A matrisen)
    matrixA = generateA(n)
    # Genererer b-matrisen.
    matrixB = np.array([[(pow(h, 4) / (E * I)) * f]] * n)
    # Loser for Y basert paa A og B.
    matrixY = solve(matrixA, matrixB)
    # Finner kondisjonstallet.
    kondisjonstall = np.linalg.cond(matrixA)
    # Returnerer losningen og kondisjonstallet.
    return matrixY, kondisjonstall


if __name__ == '__main__':
    matrixY, kondisjonstall = solve_3(10)
    print(matrixY)

485bc175ffbca72706f0eeb92bef54ce8edfe39f
<mask token>


def solve_3(n):
    h = length / n
    matrixA = generateA(n)
    matrixB = np.array([[pow(h, 4) / (E * I) * f]] * n)
    matrixY = solve(matrixA, matrixB)
    kondisjonstall = np.linalg.cond(matrixA)
    return matrixY, kondisjonstall


<mask token>

<mask token>


def solve_3(n):
    h = length / n
    matrixA = generateA(n)
    matrixB = np.array([[pow(h, 4) / (E * I) * f]] * n)
    matrixY = solve(matrixA, matrixB)
    kondisjonstall = np.linalg.cond(matrixA)
    return matrixY, kondisjonstall


if __name__ == '__main__':
    matrixY, kondisjonstall = solve_3(10)
    print(matrixY)

<mask token>
length = 2.0
width = 0.3
thickness = 0.03
density = 480.0
g = 9.8
E = 1.3 * pow(10, 10)
I = width * pow(thickness, 3) / 12.0
f = -density * width * thickness * g


def solve_3(n):
    h = length / n
    matrixA = generateA(n)
    matrixB = np.array([[pow(h, 4) / (E * I) * f]] * n)
    matrixY = solve(matrixA, matrixB)
    kondisjonstall = np.linalg.cond(matrixA)
    return matrixY, kondisjonstall


if __name__ == '__main__':
    matrixY, kondisjonstall = solve_3(10)
    print(matrixY)

from oppgave2 import generateA
import numpy as np
from scipy.linalg import solve
length = 2.0
width = 0.3
thickness = 0.03
density = 480.0
g = 9.8
E = 1.3 * pow(10, 10)
I = width * pow(thickness, 3) / 12.0
f = -density * width * thickness * g


def solve_3(n):
    h = length / n
    matrixA = generateA(n)
    matrixB = np.array([[pow(h, 4) / (E * I) * f]] * n)
    matrixY = solve(matrixA, matrixB)
    kondisjonstall = np.linalg.cond(matrixA)
    return matrixY, kondisjonstall


if __name__ == '__main__':
    matrixY, kondisjonstall = solve_3(10)
    print(matrixY)

from oppgave2 import generateA
import numpy as np
from scipy.linalg import solve

# Gitte bjelkeparametre.
length = 2.0
width = 0.3
thickness = 0.03
density = 480.0  # kg/m^3
# Andre konstanter.
g = 9.8  # gravity, m/s^2
E = 1.3 * pow(10, 10)  # Materialkonstanten N/m^2
I = (width * pow(thickness, 3)) / 12.0 #Arealmomentet.
f = - density * width * thickness * g


# Losning paa Ay=b med storrelse (n)
def solve_3(n):
    # Deler bjelkens lengde saa alle segmentene er like.
    h = length / n
    # Genererer baandmatrise (A matrisen)
    matrixA = generateA(n)
    # Genererer b-matrisen.
    matrixB = np.array([[(pow(h, 4) / (E * I)) * f]] * n)
    # Loser for Y basert paa A og B.
    matrixY = solve(matrixA, matrixB)
    # Finner kondisjonstallet.
    kondisjonstall = np.linalg.cond(matrixA)
    # Returnerer losningen og kondisjonstallet.
    return matrixY, kondisjonstall


if __name__ == '__main__':
    matrixY, kondisjonstall = solve_3(10)
    print(matrixY)
