matlab解线性方程组误差,solve和vpasolve求方程组的解，结果回代误差太大

clear all

syms  zeta tau

%%

equation1 =(168514287762519778740904808620597809997925488331217960960*tau + 4215970620146567303196271594494894273792416453160010252288*zeta + ...

49069832989606748928554076555925565227996506547203558784208597521633640448*tau*zeta - 18889465931478580854784*((11639035802365446*tau + 1)^2*(49091754281746008067561525039722523508696171273077134502560478654368580*zeta^2 +...

338115057091598301085669127819556414523581734361003372*zeta + 582184304506291015760106403571864841))^(1/2) - ....

1452722665638049521713711129503654072625069044793555211229682833293312*tau^2 + 14489095601885742910250089179264550305792)/(45774676773659427*(11639035802365446*tau + 1)^2);

%%

equation2 =(551118856547494093354321616239628439463804654476565957656576*tau + 118031194942511918809533500831098156848620806375859265130725376*zeta +...

1373769303731871582902838845772816982282030556758436934576755126498874617757696*tau*zeta - 9671406556917033397649408*2^(1/2)*((11639035802365446*tau + 1)^2 ....

*(74453248123958731088127397075069478101157125160850239996572852582652272756*zeta^2 + 59739372207277027450843719764156381913455164557205022844*zeta +...

11983334111163112989777939894185573621))^(1/2) - 52400849408184445829884242460628371886872521214260600984324000000000000*tau^2 +...

47351289923827670190652630617001582449094656)/(18992129006496957881*(11639035802365446*tau + 1)^2);

[tau,zeta]=vpasolve(equation1==0,equation2==0,[tau,zeta]);

%[tau,zeta]=solve(equation1==0,equation2==0,[tau,zeta]);

zeta=double(zeta)

tau=double(tau)