Hi all,

Hi

Just started using the global optimization toolbox. Not sure if I am doing something wrong or if it is genuinely awful. After about 20mins it produces a fit that is dreadful,an alternative cheaper package (scop) produces a great fit in about 2 mins. Can anyone see what I'm doing wrong?

Hi

Just started using the global optimization toolbox. Not sure if I am doing something wrong or if it is genuinely awful. After about 20mins it produces a fit that is dreadful,an alternative cheaper package (scop) produces a great fit in about 2 mins. Can anyone see what I'm doing wrong?

I have the following nonlinear ODE symmetric in the functions A(r) < - > B(r):

eq:= diff(A(r),r)^2/A(r)^2 + diff(A(r),r)*diff(B(r),r)/A(r)/B(r) + diff(B(r),r)^2/B(r)=k;

The task is to solve for A(r) in terms of B(r) and it goes like this: solve the quadratic equation for A'/A in terms of B'/B (or vise versa cause equation is symmetric) and then integrate.

dsolve(eq, A(r));  produces the required solution for A(r) in terms of integral over B and B'.

dsolve(eq, B(r)); returns nothing.

Should it really take 103 seconds to build a list of 95850 floats on a 2GHz machine?

I am using the data from the attached .m file and using the following code (I am not doing anything here but I will be) on the list, angles1_s which is read in from the .m file

Hello,

I need to plot a solid circle on a line. The problem is that you can see the line through the circle:

with(plots):
with(plottools): 
l := plot(x,color=red,thickness=3): 
d:=disk([0.5, 0.5],0.05,color=black): 
display({l, d},view=[0..1,0..1]);

I think I need something which is called in Microsoft Office products "Bring Forward (or to Front)" and "Send Backward (to Back)"

I suppose it's something very simple but I just cannot find it!

Thank you!

I have this function:

Eig1 := proc (R) options operator, arrow; (1/9*(-27*deltaE^3*R^9-10*M^6+3*sqrt(-189*deltaE^4*R^12*M^4-147*deltaE^2*R^6*M^8-27*M^12+60*deltaE^3*R^9*M^6))^(1/3)-(-deltaE^2*R^6-7/9*M^4)/(-27*deltaE^3*R^9-10*M^6+3*sqrt(-189*deltaE^4*R^12*M^4-147*deltaE^2*R^6*M^8-27*M^12+60*deltaE^3*R^9*M^6))^(1/3)+2/3*deltaE*R^3+2/9*M^2)/R^3 end proc

where:

> deltaE := 0.355 , and

M := sqrt(540000.0000)

Hi,

The eigenvalue command gives me a vector like: [ [5] , [4] , [3] ].

Is there a quick way I could grap the first eigenvalue and make that a function without copying and pasting?

I have the following nonlinear ODE symmetric in the functions A(r) < - > B(r):

eq:= diff(A(r),r)^2/A(r)^2 + diff(A(r),r)*diff(B(r),r)/A(r)/B(r) + diff(B(r),r)^2/B(r)=k;

The task is to solve for A(r) in terms of B(r) and it goes like this: solve the quadratic equation for A'/A in terms of B'/B (or vise versa cause equation is symmetric) and then integrate.

dsolve(eq, A(r));  produces the required solution for A(r) in terms of integral over B and B'.

Hi I used the "search" option of this forum. I have seen old messages about  the continuation of line "#" of fortran77 that is different from the ones in fortran90. I have thousands lines of fortran code created by codegeneration and I can't use "replace all # with &", because # are at the beginning and "&" at the end of the line. It takes a while to substitute all of them by hand every time I generate the fortran subroutine. Do you know if they are going to introduce the continuation of line with & instead #?

thanks

Hello

This works:

zz:= s/sqrt(s) ;

produces sqrt(s);

-------------

but zz2:=s/sqrt(s*m);

produces s/sqrt(s*m)  instead of   sqrt(s/m);

---------------

Also, how do I get sqrt(s) to display as (s)^(1/2) ??

Thanks!

Hi everyone,  I am using dsolve with a stop condition and want to suppress the error message "warning, cannot evaluate the solutin further right of.... stop condition 1 violated#

I believe that _Env_dsolve_nowarnstop:=true  will help, but I can't get it to work.

Any ideas please?  and any general info about environmental variables in general?

Thanks

Judith

Hi,the following is my confused problem

the effective statements 
   restart;
   with(group):with(StringTools):with(Maplets[Examples]):
   a:=GetInput("Input the moves.");#input: IMIM10IM2IM9IM7IM8IM10
   i:=0:b:=Explode(a):
   M10:=[[2, 12, 7, 4, 11, 6, 10, 8, 9, 5, 3]]:
   M2:=[[2, 5, 8, 6, 4, 12, 3, 9, 10, 11, 7]]:
   M8:=[[2, 4, 10, 5, 12, 11, 8, 3, 7, 6, 9]]: