@roman_pearce thanks for trying. Helps me in deciding whether to focus on simplifying the problem or recoding it elsewhere (prob wont do this now)... Cheers, R

@roman_pearce appreciate you taking the time to ponder this...

All constants are set apart from 'a' and 'b'. They are set for exact values greater than 0.01 smaller than 10. I keep the significant digits to 5 for the sake of getting the code to finish.

In other words the solutions are to be be expressed in terms of 'a's and 'b's.

Great info on techniques, thanks.

q:=0.1; be1R:=2; be2R:=8; nu1R:=0.1; nu2R:=0.1; al1:=1.1; al2:=5; sig1:=1; sig2:=1; gamma2R:=1.1; 

gamma1R:=1.1*theta1R;

...with theta1R:=1+2.5*(1-(a/1.9)^4);