matlab subs函数不能用,求助SUBS函数不能替换的问题

f=(6*pi - 6*asin(4*a - (2*b*(h - 3/2))/3) + 6*asin((4*a)/3 + (2*b*(h - 3/2))/3) + (8*(-((9*asin(4*a - (2*b*(h - 3/2))/3))/4 - (9*pi)/8 + (9*(1 - (4*a - 2/3*b*(h - 3/2))^2)^(1/2)*(4*a - (2*b*(h - 3/2))/3))/4)*((9*pi)/8 + (9*asin((4*a)/3 + (2*b*(h - 3/2))/3))/4 + (9*(1 - (4/3*a + 2/3*b*(h - 3/2))^2)^(1/2)*((4*a)/3 + (2*b*(h - 3/2))/3))/4))^(1/2))/3 + (9*pi*asin((666666666666666666666666666666677.0*a - 500000000000000000000000000000000.0*b + 333333333333333333333333333333333.0*b*h + 10540925533894597.773329645148109*a*(- 1.3333333333333333333333333333333*10^33*a^2 - 1.7777777777777777777777777777778*10^33*a*b*h + 2.6666666666666666666666666666667*10^33*a*b - 444444444444444444444444444444444.0*b^2*h^2 + 1.3333333333333333333333333333333*10^33*b^2*h - 1.0*10^33*b^2 + 1.0*10^33)^(1/2) + 222222222222222222222222222222222.0*a^2 + 500000000000000000000000000000000.0)/(222222222222222222222222222222222.0*a^2 + 500000000000000000000000000000000.0) - 1))/16 - 6*(1 - (4*a - 2/3*b*(h - 3/2))^2)^(1/2)*(4*a - (2*b*(h - 3/2))/3) + 6*(1 - (4/3*a + 2/3*b*(h - 3/2))^2)^(1/2)*((4*a)/3 + (2*b*(h - 3/2))/3) + (9*pi*asin((333333333333333333333333333333333.0*b*h - 500000000000000000000000000000000.0*b - 2.0*10^33*a + 10540925533894597.773329645148109*a*(- 1.5555555555555555555555555555556*10^34*a^2 + 5.3333333333333333333333333333333*10^33*a*b*h - 8.0*10^33*a*b - 444444444444444444444444444444444.0*b^2*h^2 + 1.3333333333333333333333333333333*10^33*b^2*h - 1.0*10^33*b^2 + 1.0*10^33)^(1/2) + 222222222222222222222222222222222.0*a^2 + 500000000000000000000000000000000.0)/(222222222222222222222222222222222.0*a^2 + 500000000000000000000000000000000.0) - 1))/16 - (pi*(6*a + (500000000000000000000000000000000.0*b*h - 750000000000000000000000000000000.0*b - 3.0*10^33*a + 15811388300841896.659994467722164*a*(- 1.5555555555555555555555555555556*10^34*a^2 + 5.3333333333333333333333333333333*10^33*a*b*h - 8.0*10^33*a*b - 444444444444444444444444444444444.0*b^2*h^2 + 1.3333333333333333333333333333333*10^33*b^2*h - 1.0*10^33*b^2 + 1.0*10^33)^(1/2) + 333333333333333333333333333333333.0*a^2 + 750000000000000000000000000000000.0)/(222222222222222222222222222222222.0*a^2 + 500000000000000000000000000000000.0) - b*(h - 3/2) - 3/2))/2 + (9*pi^2)/16 + (3*pi*(9/4 - (1/(222222222222222222222222222222222.0*a^2 + 500000000000000000000000000000000.0)*(1.0e33*a - 750000000000000000000000000000000.0*b + 500000000000000000000000000000000.0*b*h + 15811388300841896.659994467722164*a*(- 1.3333333333333333333333333333333e33*a^2 - 1.7777777777777777777777777777778e33*a*b*h + 2.6666666666666666666666666666667e33*a*b - 444444444444444444444444444444444.0*b^2*h^2 + 1.3333333333333333333333333333333e33*b^2*h - 1.0e33*b^2 + 1.0e33)^(1/2) + 333333333333333333333333333333333.0*a^2 + 750000000000000000000000000000000.0) - 3/2)^2)^(1/2)*((1.0*10^33*a - 750000000000000000000000000000000.0*b + 500000000000000000000000000000000.0*b*h + 15811388300841896.659994467722164*a*(- 1.3333333333333333333333333333333*10^33*a^2 - 1.7777777777777777777777777777778*10^33*a*b*h + 2.6666666666666666666666666666667*10^33*a*b - 444444444444444444444444444444444.0*b^2*h^2 + 1.3333333333333333333333333333333*10^33*b^2*h - 1.0*10^33*b^2 + 1.0*10^33)^(1/2) + 333333333333333333333333333333333.0*a^2 + 750000000000000000000000000000000.0)/(3*(222222222222222222222222222222222.0*a^2 + 500000000000000000000000000000000.0)) - 1/2))/4 + (pi*(1 - (4*((1.0*10^33*a - 750000000000000000000000000000000.0*b + 500000000000000000000000000000000.0*b*h + 15811388300841896.659994467722164*a*(- 1.3333333333333333333333333333333*10^33*a^2 - 1.7777777777777777777777777777778*10^33*a*b*h + 2.6666666666666666666666666666667*10^33*a*b - 444444444444444444444444444444444.0*b^2*h^2 + 1.3333333333333333333333333333333*10^33*b^2*h - 1.0*10^33*b^2 + 1.0*10^33)^(1/2) + 333333333333333333333333333333333.0*a^2 + 750000000000000000000000000000000.0)/(222222222222222222222222222222222.0*a^2 + 500000000000000000000000000000000.0) - 3/2)^2)/9)^(1/2)*(2*a + (1.0*10^33*a - 750000000000000000000000000000000.0*b + 500000000000000000000000000000000.0*b*h + 15811388300841896.659994467722164*a*(- 1.3333333333333333333333333333333*10^33*a^2 - 1.7777777777777777777777777777778*10^33*a*b*h + 2.6666666666666666666666666666667*10^33*a*b - 444444444444444444444444444444444.0*b^2*h^2 + 1.3333333333333333333333333333333*10^33*b^2*h - 1.0*10^33*b^2 + 1.0*10^33)^(1/2) + 333333333333333333333333333333333.0*a^2 + 750000000000000000000000000000000.0)/(222222222222222222222222222222222.0*a^2 + 500000000000000000000000000000000.0) - b*(h - 3/2) + 3/2))/2 + (3*pi*(9/4 - (1/(222222222222222222222222222222222.0*a^2 + 500000000000000000000000000000000.0)*(500000000000000000000000000000000.0*b*h - 750000000000000000000000000000000.0*b - 3.0e33*a + 15811388300841896.659994467722164*a*(- 1.5555555555555555555555555555556e34*a^2 + 5.3333333333333333333333333333333e33*a*b*h - 8.0e33*a*b - 444444444444444444444444444444444.0*b^2*h^2 + 1.3333333333333333333333333333333e33*b^2*h - 1.0e33*b^2 + 1.0e33)^(1/2) + 333333333333333333333333333333333.0*a^2 + 750000000000000000000000000000000.0) - 3/2)^2)^(1/2)*((500000000000000000000000000000000.0*b*h - 750000000000000000000000000000000.0*b - 3.0*10^33*a + 15811388300841896.659994467722164*a*(- 1.5555555555555555555555555555556*10^34*a^2 + 5.3333333333333333333333333333333*10^33*a*b*h - 8.0*10^33*a*b - 444444444444444444444444444444444.0*b^2*h^2 + 1.3333333333333333333333333333333*10^33*b^2*h - 1.0*10^33*b^2 + 1.0*10^33)^(1/2) + 333333333333333333333333333333333.0*a^2 + 750000000000000000000000000000000.0)/(3*(222222222222222222222222222222222.0*a^2 + 500000000000000000000000000000000.0)) - 1/2))/4);

subs(f,{a,b,h},{0.05,0.96,2.508});