I want to fetch the trigger time.
I use rhs function.
 

Download RHS_TEST.mw

I want to find all complex numbers such that 
abs(z)*(z-4-I)+2*I = (5-I)*z.
I tried solve(abs(z)*(z-4-I)+2*I = (5-I)*z, z)
and got the answer is z = -1.
This is a question in a test with multiple choice. The key of question is three numbers.

Is my commant wrong? Where is wrong?

If I have an expression which contains an unknown function f(t) e.g.

ex:=1+exp(f(t))*f(t+1)+f(t/2)

how do I replace it with some defined function F(t)?

Do I need to work with ::anything, i.e. something like

subs(f(t::anything)=F(t::anything),ex)

?

However this doesn't work.

Download Rindler_simulation_v2.mw

Dear Maple primes,

Could you, please, help me with numerical solution of an ODE?

The ODE looks like this

dz/dx = f₁(x,z) + f₂(z)

where f₁(x,z) is some simple function of x and z (that does not create any problem), but f₂(z) is given as

f₂(z) = int(f₃(t), t = z₁..z₂)

The problem appears, when the integral cannot be solved analytically.

Below is an example of the problem (here I chose the function f₃(t)= t^t as well as other functions, intervals and initial condition only for the sake of illustration of the problem):

Download z-for_primes_dsolve.mw

Thank you in advance!

Max

I really love this feature of course because it just allows you to learn so much about function properties via comparative learning processes, but I would also love to have the ability to:

1) Add elements to the list that is stored under known_functions:

FunctionAdvisor(known_functions);

2) Add my own categories of information such as the ones that exist like DE, branch_cuts, etc

I want to know about some algorithms used by Maple, if there is a way to dig through the code.

Hello!

I am currently trying to plot the function: sqrt(x^2+y^2-c), where x & y are variables in cartesian coordinates, and c is an arbitrary real constant. I came across a jagged edged plot when I asked maple to plot:

sqrt(x^2+y^2-c)       *** Code below ***

Explore(plot3d([sqrt(x^2+y^2-c),-sqrt(x^2+y^2-c)],x=-2..2,y=-2..2,numpoints=10000),c=-1..1)

So I entered in the "Re" command to take the real part:

Explore(plot3d([Re(sqrt(x^2+y^2-c)),Re(-sqrt(x^2+y^2-c))],x=-2..2,y=-2..2,numpoints=10000),c=-1..1)

This smoothed the plot out at the ring x^2+y^2=1, but however also added what looks like a solid disc as well. Is there a way I can plot the function correctly without this disc in the middle?

Thanks!

solve_V_alfa.mws

I'm trying to solve two simultaneous equations for velocity, V, and angle of inclination, alfa, given the range, ranj and maximum height, H.  The solve part works fine, but gives me several unwanted complex solutions.  I was trying to put the solutions in a list so that I could access them.  The text below is rather verbose, but I've left comments in as a possible educational aid for people understanding expression sequences, and sets. 

restart:
#interface(imaginaryunit, errorbreak=0, errorcursor=true, #labelling=false):
#with(SolveTools);

# Investigating sets and expression sequences
#S:={q, w,e}:
S:=7,2,5:
whattype(S);
L4:=convert(S,list): #Doesn't work for expression sequences
L4:=[S]; # this converts the expression sequence to a list, L4
L4[3]; # outputs 3rd element of L4 which is 5
# End of investigation

print(`Next part - solving two equations: Ranj and height H given`);
print(`To find angle of inclination, alfa, and velocity, V`);
g:=9.81:Ranj:=80:H:=15:
L:=solve({4*V^2*cos(alfa)^2-g*Ranj^2/(Ranj*tan(alfa)-2*H), Ranj*g-V*sin(2*alfa)},{V, alfa}): #::real}):
#L:=solve([4*V^2*cos(alfa)^2-g*Ranj^2/(Ranj*tan(alfa)-2*H), #Ranj*g-V*sin(2*alfa)],{V, alfa}): #::real}):
# ::real, alfa::real:
printf("L is of type %a", whattype(L));
#whattype(L);
#L2:=convert(L,list); #gives error message
L2:=[L]:
L2;
print(`Next: To distinguish between real and complex solutions`);
print(`Output in if statement is not working as expected.`);
#nops([L]);
for i from 1 to nops([L]) do
if type(L[i],complex(float)) then
#if type(L[i],nonreal) then
printf("Complex solution for i=%d\n",i);
else
printf("else part of if: i=%d\n",i);
print( L[i]);
end if;
end do; #for

Download solve_V_alfa.txt

Hi, having two different arrays of the same length,  is there a command to associate every element of one to the other of the other array? Thanks



It means, to have an array compsed by ex := [[2, 3.584328599],[3, ...] ]
Attention : there's no correlation among elements!
Thanks a lot

I'm attempting to calculate eigenvalues of a 16 x 16 matrix (Vibronic_Model_Worksheet.mw) which contains a single variable "x".

However, this matrix contains 24 instances of imaginary numbers.

When I try to solve for any number other than x = 0, I get: "Error, (in LinearAlgebra:-Eigenvalues) expecting either Matrices of rationals, rational functions, radical functions, algebraic numbers, or algebraic functions, or Matrices of complex(numeric) values"

Eventually I want to be able to solve for multiple values of x simultaneously.

I'm still fairly new to maple, so I may have overlooked something.

Hello
What is the Problem

Hi,

Here is a set of simple instructions to illustrate the first problem :

restart:
fichier := TheNameOfSomeFile:
writeto(fichier):
printf("Hello\n"):
writeto(termional):

When executed from Maple 2015 or Maple 2016, the content of TheNameOfSomeFile is

Hello

When executed from Maple 2018 it becomes

> printf("Hello\n"):
Hello
> writeto(terminal):

Using interface(echo=0) changes nothing to the result obtained with Maple 2018.

Has anyone ever encountered this problem?
How can I fix it?


Thanks in advance

Download 418.mw