Помогите пож-та решить

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