Matematické Fórum

Zdravím, lámu si hlavu s řešením příkladu do školy..Mám napsat program v matlabu pro řešení soustavy lineárních algebraických rovnic s maticí A s rozměry n x n. (má být 500x500 a všude mají být 2ky a na diagonále mají být 4ky) Vstupním parametrem programu má být vektor pravé strany dané soustavy, výstupem bude řešení soustavy. K řešení soustavy mám použí LU rozklad. S matlabem kamarád zrovna nejsem, můžete mi někdo poradit,jak dostanu do diagonály 4ky a případně jestli mám zbytek vůbec dobře??? Díky moc...

function [ x ] = LUrozklad( b )
A=2*(ones(500,500))
b
x=A\b
[L,U]=lu(A)
y=L\b
z=U\y
end