интервальная система линейных уравнений

[Mirik!]

программа должна содержать решение интервальных систем линейных уравнений Ax=b методом Гаусса.(без отыскания максимального элемента- всегда на главной диагонали)
a i-тое j-тое=0.01*N+ln(i+j), a i-тое j-тое=а j-тое i-тое
а i-тое i-тое=31+N*sin(i), i=j
b i-тое = (12*N)/ln(6+i)

N=2
5 уравнений
δ=0.01

a i-тое j-тое=> [a i-тое j-тое - δ, a i-тое j-тое+δ]
b i-тое => [b i-тое - δ, b i-тое + δ ]