restart; with(Student[NumericalAnalysis]); with(PDEtools, casesplit, declare); x := .5; a := 1; b := .5; d := 1; Q := 2; omega := .4; h1 := 1+a*cos(x); h2 := -d-b*cos(x+omega); F := Q-1-d; de := diff(alpha*f(y), y, y, y, y)+G*(diff(theta(y), y, y))+B*(diff(phi(y), y, y))+6*beta*(diff(f(y), y, y))*(diff(f(y), y, y, y))^2+3*beta*(diff(f(y), y, y, y, y))*(diff(f(y), y, y))^2 = 0, diff(theta(y), y, y)+Nb*(diff(theta(y), y))*(diff(phi(y), y))+Nt*(diff(theta(y), y))^2 = 0, diff(phi(y...

Hi 

I have a big maple sheet on 1000+ lines writing in 2D math in maple 16 that does som calculations on data sets. The sheet can do calculations on 1 data set at the time so I would like to reuse this sheet i a loop. So far I have not found a solution to get a do loop to work over several lines/inputs in a way similar to:

> for i from 1 to 10 do
>   exp1 
>   exp2 

....
>   exp1000+
> od:

Does anyone know how this can be done? 

Program:

restart;

printlevel:=3:

for a from 3 to 30 do

   for b from 23 to 30 do

       for k from b to 30 do

         for d from 4 to 30 do

            for e from 10 to 30 do

  big:=max(sqrt(a^2+b^2), sqrt(k^2+d^2), b+k,e,d+e-a):

  small:=min(sqrt(a^2+b^2), sqrt(k^2+d^2), b+k,e,d+e-a):

I have a somewhat time-consuming numerical integration which is needed for many different parameters so I'd like to parallelize the various integrations. Threads:-Seq should be just the ticket for this; so I write:

respr128:=[Threads:-Seq(evalf(Int(S(ener,Ecrit)*rangeexp(amp,depth/4,t0[128]),ener=0.001..10*Ecrit,method=_d01ajc,epsilon=1.0e-5)),depth=1..36)];

The functions S and rangeexp are defined thus:

S:=(omega,omega_c)->(9*sqrt(3)/(8*Pi)*omega/omega_c...

restart; b := 20; de1 := diff(f(eta), eta, eta, eta)+f(eta)*(diff(f(eta), eta, eta))+1-(diff(f(eta), eta))^2+k*(f(eta)^2*(diff(f(eta), eta))-2*f(eta)*(diff(f(eta), eta))*(diff(f(eta), eta, eta))) = 0, diff(theta(eta), eta, eta)+Pr*f(eta)*(diff(theta(eta), eta)) = 0, f(0) = 0, (D(f))(0) = 1, theta(0) = 1, (D(f))(b) = 0, theta(b) = 0; d1 := subs(Pr = .5, k = 1, [de1]); d2 := subs(Pr = 1, k = 1, [de1]); d3 := subs(Pr = 1.5, k = 1, [de1]); d4 := subs(Pr = 2, k = 1, [de1]); da1 := dsolve(d1, numeric...

Hi maplers

I was wondering if it is possible to change default document template and startup code?

/Delkrak

I am facing the problem below

                 i have a sequence q_n: n=0 to ∞

                 another sequence is defined as p_n=q_{n-1 ,} n=1 to ∞

                  ...

How to construct the decomposition-field of the polynomial x^4 + x + 2 over GF(11,1), where GF(11,1)
is Z[11] = Z mod 11 (see ?GF for more info), in Maple?

restart;
b := 7; de1 := diff(f(eta), eta, eta, eta)-M^2*(diff(f(eta), eta))-(diff(f(eta), eta))^2+m1*f(eta)*(diff(f(eta), eta, eta)) = 0, diff(theta(eta), eta, eta)+m1*Pr*f(eta)*(diff(theta(eta), eta)) = 0, f(0) = S, (D(f))(0) = -1, (D(theta))(0) = -1, (D(f))(b) = 0, theta(b) = 0; d1 := subs(m1 = 1, S = 1.75, Pr = 0.7, M = 2, [de1]); d2 := subs(m1 = 1, S = 1.80, Pr = 0.7, M = 2, [de1]); d3 := subs(m1 = 1, S = 1.85, Pr = 0.7, M = 2, [de1]); d4 := subs(m1 = 1, S = 1.90, Pr = 0.7, M = 2, [de1...

hi

I have two sets of data namely x[i] and y[i], at first i want to pass a spline from all points (points are in 2d x-y diagram) and then extract some point from that. how can i write the suitable code  in maple? 

thanks

Hi There,

Wondering if anyone was able to assist...

I'm trying to differentiate the following equation with respect to Qin:

f(T,Qin,Qp1,Qp2)=T(Qin-Qp1)*(Qp2-Qin)/(Qp2-Qp1)

I have tried using the diff command but with no luck, Maple is returning D(T)(Qin-Qp1)*(Qp2-Qin)/(Qp2-Qp1)-T(Qin-Qp1)/(Qp2-Qp1) which is not complete. 

Any assistance would be much appreciated.

http://www.file-extensions.org/imgs/app-picture/3023/maple.gif

I want to make an application like the picture in Maple 16. Cuold you help me how to make it ?

Do you have any applications about physics nn maple 16?

ok I want to extract data from http://online.wsj.com/mdc/public/page/2_3024-AsianStocks_HONG_KONG.html
So in maple I type:

HTTP[Get]("http://online.wsj.com/mdc/public/page/2 _3024-AsianStocks_HONG_KONG.html", timeout = 100)

but this only give me a bunch of crap and no data and is extremly slow.

What can I do to get the data in a much more efficient way?

hi

i have two row data in excell and i import in maple as the following

Cl_alpha:=<<0.0887574|0.35157>,<0.502959|0.401794>,<0.91716|0.452018>,<1.30178|0.502242>,<1.71598|0.552466>,<2.13018|0.602691>,<2.54438|0.652915>,<2.95858|0.703139>,<3.37278|0.753363>,<3.78698|0.803587>,<4.20118|0.853812>,<4.5858|0.900448>,<5|0.950673>,<5.4142|1.0009>,<5.8284|1.05112>,<6.21302|1.09776>,<6.62722|1.14439>,<7.01183|1.19462>,<7.42604|1.24126>,<7.84024|1.28789>,<8.25444|1.33453>,<8.63905|1.37758>,<9.05325|1.41345>,<9.46746|1.44574>,<9.88166|1.47444>,<10.2959|1.49596>,<10.7101|1.50673>,<11.0947|1.50314>,<11.5089|1.47803>,<11.8935|1.43498>,<12.2781|1.42422>,<12.6923|1.42063>,<13.0769|1.41704>,<13.4615|1.40628>,<13.8757|1.39552>,<14.2604|1.38117>,<14.645|1.38117>,<15.0592|1.3704>,<15.4438|1.35605>,<15.858|1.33812>,<16.2722|1.32735>,<16.6864|1.29148>,<17.071|1.29507>,<17.4556|1.27354>,<17.8698|1.19462>,<18.2544|1.23767>,<18.6391|1.18027>,<19.0533|1.24843>,<19.4675|1.26996>,<19.8521|1.21973>,<20.2367|1.23049>,<20.6509|1.24843>>;

where first row is x and second row is y. i want to interpolate the y values for new x values at the following

xnew:=<<0.119048>,<0.505952>,<0.922619>,<1.30952>,<1.72619>,<2.14286>,<2.52976>,<2.94643>,<3.3631>,<3.77976>,<4.16667>,<4.58333>,<4.97024>,<5.3869>,<5.77381>,<6.19048>,<6.54762>,<6.99405>,<7.38095>,<7.79762>,<8.18452>,<8.60119>,<8.9881>,<9.40476>,<9.82143>,<10.2381>,<10.625>,<11.0417>,<11.4286>,<11.8155>,<12.2024>,<12.5893>,<13.006>,<13.3631>,<13.7798>,<14.1667>,<14.5536>,<14.9405>,<15.3571>,<15.744>,<16.1607>,<16.5774>,<16.9643>,<17.3512>,<17.7381>,<18.125>,<18.5417>,<18.9286>,<19.7024>,<19.3155>,<20.119>,<20.506>>;

please help me

thanks for your attention in advance

I thought the assume command would force Maple to not return solutions that violated the assumption.  But it does not.  How do I make Maple only return those values according to an assumption?

Here is what I did: