I use WINDOW 10 

for j from 1 to 3 do
pout := cat("C:/Users/Eli/Documents/Animation/", "file", j, ".bmp"):
print(pout):
plotsetup(bmp, plotoutput = pout):
# plot (...)
od:

Hi, I have a problem related to the cancellation of the square root with a square in its argument, see the image below.

Thanks in advance,

Santiago

i have a task who is like this:

let

and the function f meets the following criteria

and f '(-1)=2

what is h'(-1)

how to i solve this, i am thinking of ordinary differential equations (ODEs)

and dsolve funktion in maple. but how do i do this.

or is there another way to solve it

 Hello , I need to plot the following on the same axes of implicit plot where discrete wavenumbers k1, k2, k3 on the same set of axes (the wavenumbers are points [k1,0], [k2,0], [k3,0], ...) , where k_{n}=2*n*Pi/L and the domain [0,L] 

Could you please help , I couldn't do. Thank you

restart;
with(plots, implicitplot);
alpha := 0.1;
beta := 0.1;
mu := 0.5;
u := 0.5;
v := 1;
gamma = 0.1;
                          alpha := 0.1

                          beta := 0.1

                           mu := 0.5

                            u := 0.5

                             v := 1

                          gamma = 0.1

Eq := -2*sigma*alpha*beta*mu*u - sigma*alpha*beta^2*u*v - 2.0*mu*alpha*beta^2*u*v - gamma^2*k^2*alpha*beta*u + beta*sigma^2*v + sigma^3 + 2*mu*sigma^2 - alpha*beta*sigma^2*u + sigma*beta^2*u*v + gamma^2*k^2*sigma + 2.0*mu*beta*sigma*v + 2.0*mu*beta^2*u*v;
                                            2  2              2
 Eq := 0.0995 sigma + 0.004500 - 0.005 gamma  k  + 1.095 sigma 

           3        2  2      
    + sigma  + gamma  k  sigma

implicitplot(Eq, k = 0 .. 10, sigma = -0.1 .. 0.1);
 

Please can someone help me out, have been trying for a very long time now on how to differentiate the lines on the graph.......please find attachment below and correct..................................................................

 restart;
  with(plots):
#
# Define the ODE system
#
  odeSys:= { (diff(F(eta), eta, eta, eta))*(1+epsilon-alpha((diff(F(eta), eta, eta))^2))+F(eta)*(diff(F(eta), eta, eta))+S*(diff(F(eta), eta))-(1/2)*S*eta*(diff(F(eta), eta, eta))-(diff(F(eta), eta))^2-M*(diff(F(eta), eta)), (diff(theta(eta), eta, eta))*(1+R)-delta*(diff(F(eta), eta))^2-Pr((3/2)*S*theta(eta)+(1/2)*S*eta*(diff(theta(eta), eta))-2*(diff(F(eta), eta))*theta(eta)+F*(diff(theta(eta), eta)))};
#
# Define the first set of boundary conditions
#
  bcs1:= { F(0) = 0, (D(F))(0) = 1, (D(F))(inf) = 0, theta(0) = 1, theta(inf) = 0
         }:

  RVals:=[0.1, 0.5, 1]:
  for k from 1 by 1 to numelems(RVals) do
      pList:=[ epsilon = 0.18, M = 0.5, S = 1.5, delta = 0.3, Pr = 1.5, alpha = 0.4, R = RVals[k],inf=1
]:
      sol1[k]:= dsolve( eval
                        ( `union`( odeSys, bcs1),
                           pList
                       ),
                       numeric
                     );
     od:
  display
  ( [ seq
      ( odeplot
        ( sol1[i],
          [eta, theta(eta)],
          eta=0..2
        ),
        i=1..numelems(RVals)
      )
    ],
    color = [red, green, blue],
    title = typeset( theta(eta), " versus ", eta),
    titlefont = [times, bold, 20]
  );
 

I want to do a organised introduction in Maple  with the student and main vectorcalculus package
Is there some studymatrial to find for this ?
 

I am using Maple to solve a system of ODEs numerically. Right now, I want to find the integration of the output of the system of ODEs. How it is possible to do this? 

F := dsolve(ODESys union ICs, {y0(t), y1(t), y2(t), y3(t)}, type = numeric)

Y0 := t -> rhs(op(2, F(t)))

Now, I want to find int(Y0,t=0..1).

Why this simple case cannot be simplified to Complete Square form

restart:

EQ := 1/4*n[1]^4 - 1/2*n[1]^2 + 1/4

simplify(1/4*n[1]^4 - 1/2*n[1]^2 + 1/4, 'size')

with(Student[Precalculus])
CompleteSquare(EQ);
                   

Dear all, 

Would you tell me a way to apply the identity relationship in Maple, to rewrite 'subexp' below to 'subexp2'?

Download Q20200901.mw

In the plane, an ABC triangle is considered for the vertices B and C are fixed, A being variable so that b +c remains constant and equal to a given length l. (b=distance(A,C), (c=distance(A,B)
How to show that the product tan(B/2=*tan(C/2) remains constant ?

Hello,

I would like to solve this equation and keep the parameters a and eta unknown and assign them a value later on. Is it possible ?

I know I can solve this equation numerically, it works, but I want to change the parameters to fit experimental data.

I tried this :

ode:=diff(y(t),t)=-(4/3)*(a/eta)*(y(t)^3+1):
lambdaE:=dsolve({ode,y(0)=1},y(t),explicit);

But the solution appears with RootOf expression and I can't get a numerical solution afterwards.

lambdaE := y(t) = RootOf(6*eta*sqrt(3)*arctan(((2*_Z - 1)*sqrt(3))/3) - eta*sqrt(3)*Pi - 3*eta*ln(_Z^2 - _Z + 1) + 6*eta*ln(_Z + 1) - 6*eta*ln(2) + 24*t*a)

I tried some things from forums but nothing works.

Do you have any leads on how to do it please ?

Thank you in advance

Manon

is there easy way to simulate lines rotation like robotic arm following point P moving by adjusting angle q1 and q2?

i feel point two lines and simulate static two lines in animate is also not easy.

with(plots):
with(geometry):
#Joint at (5,5)
#length a  2, y=7, x, solve(subs(y=7,2^2 = (y-5)^2+(x-5)^2), x)[1];
#length b 2, y=1, x, solve(subs(y=4,2^2 = (y-5)^2+(x-5)^2), x)[1];
point(A, 5, 5), point(B, solve(subs(y=7,2^2 = (y-5)^2+(x-5)^2), x)[1], 7):
point(C, 5, 5), point(D, solve(subs(y=4,2^2 = (y-5)^2+(x-5)^2), x)[1], 4):
line(L1, [A, B]);
line(L2, [C, D]);
_EnvHorizontalName := x;
_EnvVerticalName := y;
Equation(L1);
Equation(L2);
solve(Equation(L1), y);
solve(Equation(L2), y);
plot([solve(Equation(L1), y),solve(Equation(L2), y)], x = 1 .. 2);
animate(plot, [Equation(L1),Equation(L2), x = 0 .. 2], A = -3 .. 3);

Angle Reference from below

angle of L2 , q1a

angle of L1,  q2a

p1 := vector([a*cos(q1), a*sin(q1)]);
P := evalm(p1 + b*vector([cos(q1+q2), sin(q1+q2)]));
Pe := map(expand, P);
A := {cos(q1) = c1, sin(q1) =s1, cos(q2)=c2, sin(q2)=s2};
P := subs(A, op(Pe));
F1 := [x - P[1], y - P[2], s1^2+c1^2-1, s2^2+c2^2-1 ];
F2 := subs({a=1, b=1}, F1);
with(Groebner):
g1 := Basis(F2, plex(c2, s2, c1, s1, y, x));
g1 := [-4*x^2+x^4+2*x^2*y^2+y^4-4*y*x^2*s1-4*y^3*s1+4*s1^2*x^2+4*s1^2*y^2, 2*y*s1+2*x*c1-x^2-y^2, 4*x-x^3-y^2*x+2*y*x*s1-2*y^2*c1-4*s1^2*x+4*y*s1*c1, s1^2+c1^2-1, x*s1-y*c1+s2, -y^2+2-x^2+2*c2];
map(indets, g1);
 

g2 := Basis(F2, plex(c2, s2, c1, s1));
 

assume point x = 1, y = 0.5
g2spec := subs({x=1, y=1/2}, [g2[1],g2[3]]);
S1 := [solve(g2spec[1])];
q1a := evalf(arcsin(S1[1]));
q1b := evalf(arcsin(S1[2]));

S2 := [solve(subs(s1=S1[1], g2spec[2])), solve(subs(s1=S1[2], g2spec[2])) ];
q2a := evalf(arcsin(S2[1]));
q2b := evalf(arcsin(S2[2]));

q2a := subs({a=2,b=2}, q2a);
if q2a < 0 then
 q2a_pi := evalf(eval(-q2a-Pi));
else
 q2a_pi := evalf(eval(Pi-q2a));
end if:
q2a_pi;

q2b := subs({a=2,b=2}, q2b);
if q2b < 0 then
 q2b_pi := evalf(eval(-q2b-Pi));
else
 q2b_pi := evalf(eval(Pi-q2b));
end if:
q2b_pi;

I've been using Physics:-Latex in my large program, and checking the output visually. I noticed this small issue on one page.

When asking Physics:-Latex to generate Latex for this expression

It gives

While when using standard latex() command, the result is

Below is the worksheet attached. Using Latest Physics and Maple.
 

THe latex file used to compile the Maple output Latex using TeXLive distribution on Linux is

\documentclass[12pt]{article}
\usepackage{amsmath}
%\usepackage{maplestd2e} %not needed for this example
\begin{document}

%output of Physics:-Latex()
\[
\frac{\left(\int {\rm e}^{-t} f \left(t \right)d t \right) {\rm e}^{t}}{2}+-\frac{1}{2} {\rm e}^{-t} \left(\int {\rm e}^{t} f \left(t \right)d t \right)
\]

%output of maple latex()
\[
{\frac {\int \!{{\rm e}^{-t}}f \left( t \right) \,{\rm d}t{{\rm e}^{t}}}{2}}-{\frac {{{\rm e}^{-t}}\int \!{{\rm e}^{t}}f \left( t \right) \,{\rm d}t}{2}}
\]

\end{document}

Download latex_issue_5.mw

Hello everyone, my computer installed Maple 2016, 2020 and set Maple 2016 as default opening files .mw, .mws. How to change default to Maple 2020 . I tried searching on Google but not found.

Thank you very much.

Hello!

How to calculate the residual sum of squares of a Nonlinear Fitting?

I tried but was unsuccessful.