optimization - Monte Carlo simulation with condition loop in MATLAB -


i physicist.

i have written code on monte carlo simulation condition loop.

i getting error in running simulation code. want consider positive value of simulation result. when run code, error.

i running 10000 iterations , have 5 parameters such a,b,c,d , e. generating random number each parameter using variance , mean of every parameter of normal distribution.

the code follows: 

    n = 10000;     constant = 5;     arand = (3*10^(12)*randn(1,n)) + 7*10^(6)*ones(1,n);     brand = (9*randn(1,n)) + 17*ones(1,n);     crand = (2*10^(-4)*randn(1,n)) + 0.2*ones(1,n);     drand = (0.0017*randn(1,n)) + 0.50*ones(1,n);     erand = (0.00004*randn(1,n)) + 1.5*ones(1,n);     if p1 > 0        p1 =  constant*arand.*brand.*crand.*drand.*(1/erand)     end      plot(p1); 

p1 =  constant*arand.*brand.*crand.*drand./erand;

it's not clear p1 before if statement.

note, "if p1 > 0" means "if all(p1 > 0)".


Comments

Post a Comment

Popular posts from this blog

linux - xterm copying to CLIPBOARD using copy-selection causes automatic updating of CLIPBOARD upon mouse selection -

i've tried various ways coerce xterm (versions 285 , 292) copy selection clipboard clipboard whenever press ctrl-shift-c. promising way far has been, in ~/.xresources, put this: xterm*vt100.translations: #override \ ctrl shift <keypress> c: copy-selection(clipboard) \n\ ctrl shift <keypress> v: insert-selection(clipboard) ctrl-shift-v works perfectly, copying has nuance... if restart xterm, highlighting text puts things in primary clipboard; expected, proper, default behavior. if hit ctrl-shift-c, copies current selection clipboard clipboard. the bug, however, if highlight text after first time pressed ctrl-shift-c, you'll see highlighting copies both primary and clipboard. cannot figure out how stop xterm updating clipboard upon selection... , doesn't make sense me in first place. i said, @ specific point in time, copy selection clipboard... yet starts updating automatically upon selection after doing once... anyone have workaround...
Read more

c++ - qgraphicsview horizontal scrolling always has a vertical delta -

Image
i've subclassed qgraphicsview custom canvas used in cad application. i've reimplemented qgraphicsview::wheelevent check keyboard modifiers control key and, if control key pressed, zoom. i'm trying implement horizontal scroll when user holds shift , uses wheel. the problem i'm having horizontal scrolling always scrolls 0.279. not huge problem, hugely annoying , points else being wrong. so, here questions: is right way implement horizontal scrolling? if not, is? how eliminate delta of 0.279? thanks in advance. code , sample output below void myview::zoom(int delta) { double factor = pow(1.2, delta/abs(delta)); this->scale(factor, factor); } void myview::scrollhorizontal(int level) { qpointf center = maptoscene(viewport()->rect().center()); qdebug() << "center: " << center.x() << ", " << center.y(); centeron(qpointf(center.x() - level, center.y())); } void myview::wheelevent(qwheeleve...
Read more