Plot 如何将此方程转换为倍频程代码和绘图

我在Desmos网站上创建了一个方程式

我使用下面的代码尝试在八度音阶中重新创建它。但当我画出来的时候，结果就不一样了。如何将代码固定在八度（如果可能的话，不改变主方程），使其看起来像Desmos图像

x = linspace(0,1,20);
y = linspace(0,1,20);
S=[13.2];

T=1.12;
for zz=1:1:length(S)
  eq1=exp(S(1,zz)*T*log(x))+exp(S(1,zz)*T*log(y));
  hold on
  plot(x,eq1)
  zz
end

PS：我使用的是倍频程4.2.2Desmos.com不绘制（x，eq1），而是（x，y），约束条件是x，y满足给定的方程。因此，您可以为x的每个值求解y，并绘制对（x，y）。
由于log（x），log（y）存在，x和y都大于0（否则您将不得不为x绘图Desmos.com不绘制（x，eq1），而是（x，y），约束条件是x，y满足给定的方程。因此，您为x的每个值求解y，并绘制对（x，y）。
因为log（x），log（y）存在，所以x和y都大于0（否则您必须为x

clear; clc;
x = linspace(0,1,150);

S = 13.2;
T = 1.12;

y = zeros(size(x));
for i = 1:length(x)
  y(i) = (1-exp(S*T*log(x(i))))^(1/S/T);
end
plot(x,y)  

 S = 13.2;
 T = 1.12;
 f = @(x)exp(log(1-exp(S*T*log(x)))./(S*T));
 fplot(f, [0, 1])