exm050242_1.m

来自「MATLAB 6.5 从基础讲解其功能」· M 代码 · 共 11 行

digits									%显示省缺符号数值计算相对精度  

p0=sym('(1+sqrt(5))/2');	%p0为(1+sqrt(5))/2	准确值		  
p1=sym((1+sqrt(5))/2)			%p1是(1+sqrt(5))/2在数值环境下的近似值
e01=vpa(abs(p0-p1))				%在符号环境32位省缺精度下，观察p0, p1间误差  

p2=vpa(p0)							%p2是p0在32位相对精度下的近似
e02=vpa(abs(p0-p2),64)		%在64位相对精度下，观察p0, p2间的误差  

digits									%验证vpa变精度计算不影响全局的符号数值计算精度