明年六月一定毕业
2008-11-19, 22:10
:heart: 非线性方程组50个未知数49个方程，据说最好是broyden秩一法
(1/4)+(x51)^2+(x52)^2+(x53)^2+(x54)^2+(x55)^2+(x41)^2+(x42)^2+(x43)^2+(x44)^2+(x45)^2+(x31)^2+(x32)^2+(x33)^2+(x34)^2+(x35)^2+(x21)^2+(x22)^2+(x23)^2+(x24)^2+(x25)^2+(x11)^2+(x12)^2+(x13)^2+(x14)^2+(x15)^2=1/2;
(1/4)+(y51)^2+(y52)^2+(y53)^2+(y54)^2+(y55)^2+(y41)^2+(y42)^2+(y43)^2+(y44)^2+(y45)^2+(y31)^2+(y32)^2+(y33)^2+(y34)^2+(y35)^2+(y21)^2+(y22)^2+(y23)^2+(y24)^2+(y25)^2+(y11)^2+(y12)^2+(y13)^2+(y14)^2+(y15)^2=1/2;
(x51)*(y51)+(x52)*(y52)+(x53)*(y53)+(x54)*(y54)+(x55)*(y55)+(x41)*(y41)+(x42)*(y42)+(x43)*(y43)+(x44)*(y44)+(x45)*(y45)+(x31)*(y31)+(x32)*(y32)+(x33)*(y33)+(x34)*(y34)+(x35)*(y35)+(x21)*(y21)+(x22)*(y22)+(x23)*(y23)+(x24)*(y24)+(x25)*(y25)+(x11)*(y11)+(x12)*(y12)+(x13)*(y13)+(x14)*(y14)+(x15)*(y15)=0;
(x51)*(x15)=0;
(x51)*(x13)+(x52)*(x14)+(x53)*(x15)=0;
(x51)*(x11)+(x52)*(x12)+(x53)*(x13)+(x54)*(x14)+(x55)*(x15)=0;
(x53)*(x11)+(x54)*(x12)+(x55)*(x13)=0;
(x55)*(x11)=0;
(x51)*(x24)+(x52)*(x25)+(x41)*(x14)+(x42)*(x15)=0;
(x51)*(x22)+(x52)*(x23)+(x53)*(x24)+(x54)*(x25)+(x41)*(x12)+(x42)*(x13)+(x43)*(x14)+(x44)*(x15)=0;
(x52)*(x21)+(x53)*(x22)+(x54)*(x23)+(x55)*(x24)+(x42)*(x11)+(x43)*(x12)+(x44)*(x13)+(x45)*(x14)=0;
(x54)*(x21)+(x55)*(x22)+(x44)*(x11)+(x45)*(x12)=0;
(x51)*(x35)+(x41)*(x25)+(x31)*(x15)=0;
(x51)*(x33)+(x52)*(x34)+(x53)*(x35)+(x41)*(x23)+(x42)*(x24)+(x43)*(x25)+(x31)*(x13)+(x32)*(x14)+(x33)*(x15)=0;
(x51)*(x31)+(x52)*(x32)+(x53)*(x33)+(x54)*(x34)+(x55)*(x35)+(x41)*(x21)+(x42)*(x22)+(x43)*(x23)+(x44)*(x24)+(x45)*(x25)+(x31)*(x11)+(x32)*(x12)+(x33)*(x13)+(x34)*(x14)+(x35)*(x15)=0;
(x53)*(x31)+(x54)*(x32)+(x55)*(x33)+(x43)*(x21)+(x44)*(x22)+(x45)*(x23)+(x33)*(x11)+(x34)*(x12)+(x35)*(x13)=0;
(x55)*(x31)+(x45)*(x21)+(x35)*(x11)=0;
(x51)*(x44)+(x52)*(x45)+(x41)*(x34)+(x42)*(x35)+(x31)*(x24)+(x32)*(x25)+(x21)*(x14)+(x22)*(x15)=0;
(x51)*(x42)+(x52)*(x43)+(x53)*(x44)+(x54)*(x45)+(x41)*(x32)+(x42)*(x33)+(x43)*(x34)+(x44)*(x35)+(x31)*(x22)+(x32)*(x23)+(x33)*(x24)+(x34)*(x25)+(x21)*(x12)+(x22)*(x13)+(x23)*(x14)+(x24)*(x15)=0;
(x52)*(x41)+(x53)*(x42)+(x54)*(x43)+(x55)*(x44)+(x42)*(x31)+(x43)*(x32)+(x44)*(x33)+(x45)*(x34)+(x32)*(x21)+(x33)*(x22)+(x34)*(x23)+(x35)*(x24)+(x22)*(x11)+(x23)*(x12)+(x24)*(x13)+(x25)*(x14)=0;
(x54)*(x41)+(x55)*(x42)+(x44)*(x31)+(x45)*(x32)+(x34)*(x21)+(x35)*(x22)+(x24)*(x11)+(x25)*(x12)=0;
(x51)*(x55)+(x41)*(x45)+(x31)*(x35)+(x21)*(x25)+(x11)*(x15)=0;
(x51)*(x53)+(x52)*(x54)+(x53)*(x55)+(x41)*(x43)+(x42)*(x44)+(x43)*(x45)+(x31)*(x33)+(x32)*(x34)+(x33)*(x35)+(x21)*(x23)+(x22)*(x24)+(x23)*(x25)+(x11)*(x13)+(x12)*(x14)+(x13)*(x15)=0;
(y51)*(y15)=0;
(y51)*(y13)+(y52)*(y14)+(y53)*(y15)=0;
(y51)*(y11)+(y52)*(y12)+(y53)*(y13)+(y54)*(y14)+(y55)*(y15)=0;
(y53)*(y11)+(y54)*(y12)+(y55)*(y13)=0;
(y55)*(y11)=0;
(y51)*(y24)+(y52)*(y25)+(y41)*(y14)+(y42)*(y15)=0;
(y51)*(y22)+(y52)*(y23)+(y53)*(y24)+(y54)*(y25)+(y41)*(y12)+(y42)*(y13)+(y43)*(y14)+(y44)*(y15)=0;
(y52)*(y21)+(y53)*(y22)+(y54)*(y23)+(y55)*(y24)+(y42)*(y11)+(y43)*(y12)+(y44)*(y13)+(y45)*(y14)=0;
(y54)*(y21)+(y55)*(y22)+(y44)*(y11)+(y45)*(y12)=0;
(y51)*(y35)+(y41)*(y25)+(y31)*(y15)=0;
(y51)*(y33)+(y52)*(y34)+(y53)*(y35)+(y41)*(y23)+(y42)*(y24)+(y43)*(y25)+(y31)*(y13)+(y32)*(y14)+(y33)*(y15)=0;
(y51)*(y31)+(y52)*(y32)+(y53)*(y33)+(y54)*(y34)+(y55)*(y35)+(y41)*(y21)+(y42)*(y22)+(y43)*(y23)+(y44)*(y24)+(y45)*(y25)+(y31)*(y11)+(y32)*(y12)+(y33)*(y13)+(y34)*(y14)+(y35)*(y15)=0;
(y53)*(y31)+(y54)*(y32)+(y55)*(y33)+(y43)*(y21)+(y44)*(y22)+(y45)*(y23)+(y33)*(y11)+(y34)*(y12)+(y35)*(y13)=0;
(y55)*(y31)+(y45)*(y21)+(y35)*(y11)=0;
(y51)*(y44)+(y52)*(y45)+(y41)*(y34)+(y42)*(y35)+(y31)*(y24)+(y32)*(y25)+(y21)*(y14)+(y22)*(y15)=0;
(y51)*(y42)+(y52)*(y43)+(y53)*(y44)+(y54)*(y45)+(y41)*(y32)+(y42)*(y33)+(y43)*(y34)+(y44)*(y35)+(y31)*(y22)+(y32)*(y23)+(y33)*(y24)+(y34)*(y25)+(y21)*(y12)+(y22)*(y13)+(y23)*(y14)+(y24)*(y15)=0;
(y52)*(y41)+(y53)*(y42)+(y54)*(y43)+(y55)*(y44)+(y42)*(y31)+(y43)*(y32)+(y44)*(y33)+(y45)*(y34)+(y32)*(y21)+(y33)*(y22)+(y34)*(y23)+(y35)*(y24)+(y22)*(y11)+(y23)*(y12)+(y24)*(y13)+(y25)*(y14)=0;
(y54)*(y41)+(y55)*(y42)+(y44)*(y31)+(y45)*(y32)+(y34)*(y21)+(y35)*(y22)+(y24)*(y11)+(y25)*(y12)=0;
(y51)*(y55)+(y41)*(y45)+(y31)*(y35)+(y21)*(y25)+(y11)*(y15)=0;
(y51)*(y53)+(y52)*(y54)+(y53)*(y55)+(y41)*(y43)+(y42)*(y44)+(y43)*(y45)+(y31)*(y33)+(y32)*(y34)+(y33)*(y35)+(y21)*(y23)+(y22)*(y24)+(y23)*(y25)+(y11)*(y13)+(y12)*(y14)+(y13)*(y15)=0;
(x51)+(x53)+(x55)+(x42)+(x44)+(x31)+(x33)+(x35)+(x22)+(x24)+(x11)+(x13)+(x15)+(y51)+(y53)+(y55)+(y42)+(y44)+(y31)+(y33)+(y35)+(y22)+(y24)+(y11)+(y13)+(y15)=1/2;
(x52)+(x54)+(x41)+(x43)+(x45)+(x32)+(x34)+(x21)+(x23)+(x25)+(x12)+(x14)+(y52)+(y54)+(y41)+(y43)+(y45)+(y32)+(y34)+(y21)+(y23)+(y25)+(y12)+(y14)=1/2;
2*(x51)-2*(x55)+(x42)-(x44)+2*(x31)-2*(x35)+(x22)-(x24)+2*(x11)-2*(x15)+3*(y51)+(y53)-(y55)+2*(y42)+3*(y31)+(y33)-(y35)+2*(y22)+3*(y11)+(y13)-(y15)=1/4;
(x52)-(x54)+2*(x41)-2*(x45)-2*(x21)-2*(x25)+(x12)-(x14)+2*(y52)+3*(y41)-(y43)-(y45)+3*(y21)+(y23)-(y25)+2*(y12)=1/4;
2*(x51)+2*(x53)+2*(x55)+(x42)+(x44)-(x22)-(x24)-2*(x11)-2*(x13)-2*(x15)+2*(y51)+2*(y53)+2*(y55)+(y42)+(y44)-(y22)-(y24)-2*(y11)-2*(y13)-2*(y15)=1/4;
2*(x52)+2*(x54)+(x41)+(x43)+(x45)-(x21)-(x23)-(x25)-2*(x12)-2*(x14)+2*(y52)+2*(y54)+(y41)+(y43)+(y45)-(y21)-(y23)-(y25)-2*(y12)-2*(y14)=1/4;
(1/4)+(x51)^2+(x52)^2+(x53)^2+(x54)^2+(x55)^2+(x41)^2+(x42)^2+(x43)^2+(x44)^2+(x45)^2+(x31)^2+(x32)^2+(x33)^2+(x34)^2+(x35)^2+(x21)^2+(x22)^2+(x23)^2+(x24)^2+(x25)^2+(x11)^2+(x12)^2+(x13)^2+(x14)^2+(x15)^2=1/2;
(1/4)+(y51)^2+(y52)^2+(y53)^2+(y54)^2+(y55)^2+(y41)^2+(y42)^2+(y43)^2+(y44)^2+(y45)^2+(y31)^2+(y32)^2+(y33)^2+(y34)^2+(y35)^2+(y21)^2+(y22)^2+(y23)^2+(y24)^2+(y25)^2+(y11)^2+(y12)^2+(y13)^2+(y14)^2+(y15)^2=1/2;
(x51)*(y51)+(x52)*(y52)+(x53)*(y53)+(x54)*(y54)+(x55)*(y55)+(x41)*(y41)+(x42)*(y42)+(x43)*(y43)+(x44)*(y44)+(x45)*(y45)+(x31)*(y31)+(x32)*(y32)+(x33)*(y33)+(x34)*(y34)+(x35)*(y35)+(x21)*(y21)+(x22)*(y22)+(x23)*(y23)+(x24)*(y24)+(x25)*(y25)+(x11)*(y11)+(x12)*(y12)+(x13)*(y13)+(x14)*(y14)+(x15)*(y15)=0;
(x51)*(x15)=0;
(x51)*(x13)+(x52)*(x14)+(x53)*(x15)=0;
(x51)*(x11)+(x52)*(x12)+(x53)*(x13)+(x54)*(x14)+(x55)*(x15)=0;
(x53)*(x11)+(x54)*(x12)+(x55)*(x13)=0;
(x55)*(x11)=0;
(x51)*(x24)+(x52)*(x25)+(x41)*(x14)+(x42)*(x15)=0;
(x51)*(x22)+(x52)*(x23)+(x53)*(x24)+(x54)*(x25)+(x41)*(x12)+(x42)*(x13)+(x43)*(x14)+(x44)*(x15)=0;
(x52)*(x21)+(x53)*(x22)+(x54)*(x23)+(x55)*(x24)+(x42)*(x11)+(x43)*(x12)+(x44)*(x13)+(x45)*(x14)=0;
(x54)*(x21)+(x55)*(x22)+(x44)*(x11)+(x45)*(x12)=0;
(x51)*(x35)+(x41)*(x25)+(x31)*(x15)=0;
(x51)*(x33)+(x52)*(x34)+(x53)*(x35)+(x41)*(x23)+(x42)*(x24)+(x43)*(x25)+(x31)*(x13)+(x32)*(x14)+(x33)*(x15)=0;
(x51)*(x31)+(x52)*(x32)+(x53)*(x33)+(x54)*(x34)+(x55)*(x35)+(x41)*(x21)+(x42)*(x22)+(x43)*(x23)+(x44)*(x24)+(x45)*(x25)+(x31)*(x11)+(x32)*(x12)+(x33)*(x13)+(x34)*(x14)+(x35)*(x15)=0;
(x53)*(x31)+(x54)*(x32)+(x55)*(x33)+(x43)*(x21)+(x44)*(x22)+(x45)*(x23)+(x33)*(x11)+(x34)*(x12)+(x35)*(x13)=0;
(x55)*(x31)+(x45)*(x21)+(x35)*(x11)=0;
(x51)*(x44)+(x52)*(x45)+(x41)*(x34)+(x42)*(x35)+(x31)*(x24)+(x32)*(x25)+(x21)*(x14)+(x22)*(x15)=0;
(x51)*(x42)+(x52)*(x43)+(x53)*(x44)+(x54)*(x45)+(x41)*(x32)+(x42)*(x33)+(x43)*(x34)+(x44)*(x35)+(x31)*(x22)+(x32)*(x23)+(x33)*(x24)+(x34)*(x25)+(x21)*(x12)+(x22)*(x13)+(x23)*(x14)+(x24)*(x15)=0;
(x52)*(x41)+(x53)*(x42)+(x54)*(x43)+(x55)*(x44)+(x42)*(x31)+(x43)*(x32)+(x44)*(x33)+(x45)*(x34)+(x32)*(x21)+(x33)*(x22)+(x34)*(x23)+(x35)*(x24)+(x22)*(x11)+(x23)*(x12)+(x24)*(x13)+(x25)*(x14)=0;
(x54)*(x41)+(x55)*(x42)+(x44)*(x31)+(x45)*(x32)+(x34)*(x21)+(x35)*(x22)+(x24)*(x11)+(x25)*(x12)=0;
(x51)*(x55)+(x41)*(x45)+(x31)*(x35)+(x21)*(x25)+(x11)*(x15)=0;
(x51)*(x53)+(x52)*(x54)+(x53)*(x55)+(x41)*(x43)+(x42)*(x44)+(x43)*(x45)+(x31)*(x33)+(x32)*(x34)+(x33)*(x35)+(x21)*(x23)+(x22)*(x24)+(x23)*(x25)+(x11)*(x13)+(x12)*(x14)+(x13)*(x15)=0;
(y51)*(y15)=0;
(y51)*(y13)+(y52)*(y14)+(y53)*(y15)=0;
(y51)*(y11)+(y52)*(y12)+(y53)*(y13)+(y54)*(y14)+(y55)*(y15)=0;
(y53)*(y11)+(y54)*(y12)+(y55)*(y13)=0;
(y55)*(y11)=0;
(y51)*(y24)+(y52)*(y25)+(y41)*(y14)+(y42)*(y15)=0;
(y51)*(y22)+(y52)*(y23)+(y53)*(y24)+(y54)*(y25)+(y41)*(y12)+(y42)*(y13)+(y43)*(y14)+(y44)*(y15)=0;
(y52)*(y21)+(y53)*(y22)+(y54)*(y23)+(y55)*(y24)+(y42)*(y11)+(y43)*(y12)+(y44)*(y13)+(y45)*(y14)=0;
(y54)*(y21)+(y55)*(y22)+(y44)*(y11)+(y45)*(y12)=0;
(y51)*(y35)+(y41)*(y25)+(y31)*(y15)=0;
(y51)*(y33)+(y52)*(y34)+(y53)*(y35)+(y41)*(y23)+(y42)*(y24)+(y43)*(y25)+(y31)*(y13)+(y32)*(y14)+(y33)*(y15)=0;
(y51)*(y31)+(y52)*(y32)+(y53)*(y33)+(y54)*(y34)+(y55)*(y35)+(y41)*(y21)+(y42)*(y22)+(y43)*(y23)+(y44)*(y24)+(y45)*(y25)+(y31)*(y11)+(y32)*(y12)+(y33)*(y13)+(y34)*(y14)+(y35)*(y15)=0;
(y53)*(y31)+(y54)*(y32)+(y55)*(y33)+(y43)*(y21)+(y44)*(y22)+(y45)*(y23)+(y33)*(y11)+(y34)*(y12)+(y35)*(y13)=0;
(y55)*(y31)+(y45)*(y21)+(y35)*(y11)=0;
(y51)*(y44)+(y52)*(y45)+(y41)*(y34)+(y42)*(y35)+(y31)*(y24)+(y32)*(y25)+(y21)*(y14)+(y22)*(y15)=0;
(y51)*(y42)+(y52)*(y43)+(y53)*(y44)+(y54)*(y45)+(y41)*(y32)+(y42)*(y33)+(y43)*(y34)+(y44)*(y35)+(y31)*(y22)+(y32)*(y23)+(y33)*(y24)+(y34)*(y25)+(y21)*(y12)+(y22)*(y13)+(y23)*(y14)+(y24)*(y15)=0;
(y52)*(y41)+(y53)*(y42)+(y54)*(y43)+(y55)*(y44)+(y42)*(y31)+(y43)*(y32)+(y44)*(y33)+(y45)*(y34)+(y32)*(y21)+(y33)*(y22)+(y34)*(y23)+(y35)*(y24)+(y22)*(y11)+(y23)*(y12)+(y24)*(y13)+(y25)*(y14)=0;
(y54)*(y41)+(y55)*(y42)+(y44)*(y31)+(y45)*(y32)+(y34)*(y21)+(y35)*(y22)+(y24)*(y11)+(y25)*(y12)=0;
(y51)*(y55)+(y41)*(y45)+(y31)*(y35)+(y21)*(y25)+(y11)*(y15)=0;
(y51)*(y53)+(y52)*(y54)+(y53)*(y55)+(y41)*(y43)+(y42)*(y44)+(y43)*(y45)+(y31)*(y33)+(y32)*(y34)+(y33)*(y35)+(y21)*(y23)+(y22)*(y24)+(y23)*(y25)+(y11)*(y13)+(y12)*(y14)+(y13)*(y15)=0;
(x51)+(x53)+(x55)+(x42)+(x44)+(x31)+(x33)+(x35)+(x22)+(x24)+(x11)+(x13)+(x15)+(y51)+(y53)+(y55)+(y42)+(y44)+(y31)+(y33)+(y35)+(y22)+(y24)+(y11)+(y13)+(y15)=1/2;
(x52)+(x54)+(x41)+(x43)+(x45)+(x32)+(x34)+(x21)+(x23)+(x25)+(x12)+(x14)+(y52)+(y54)+(y41)+(y43)+(y45)+(y32)+(y34)+(y21)+(y23)+(y25)+(y12)+(y14)=1/2;
2*(x51)-2*(x55)+(x42)-(x44)+2*(x31)-2*(x35)+(x22)-(x24)+2*(x11)-2*(x15)+3*(y51)+(y53)-(y55)+2*(y42)+3*(y31)+(y33)-(y35)+2*(y22)+3*(y11)+(y13)-(y15)=1/4;
(x52)-(x54)+2*(x41)-2*(x45)-2*(x21)-2*(x25)+(x12)-(x14)+2*(y52)+3*(y41)-(y43)-(y45)+3*(y21)+(y23)-(y25)+2*(y12)=1/4;
2*(x51)+2*(x53)+2*(x55)+(x42)+(x44)-(x22)-(x24)-2*(x11)-2*(x13)-2*(x15)+2*(y51)+2*(y53)+2*(y55)+(y42)+(y44)-(y22)-(y24)-2*(y11)-2*(y13)-2*(y15)=1/4;
2*(x52)+2*(x54)+(x41)+(x43)+(x45)-(x21)-(x23)-(x25)-2*(x12)-2*(x14)+2*(y52)+2*(y54)+(y41)+(y43)+(y45)-(y21)-(y23)-(y25)-2*(y12)-2*(y14)=1/4;