how to get numerical values of the following:

cos(pi/2)^2  where square should be on cos not on angle

 i need a help on how to find the local truncation error of a hybird block method of three order IVP. here is the one of the discrete methods i need to find the local truncation error. please any solution should put here thanks 

I want to obtain a report in a particular table format in TextArea component.REPORT (TextArea)area and ComboBox are not visble here.

Can any one help me? TextBox_Format_Output_Request.mw
 

Download TextBox_Format_Output_Request.mw

I attach the doc for clarity of my odoubt. I have made all the required compoennets commands except for the formatting of printf formaat. Thanks for help.

Ramakrishnan V

Hi

trying to make table of points of the pair(lambda1,lambda2) to get the approximate function this technique  is called interpolation

contourplot(4*lambda2*result^2/(Pi*(lambda2+1)^2)-lambda1,lambda1=0..1,lambda2=0..1,contours=[0],axes=boxed,title=tit,titlefont=[SYMBOL,16],thickness=1,color=black,font=[1,1,18],tickmarks=[2, 4],linestyle=1,view=[0.002..1,0.002..1]);

any one knows how to find it

Hello everyone,

I need your helps .I have this code to solve equation numerically with random coefficients, but sometimes the solution find a singularity, for that I create conditional if to avoid it , but its does not work.

Thank you a lot.

Download g.mw

Hi all,

I am hoping to use Maple to develop a list of 4 distinct linear factors with the following conditions...

1)  each individual factor does not have a common integer factor...  ie...  3x+9 is not permitted

2)  the coefficients in each individual factor has an integer between (-9 and 9) without the possibility of zero...  ie x + 0 and 0*x +3 are not permitted

3) each seperate factor is prime to one-another.  ie.  factor1=3x+4  and factor2=-3x-4 is not permitted.

I was hoping to use something short and sweet, using a sequence (if possible), but continue to run into problems.

Any help out there would be fantastic

Thanks,

Mark

Would someone be able to explain this strange output?

restart;
sol:=int(1/(x*sin(x)),x);
algsubs(x=t,sol);

How did Maple manage to obtain this? I tried tracing but it did not help me figure it out.  If it is of any help:

int(1/sin(x),x);

It seems algsubs got confused somewhere?

Maple 2017.2, windows.

I am learning how to use Threads. I looked for an example code in Start / Programming / Example: Task Model and Multithreaded Programming. I found a MandelBrot code, which I was able to change a value in the If -then statement to run in Threads or no Threads mode. I added a time() function to get a sense of the execution time. The results are about 140s with Threads and 91s with no Threads. I would have thought that the execution time would be smaller in the Threads mode than in the no Threads mode. Did I assume incorrectly that I was creating a Threads and no Threads situations by modifying the if-then statement? Or is this not good example? Here is the code:

Threads_Ex3a.mw

with(PDEtools); declare(u(x, y, z, t), U(X, Y, Z, T)); interface(showassumed = 0); assume(a > 0, p > 0); W := diff_table(u(x, y, z, t)); E := 6*W[]*W[x]+W[t]+W[x, y, z] = 0; InvE := proc (PDE) local Eq1, Eq2, Eq3, Eq4, tr1, tr2, tr3, tr4, term1, term2, term3, term4, sys1; tr1 := {t = T/a^beta, x = X/a^alpha, y = Y/a^mu, z = Z/a^nu, u(x, y, z, t) = U(X, Y, Z, T)/a^zeta}; tr2 := eval(tr1, zeta = 1); Eq1 := combine(dchange(tr2, PDE, [X, Y, Z, T, U])); Eq2 := map(lhs, PDE = Eq1); term1 := select(has, select(has, select(has, rhs(Eq2), a), beta), a); term2 := expand(rhs(Eq2)/term1); term3 := select(has, select(has, term2, a), a); sys1 := {select(has, op(1, term3), a) = 1, select(has, op(2, term3), a) = 1}; tr3 := solve(sys1, {alpha, beta, mu, nu}); tr4 := subs(tr3, tr2); print(tr3, tr4); Eq3 := dchange(tr4, PDE, [X, Y, Z, T, U]); term4 := select(has, op(1, lhs(Eq3)), a); Eq4 := expand(Eq3/term4); PDE = simplify(Eq4) end proc; InvE(E)

> restart;
> u := c-6*mu*(1+lambda*sqrt(-mu)*coth(A+sqrt(-mu)*(x+y+mu*t)))/((a+b)*sqrt(-mu)*coth(A+sqrt(-mu)*(x+y+mu*t)))+6*mu*sqrt(sigma*(1-coth(A+sqrt(-mu)*(x+y+mu*t))^2))/(sqrt(sigma)*(a+b)*sqrt(-mu)*coth(A+sqrt(-mu)*(x+y+mu*t)));
            /                (1/2)     /         (1/2)               \\
       6 mu \1 + lambda (-mu)      coth\A + (-mu)      (x + y + mu t)//
   c - ----------------------------------------------------------------
                         (1/2)     /         (1/2)               \     
            (a + b) (-mu)      coth\A + (-mu)      (x + y + mu t)/     

                                                                 (1/2)   
               /      /                                       2\\        
               |      |        /         (1/2)               \ ||        
          6 mu \sigma \1 - coth\A + (-mu)      (x + y + mu t)/ //        
      + -----------------------------------------------------------------
             (1/2)              (1/2)     /         (1/2)               \
        sigma      (a + b) (-mu)      coth\A + (-mu)      (x + y + mu t)/
> c := 1;
                                      1
> mu := -1.5;
                                     1.5
> lambda := 1.5;
                                     1.5
> a := 1;
                                      1
> b := 1;
                                      1
> A := .5;
                                     0.5
> y := 0;
                                      0
> plot3d([abs(u)], x = -10 .. 3, t = -10 .. 3);
Warning, unable to evaluate the function to numeric values in the region; see the plotting command's help page to ensure the calling sequence is correct

> restart;
> u := a[0]-12*sqrt(mu)*(1+lambda*sqrt(mu)*tan(A-sqrt(mu)*(x+y+4*mu*t)))/((a+b)*tan(A-sqrt(mu)*(x+y+4*mu*t)));
              (1/2) /             (1/2)    /      (1/2)                 \\
         12 mu      \1 + lambda mu      tan\A - mu      (x + y + 4 mu t)//
  a[0] - -----------------------------------------------------------------
                                /      (1/2)                 \            
                     (a + b) tan\A - mu      (x + y + 4 mu t)/            
> a[0] := 2;
                                      2
> mu := 1.5;
                                     1.5
> lambda := 2;
                                      2
> a := .5;
                                     0.5
> b := .5;
                                     0.5
> A := 1.5;
                                     1.5
> y := 0;
                                      0
> plot3d([abs(u)], x = -3 .. 3, t = -3 .. 3);
Warning, unable to evaluate the function to numeric values in the region; see the plotting command's help page to ensure the calling sequence is correct

i have a problem with readstat, when using proc , the readstat command returns some errors and do not complete the reading procedure, but without using proc, it is ok, what should i do ? is this a bug?
 

Download read.mw

Hi. I would like to learn to use your package. The example  document for it not running (for me). Included a screen shot and my test document. Basicially It would like to import .mpl procedures and use them as the exports of a module/package. I have no experience with code editors and this is the first time I have tried to use the code edit region. Also whathat does the (**) mean/do?

Edit: I redid the documet.

Download 1-CodeBuilder_Pkg_test.mw

@Joe Riel 

Hello,

I am plotting a three dimensional plot within a range in maple. There are two variables x and y. Third dimension is the objective function ( Z).

On this three dimensional plot, I want to show the point ( preferably by some colour, say red) which have maximum objective value. How it can be done in maple?

Also, Can I show the maximum objective function value (Z) and its coordinates (x,y) on the graph.

Thanks and regard,

Nilesh

I want ot add an annotation to a plot using the drawing facility. Then I export the plot. My problem is that I cannot increase the viewport of the plot without also increasing the size of the plot. For example the code:

produces this plot:

I have constrained the plot and used the size option which together have the side effect of giving me more horizontal space. But this is a wretched hack! I want to size the viewport of the plot so I can write on it.