循环 fsolve Scilab

Question

举一个愚蠢的例子，假设我希望求解以下非线性方程 x^2 - F(c)=0，，其中 c 可以取零之间的不同值1 和 F 是标准正态 CDF。如果我希望求解 c 的一个特定值，我将使用以下代码：

c = linspace(0,1,100);
L = length(c);
x0 = c;

function Y = eq(x)
Y = x^2 - cdfnor("PQ",x-c(1),0,1)
endfunction

xres = fsolve(x0(1),eq);

我的问题是：是否有一种方法可以求解 c 的每个值的方程> （不仅是 c(1)）？具体来说，我是否可以在 fsolve 上使用循环？如果是这样，怎么办？

Answer 1

只需像这样修改你的脚本：

c = linspace(0,1,100);
L = length(c);
x0 = c;

function Y = eq(x)
Y = x^2 - cdfnor("PQ",x-c,zeros(c),ones(c))
endfunction

xres = fsolve(x0,eq);