%exm052_4.m

p0=sym('(1+sqrt(5))/2')			%p0為(1+sqrt(5))/2	準確值		  
p1=sym((1+sqrt(5))/2)			%p1是(1+sqrt(5))/2在數值環境下的近似值
pd=double(p0)						%雙精度數值
e01=vpa(abs(p0-p1))				%在符號環境32位省缺精度下，觀察p0, p1間誤差  
e0d=vpa(abs(p0-pd))				%在符號環境32位省缺精度下，觀察雙精度表達的誤差  

p2=vpa(p0)							%p2是p0在32位相對精度下的近似
e02=vpa(abs(p0-p2),64)			%在64位相對精度下，觀察p0, p2間的誤差  

digits								%驗證vpa變精度計算不影響全局的符號數值計算精度