a?=?[523;121;2-35];
b?=?[24;12;30];
[rowcol]=size（a）?;
x0?=?zeros（row1）;
d=diag（a）;
l=-tril（a）+diag（a）;
u=-triu（a）+diag（a）;
m=inv（d）*（l+u）;