Recently, my post is deleted without my knowledge.

https://www.mapleprimes.com/posts/209617-How-To-Plot-Stream-Line-contour?sp=209617

I have 2 same integrals, but there is difference between  them at the end.WHy?

Download Ajib.mw

I am trying to solve a set of differential equations in maple. The maple code:

de1:=diff(u1(x,t),x)=-h*(u1(x,t))+i*exp(-i*t)*(u2(x,t));
              d
       de1 := -- u1(x, t) = -h u1(x, t) + i exp(-i t) u2(x, t)
              dx

> de2:=diff(u1(x,t),t)=(-h-i/2)*u1(x,t)-h*exp(-i*t)*u2(x,t);

         d
  de2 := -- u1(x, t) = (-h - 1/2 i) u1(x, t) - h exp(-i t) u2(x, t)
         dt

> dsolve({de1,de2},{u1(x,t),u1(x,t)});
 

However, when I execute it, maple show "Warning: system is consistent". How to remove it

I have a matrix 50*500000 was computed through a loop. when I double click on it I faced with the following error.

Rtable browser error

unable to display structured data

What the problem is?

I asked a question about importing a matrix with symbolic entries the other day.  Just got email saying there had been a reply.  But I can't find my question or that reply!  ............Edgar

Hi,

I would like to translate a set of ODEs from Maple to FORTRAN and optimize them. The problem is that even I have all the equations in one file the maple uses some names in optimizations same as what it has used before. How can I fixed it?

for Example:

with(codegen, optimize, makeproc, cost, prep2trans):

Fortran(eq_dP_5, optimize, resultname = fp1);
      t2 = int(z5 ** 2)
      t3 = int(t2 ** 2)
      t5 = int(z4 ** 2)
      t6 = int(x ** 2)
      t8 = x - 1
      t9 = t8 ** 2
      t10 = 1 / t9
      t11 = t10 * t6 * t5
      t13 = int(mTh ** 2)
      t23 = 2 * z2 * x + z4 * x - 2 * z2
      t26 = 1 / z4
      t32 = Az5 ** 2
      t38 = int(z6 ** 2)
      fp1 = -dble(z4) / dble(t23) / dble(t8) * (dble(t11 * t3 * LmdSF) +
     # dble(2 * t11 * t13 * t2) + dble(2 * t26 / Az5 * t23 * Bz5 * z6 *
     #t8) - 0.16D2 * dble(t10) * t32 * z1 * dble(t6) * dble(t5) * 0.3141
     #592654D1 - dble(2 * t26 * t38 * t23 * t9 * x) + dble(4 / x * t26 *
     # t8 * z2)) * (z1 + valrho) / 0.4D1

Fortran(eq_dmu_5, optimize, resultname = fmu2);
      t1 = z6 ** 2
      t2 = z4 / 0.2D1
      t5 = x ** 2
      t6 = t5 ** 2
      t28 = z5 ** 2
      t29 = t28 ** 2
      t32 = Az5 ** 2
      t35 = mTh ** 2
      t39 = z4 ** 2
      t49 = (x - 1) ** 2
      t50 = t49 ** 2
      fmu2 = (t1 * (z2 + t2) * t6 * dble(x) - 0.5D1 * t1 * (z2 + 0.2D1 /
     # 0.5D1 * z4) * t6 + 0.10D2 * (z2 + 0.3D1 / 0.10D2 * z4) * t1 * t5
     #* dble(x) - 0.10D2 * (z2 + z4 / 0.5D1) * t1 * t5 + ((0.5D1 * z2 +
     #t2) * t1 + 0.4D1 * (LmdSF * t29 / 0.16D2 + t32 * 0.3141592654D1 *
     #valrho + t35 * t28 / 0.8D1) * t39 * z4) * dble(x) - z2 * t1) * dbl
     #e(x) / dble(t50)

I get this message everytime I am trying to open my worksheet: "Recover Corrupt File and Save As"

2.G_anden_aflevering.mw

Is there a simple (i.e. without using of op( ), coeff( ), etc.)  way to force Maple to put a common constant out of brackets in the expression like 2*a+2*b -> 2(a+b)? The task seems as very simple. But I can't find a solution

I've tried alot, but nothing seems to work. The problem occurres whenever im trying to use multiple SI Units in caculations that require division. Example:

Download Maple_prob1.mw

Hello!

I have a pretty simple problem and I feel pretty dumb just asking this. 
If I have a 2x2 Matrix and I remove one column and one row, how do I get the result to be a scalar, not 0th vector?

R := RandomMatrix(2);
R1 := DeleteRow(R,-1);
R2 := DeleteColumn(R1,-1);

                         R := [ -82  41 ]
                                  [-70  91]
                        R1 := [-82  41]
                          R2 := [-82]

I am trying to find the positive integer numbers a, b, c, d, e, f so that the function y =  (a x^2 + b x + c)/(d x^2 + e x + f) is increasing function and the equation  (a x^2 + b x + c)/(d x^2 + e x + f) = x has three integer solutions. I tried
restart;
ListTools[Categorize];
L := [];
for a to 10 do
for b to 10 do
for c to 10 do
for d to 10 do
for e to 10 do
for f to 10 do
if a*e-b*d > 0 and a^2*f^2-a*b*e*f-2*a*c*d*f+a*c*e^2+b^2*d*f-b*c*d*e+c^2*d^2 < 0 and igcd(a, b, c) = 1 and igcd(d, e, f) = 1 and igcd(a, b, c, d, e, f) = 1 and x = (a*x^2+b*x+c)/(d*x^2+e*x+f) then L := [op(L), [a, b, c, d, e, f, x]]

L:=[op(L), [a, b, c, d,e,f,x]]; fi; 
od: od: od: od:od:od: 
nops(L);  

I did not get the result. 
With Mathematica, I got like this picture (a, b, c, d, e, f, x1, x2, x3)

How can I get the correct results?

Is it possible to save table settings (and set a new default), so I dont have to manually change the "indent" number in every table i make in a document?

Its not possible to save the settings in the table menu itself, and there is no "table" option under preferences.

Almost Happy maple user on macbook pro 18

Here I put two functions defined by piece-wise command and then using int(...,numeric) to take their integration but Maple computes only one while for the other one just gives the expression back. The integration domain is bounded, the function is bounded too so I don't think the reason is being divergent in the one that Maple is not computing it.

Is there any prepared command or package for Maple which computes the integrals like this or is better than just int(...,numeric)? Or I have to write an algorithm an force Maple to compute my integral using an algorithm given by myself?
--
The codes in plain text:
 

int(piecewise(0 < -k5*x1*x2+k1*x1 < 1, 1, -k5*x1*x2+k1*x1 <= 0 and 1 <= -k5*x1*x2+k1*x1, 0)*piecewise(0 < -k5*x1*x2+k2*x2 < 1, 1, -k5*x1*x2+k2*x2 <= 0 and 1 <= -k5*x1*x2+k2*x2, 0), k1 = 0 .. 1, k2 = 0 .. 1, k5 = 0 .. 1, x1 = 0 .. 1, x2 = 0 .. 1, numeric);

int(piecewise(0 < x+y < 1, 1, x+y <= 0 and 1 <= x+y, 0)*piecewise(0 < y*x^2 < 1, 1, y*x^2 <= 0 and 1 <= y*x^2, 0), x = 0 .. 1, y = 0 .. 1, numeric);

I have a file consisting of

a b

b a

I want to read it into Maple and convert to a 2x2 matrix.  This will be trivial to most of you, but it isn't to me!  Thanks for some help.  .........Edgar

Hi

I need your help, to find all possible leading behavior of the following curves as x ->infinity of the following curves  defined using the implicit equation:

restart;
with(plots);
implicitplot(x^2*y^2-2*x*y^3+y^4-y^3+x^2 = 1, x = 100 .. 1000, y = 0 .. 1000);
 

Maybe we must change the range of y or x to conclude and obtain all possibles behavior.

I hope have any idea to conclude.

many thanks