matlab问题求解答

clear;clc
a=6-4sqrt(2);
b=sqrt(2)-1;
n=0;
m=1.0e-8;
while abs(a-pi)>m
b=(1-(1-b^4)^0.25)/(1+(1-b^4)^0.25);
n=n+1;
if n>1
a=a(b+1)^4-(2^n)*(1+b+b^2)*b;
end
end
fprintf ('%.8f',a)
结果为NaN

想知道错在哪，怎么改

是 b 的值太小了，算两次就基本上到 0 了，后面的 a 就没办法了

确定下题目有没问题