求大神！！！帮助用matlab解决问题

clear
clc
f=@(Sopt,D,H,Ra)Sopt./D.*(2+Sopt./D)./D./(1+Sopt./D).^(2/3)-2.75*(H./D).^(1/3).*Ra.^(-1/4);
g=@(D,H,Ra)fsolve(@(Sopt)f(Sopt,D,H,Ra),1);%该函数早闭可以对给定的D、H和Ra计算出Spot,例如g(1,5,300)
for n=5:5:20
    [Ra,D,Sopt]=meshgrid(linspace(300,10000),linspace(1,10),linspace(0,100));
  搜渗  figure;
    陆漏裂isosurface(Ra,D,Sopt,f(Sopt,D,n*D,Ra),0);
    grid on;
    xlabel('Ra');
    ylabel('D');
    zlabel('S_{opt}');
    title(strcat('H/D=',num2str(n)));
end