Octave - Calcular raiz pelo método Regula-Falsi

danielms

O método de regula falsi (regra falsa) ou falsa posição é um método iterativo de resolução numérica de equações não lineales. O método combina o método de bisseção e o método da secante.

Por: Daniel Moreira dos Santos

Fazendo um ajuste não linear em dados experimentais - FORTRAN 90

iagolira

Olá pessoal, como sou amante do Fortran, resolvi criar um programa que faz um ajuste não linear em dados experimentais. - Eu sei, já existe programas para tais! A grande utilidade é quando usa-se muitos parâmetros a serem determinados, o que é vantajoso em relação aos demais. No programa existe a função PLOT, onde nesta dá-se a entrada da função a fazer o ajuste. Exemplo: FUNCTION PLOT(X,A,Qp) IMPLICIT NONE INTEGER, PARAMETER :: qpl = selected_real_kind(15, 16) INTEGER :: Qp real(KIND=qpl) :: X, A(Qp), PLOT PLOT=1.0_qpl+A(1)*exp(-A(2)*X)*cos(A(3)*X)-& &A(4)*exp(-A(5)*X)*cos(A(6)*X)-& &A(7)*exp(-A(8)*X)*cos(A(9)*X)+& &A(10)*exp(A(11)*X)*cos(A(12)*X)-& &A(13)*exp(-A(14)*X)*cos(A(15)*X) END FUNCTION PLOT Percebe-se que a função PLOT têm 16 parâmetros a serem determinados, então percebe-se que é fácil entrar com os valores. COMPILANDO No terminal digite: gfortran Fit.Date.f90 -o FitDate.x -O3 O "-O3" é opcional, pois é um parâmetro de otimização. EXECUTANDO Ainda no terminal, digite: ./FitDate.x Então, aparecerá uma tela pedidos os arquivo que contém os dados a serem analisados, a quantidade de parâmetros e o erro que você quer cometer. Quanto menor o erro, mais demorado. No programa existe uma variável chamada de "tol" (PARAMETER(tol=0.000000000001)), esta é a precisão do cálculo, então ajuste para suas necessidades.

Por: Iago Lira

Octave - Sistemas lineares por Gauss-Seidel

danielms

Resolver o sistema Ax=b pelo método iterativo de Gauss-Seidel. O método de Gauss-Seidel é um método iterativo para resolução de sistemas de equações lineares. O seu nome é uma homenagem aos matemáticos alemães Carl Friedrich Gauss e Philipp Ludwig von Seidel. É semelhante ao método de Jacobi (e como tal, obedece ao mesmo critério de convergência). É condição suficiente de convergência que a matriz seja estritamente diagonal dominante, i. e., fica garantida a convergência da sucessão de valores gerados para a solução exata do sistema linear.

Por: Daniel Moreira dos Santos

Método de Jacobi em SCILAB

arieldll

Solução computacional para sistemas de equações lineares.

Por: Ariel Galante Dalla Costa

Matriz de Hilbert e resolução de sistemas lineares

danielms

Uma matriz de Hilbert é uma matriz quadrada com a seguinte forma: H_{ij} = \frac{1}{i+j-1} Aqui vamos mostrar como resolver um sistema linear HX=b utilizando o octave, onde H é uma matriz de Hilbert.

Por: Daniel Moreira dos Santos