📄 montecarlo.m
字号:
function root = Montecarlo(f,B,x0,eps)
format long;
if(nargin==3)
eps=1.0e-4;
end
Fx = subs(sym(f),findsym(sym(f)),x0);
while abs(Fx)>eps
Fx = subs(sym(f),findsym(sym(f)),x0);
Fx1 = abs(Fx) + 10;
while (abs(Fx1) >= abs(Fx))
r = 2*B*rand()-B;
x1 = x0 + r;
Fx1 = subs(sym(f),findsym(sym(f)),x1);
end
x0 = x1;
end
root = x1;
format short;
⌨️ 快捷键说明
复制代码
Ctrl + C
搜索代码
Ctrl + F
全屏模式
F11
切换主题
Ctrl + Shift + D
显示快捷键
?
增大字号
Ctrl + =
减小字号
Ctrl + -