Hi,

I get an error when I try to create a modelica block including a piecewise function with one of its expressions including csgn function.

similarly when I use the piecewise And for its conditions.

are these not supported in Modelica?

Hello,
When I try to put the Gcdex in a procedure and start maplemint, then there occurs an error.

Gcdex(x^2 - 1, x - 2, x ,'s','t') mod 3;

--> works

But:

restart;
a := proc()
    Gcdex(x^2 - 1, x - 2, x ,'s','t') mod 3;
end proc;
maplemint(a);

Then there is an error I don't understand.

Error, (in maplemint/expression) not implemented POLY

By the way I have a fundamental problem to understand, where the values s and t are saved after calling Gcdex (or Quo, Rem, etc.). Till now I thought, that variables s and t are created, but when I declare s, t at the beginning as local variables and start maplemint, then there is something like:

    These parameters have the same name as constants:
      3
    These local variables were used before they were assigned a value:
      r::name, (-x-1)::name, (x-1)::name

So the names of s and t changed, they don't assign a new value? I don't understand that.

How do I make find and replace work?  Currently the replace and find button is grayed out.  What magic gets me into a state where the button can be used?

Thanks

P.S. Is there any "package" or "mode" or way some how that emacs key bindings can be made to work (including things like find and replace)?  The user interface would be much improved if I knew how to enable that.

I having a hard time with defining a vector, in order to store in it some data, then plot it and export it to a file, I copied all what's in the help instructures but it doesn't work everytime, please it's urgent for my PhD thesis !

Hello dear Maple,

My name is Bulat, I'm student of Kazan National Research Technical University ( Russia). In our High Program we used your product ( Maple V, Release 4). Now I have two problems and I haven't no idea how I resolve their. I am forced to ask for your help. I upload PrintScreen of my two problems. Please help me to solve them. I' ll be grateful for your help. Sorry for my English :(.

Yours very truly, Bulat

I  encountered a non-integrable integral in the process of solving the following process, . How to achieve its numerical solution? Such as in a looping   code：

#######
pa[i] := pa[i-1]-(Int(subs(t = tau, Lpa[i-1]+Na1[i-1]-Na2[i-1]), tau = 0 .. t)); 

pw[i] := pw[i-1]-(Int(subs(t = tau, Lpw[i-1]+Nw1[i-1]-Nw2[i-1]), tau = 0 .. t)); u[i] := u[i-1]-(Int(subs(t = tau, Lu[i-1]+Nu1[i-1]+Nu2[i-1]), tau = 0 .. t))；

######
Detailed code see annexBC2.mw

what does it means and what it will do. Can some one help me for solving this

Shootlib := "C:/Shoot9"; libname := Shootlib, libname; with(Shoot);

while i m receiving the following message:

"Error, invalid input: with expects its 1st argument, pname, to be of type {`module`, package}, but received Shoot "

Full program is :

restart; Shootlib := "C:/Shoot9"; libname := Shootlib, libname; with(Shoot);
with(plots):
N1 := 1.0; N2 := 2.0; N3 := .5; Bt := 6; Re_m := N1*Bt; gamma1 := 1;
FNS := {f(eta), fp(eta), fpp(eta), g(eta), gp(eta), m(eta), mp(eta), n(eta), np(eta), fppp(eta)};
ODE := {diff(f(eta), eta) = fp(eta),
        diff(fp(eta), eta) = fpp(eta),
        diff(fpp(eta), eta) = fppp(eta),
        diff(g(eta), eta) = gp(eta),
        diff(gp(eta), eta) = N1*(2.*g(eta)+(eta-2.*f(eta))*gp(eta)+2.*g(eta)*fp(eta)+2.*N2*N3*(m(eta)*np(eta)-n(eta)*mp(eta))),
        diff(m(eta), eta) = mp(eta),
        diff(mp(eta), eta) = Re_m*(m(eta)+(eta-2.*f(eta))*mp(eta)+2.*m(eta)*fp(eta)),
        diff(n(eta), eta) = np(eta),
        diff(np(eta), eta) = Re_m*(2.*n(eta)+(eta-2.*f(eta))*np(eta)+2.*N2/N3*m(eta)*gp(eta)),
        diff(fppp(eta), eta) = N1*(3.*fpp(eta)+(eta-2.*f(eta))*fppp(eta)-2.*N2*N2*m(eta)*(diff(mp(eta), eta)))
       }:
   
blt := 1.0;
IC := { f(0) = 0,
        fp(0) = 0,
        fpp(0) = alpha1,
        g(0) = 1,
        gp(0) = beta1,
        m(0) = 0,
        mp(0) = beta2,
        n(0) = 0,
        np(0) = beta3,
        fppp(0) = alpha2
      };
BC := { f(blt) = .5,
        fp(blt) = 0,
        g(blt) = 0,
        m(blt) = 1,
        n(blt) = 1};
infolevel[shoot] := 1;
 

HI

please help me for dsolve this nonlinear differential equations

thanls...

HAB.mw
 

Download HAB.mw

I uploaded this file AnimationTest.mw to the Maple Cloud. One a 3d animation and one just a 3d plot. When I open the file in the Maple Cloud from my browser,  I cannot rotate either plot with the mouse. Is there a way to change the worksheets so rotation from the browser will be possible?

I Could Not Write An If Then Or Ifelse Statement. Please Help Me.

f := unapply(x^2-2, x); a := 1; b := 2; n := 10; Digits := 10;
      2    
x -> x  - 2
                               1
                               2
                               10
                               10
c := evalf(eval((a*f(b)-b*f(a))/(f(b)-f(a))));
                          1.333333333
if  f(c)*f(a)<0 then ;
          "      k:=evalf(eval(|(f(c))/(b-c)|)) and "

                          /(1 + k) a f(b) - b f(a)\
             x[i] := evalf|-----------------------|
                          \  (1 + k) f(b) - f(a)  /
            "     elif f(x[i])*f(a)<0 then b:=x[i]"
                 "     else b:=c and a:=x[i] "
                  "     if f(c)*f(a)>0 then "
                 "      k:=|(f(c))/(b-c)|and "

                          /a f(b) - b f(a) (1 + k)\
             x[i] := evalf|-----------------------|
                          \  f(b) - f(a) (1 + k)  /
            "     elif f(x[i])*f(a)>0 then a:=x[i]"
              "     else a:=c and b:=x[i] end if"

Error, unterminated 'if' statement
     Typesetting:-mambiguous(Typesetting:-mambiguous(

       if fApplyFunction(c)sdotfApplyFunction(a)lt0 then , 

       Typesetting:-merror("unterminated 'if' statement")))

Im trying to solve 12 equations with 12 variables but I can't solve. Please help and advise me to solve this problem. Iproject3.mw
project3.mw

Hi, I have a big system with 27 polynomial equations in 16 unknowns: f_1=...=f_27=0.  I can store these equations but I cannot calculate a Grobner basis of the ideal  J generated by my polynomials (allocation problem) - I use the library "with(FGb)"-  What interests me is whether my system is minimal in the following sense.

If, for example,  I remove f_1, is the ideal generated by (f_2,...f_27)  J again ? That is to say, is f_1 in the ideal generated by f_2,...,f_27 ? I would like to get an answer "yes" or "no" for each removed  f_i.

My question: can we solve the problem above  without calculating a Grobner basis of J?

Thanks in advance.

i have a system with 5 dif equations and five unknows. i have told to maple to solve it numerically with interactively solve comand (right cilck button). the window open like it normally does and i put values to my parameters, with an initial condition for the system (Q(0)=0). then i press numerically solve and that's all, the program just keep evaluating with no answer. i wait for 15 min, which i think is too much time, and got any answer yet.

hope you can help with this

thanks.. 
 

Download propuesta_transfer.mw

I assigned

before an algebraic calculation so I would like to get  or have the program print the 70 digits of the answer and not just 10 digits. Because when I press ENTER, I get only 10 digits.

I am trying to implement Subresultant p.r.s. algorithm for calculating greatest common divisor. The algorithm decribed in the book:

My code return the correct GCD, however the sub-resultant terms are different from the result of the built-in function. The last term a[i-1] is huge and involves fractions. I think my implementation is same as the algorithm described in the textbook.

I have attached the file. Could anybody spot anything wrong in the code? Why do fractions still appear? In my code, "lsr" is last subresultant term returned from the built-in function, the second one is my result.

Download subresultant.mw