>> x =[125.1200 126.6500 128.6300 129.5100 132.8900 136.7400 137.0100 144.0100 150.0300 155.0600];
>> y=[26.5000 37.9300 50.9300 53.4000 64.2000 74.6500 77.4000 87.5400 95.5000 101.8600];
>> fun=@(pb,x) 0.94*(x.^pb(1)-pb(2)^pb(1)).^(1/pb(1));
>> pb=nlinfit(x,y,fun,[1 2])
pb =
2.6227 124.0240
p=2.6227 b=124.0240