I am slightly confused as I can't apply the seemingly correct function to a sequence. It seems like modp does not like my inverse. But I am not aware of any other way of finding the modular inverse.  

Download example.mw

restart;
alias(u = u(x, y, t), f = f(x, y, t));
                              u, f
u := (f+sqrt(R))*exp(I*R*t);
                    /     (1/2)\           
                    \f + R     / exp(I R t)
pde1 := I*(diff(u, t))+diff(u, x, x)+2*lambda*u*abs(u)*abs(u)-gamma*(diff(u, x, t));
   // d   \                /     (1/2)\             \
 I ||--- f| exp(I R t) + I \f + R     / R exp(I R t)|
   \\ dt  /                                         /

      / d  / d   \\           
    + |--- |--- f|| exp(I R t)
      \ dx \ dx  //           

               /     (1/2)\                           2 
    + 2 lambda \f + R     / exp(I R t) (exp(-Im(R t)))  

               2
   |     (1/2)| 
   |f + R     | 

            // d  / d   \\                / d   \             \
    - gamma ||--- |--- f|| exp(I R t) + I |--- f| R exp(I R t)|
            \\ dx \ dt  //                \ dx  /             /
 

I am using Maple 2023.

When using PDE Numeric Help Side the text sides is not smooth moving when scrolling downwards or upwards.

Any advice what to do????

Kjell

Hi,

I am looking to randomize the coordinates of the three non-aligned points A, B, C to generate random drawings of triangles. Thank you for your insights.

CaractTriangle.mw

The so-called look-and-say sequence is something like:

Rosetta Code gives an example of how to generate such sequences in Maple; however, the code in that example is lengthy. More to the point, compared to the subsequent Mathematica analogue, it is rather inefficient (though more readable, but that's not the point). 
I slightly modify the original code, and now the modified version is significantly faster than the original one: 
 

Download LookAndSay.mw

For the sake of convenience, I paste the original Maple code from the page above into here:

generate_seq := proc(s)
	local times, output, i;
	times := 1;
	output := "";
	for i from 2 to StringTools:-Length(s) do
		if (s[i] <> s[i-1]) then
			output := cat(output, times, s[i-1]);
			times := 1; # re-assign
		else 
			times ++;
		end if;
	end do;
	cat(output, times, s[i - 1]);
end proc:

Look_and_Say :=proc(n)
	local value, i;
	value := "1";
	print(value);
	for i from 2 to n do
		value := generate_seq(value);
		print(value);
	end do;
end proc:

#Test:
Look_and_Say(10);

Below is a modified version by me: 

LookAndSay := proc(m::posint, n::nonnegint := 1, $)::'Vector'(nonnegint);
    uses StringTools;
    local delete::'equationlist' := Repeat~(Explode(ExpandCharacterClass(":digit:")), 2), aux::procedure[nonnegint](posint) := proc(_::posint)::nonnegint;
        options cache;
        description "https://rosettacode.org/wiki/Look-and-say_sequence#Maple";
        if _ = 1 then
            n
        else
            local char::'nonemptystring' := String(thisproc(_ - 1)), temp::list(string)(*, str*);
            try
                temp := RegSplit(sprintf("(?=[%s])(?<=(?!%s)[%s])", Unique(char), Join(select[2](foldr, eval(apply), delete, rcurry(OrMap, Unique(char)), curry(curry, Has)), "|"), Unique(char)), char)
            catch "cannot compile regular expression: repetition-operator operand invalid":
                local rule::'equationlist' := (str -> str = Insert(str, 1, " "))~(subs(delete =~ NULL, Generate(2, Support(char))));
                temp := StringSplit(Subs(rule, char), " ") # to be replaced
            end;
            parse(String(seq('length(str), str[1]', str in temp)))
        fi;
    end;
    forget(aux);
    Vector[column](m, aux, datatype = nonnegint)
end:

The form LookAndSay(m, n); produces the first m terms in a "look-and-say" sequence starting with n (which is by default 1). Unfortunately, even if "the modified version is significantly faster than the original one", it is still much slower (25s versus 2.5s) than the uncompiled Mathematica implementation (on the same modern computer):

So, is there a workaround to evaluate LookAndSay(50, 1): in about three seconds (instead of half a minute) in the first call (without looking up a pre-calculated table) (on the same computer) as well? 

This is old question https://www.mapleprimes.com/questions/208909-Code-For-Integer-Points-On-Sphere. Now I see this question at here 
https://mathematica.stackexchange.com/questions/288956/how-can-i-get-all-squares-on-this-sphere-so-that-its-coordinates-are-integer-num
 My idea is select all diameters which diameters are perpendicularly from the sphere (x-2)^2 + (y-4)^2 + (z-6)^2 = 15^2. How can I tell Maple to do that?

HPMthesis.mw

Advanced Engineering Mathematics with Maple (AEM) by Dr. Lopez is such an art.

Mathematics and Control Theory talks easily in Maple...

Thanks Prof. Lopez. You are the MAN !!

Dr. Ali GÜZEL

how to plot graphs for both methods and comparison of different method values for Diff(f(eta),eta, eta) at eta =0

Download DTM_practice.mw

how I can extract data from plot3d such as x, y, phi2?

thanks

tez-1.mw

Download tez-1.mw
 

Download tez-1.mw

Hello 
I have two questions if you please.
Q1)
Is there any kind of a ruler in the maple worksheet that helps me to Adjust the wideness of the line text? 0r even to  spilt the outputs in two columns 
you know as the ruler in the "Word document ".
Q2) Is there any way"method" to change the output type 
since the default adjustment of the outputs is a type of stretchable image? But I need the outputs of the text type. That may  help me to insert them in a given templet

Thank you

Hello everyone, 

I am trying to solve 6 equations with 6 unknowns in maple. I tried to simplify them and they are shown below. Unfortunately, after using solve I get:

Warning, solutions may have been lost
                            sol := ()

Can you help me with it, please? I tried to solve it with Matlab and got a warning that no explicit solution could be found and had no luck with mathematica either.  Thank you very much!!

restart;

omega_1 := 1;
beta_1 = 0.2;
omega_2 := 0.15;
beta_2 := 0.15;
mu := 0.4;
g := 9.81;
K := 5;
vb := 0.5;
tc := arccos(vb/(A*omega));

Ns := mu*(omega_2^2 - K)*B*(2*omega*tc - sin(2*omega*tc) - pi)/pi;

Nc := 4*mu*(omega_2^2 - K)*x2*sin(omega*tc)/pi + mu*(omega_2^2 - K)*C*(2*omega*tc + 2*sin(2*omega*tc) - pi)/pi;

a0 := -mu*(omega_2^2 - K)*x2*(1 - 2*omega*tc/pi) + 2*mu*(omega_2^2 - K)*C*sin(omega*tc)/pi;

eq1 := -A*omega^2 + A*omega_1^2 - B*K - Ns;

eq2 := 2*A*beta_1*omega*omega_1 - C*K - Nc;

eq3 := -2*C*beta_2*omega*omega_2 - B*omega^2 + B*omega_2^2 - A*K;

eq4 := 2*B*beta_2*omega*omega_2 - C*omega^2 + C*omega_2^2;

eq5 := omega_1^2*x1 - K*x2 - a0;

eq6 := omega_2^2*x2 - K*x1 - g;

sol := solve({eq1, eq2, eq3, eq4, eq5, eq6}, {A, B, C, omega, x1, x2});

Hello!

I have a Maple program for plotting from z plane to w plane under transformation w=1/z

I have lines in the z plane and should have had nice circles in the w plane, but the circles are not so good. It is a hole around zero.I can do it manually, but I like to do it with the Maple program.

What must be done in the program to get nice circles???

Thanks for any answers!!

Regards

Kjell

Download convert-latex.mw

Hi everybody

I have a differential equation where derivatives (for example d/dr) have been displayed by prime notation.

I want to convert this equation latex format by keeping the prime notation not d/dr.

I would appreciate it if anyone could help me.

I have attached a sample code.

Thank you

Hadi

I'm having problems trying to plot a discrete plot for a Poisson Distribution CDF.  I followed a model found elsewhere that someone had used to do a discrete plot of a sum.  The OP, in that case, was looking to plot a partial sum of:

sum(1/(n^3*sin^2(n)),n=1..400);

The solution that was suggested to the OP was the following:
 

with(DynamicSystems):
f[N1_]=sum[1/(n^3*sin^2[n]), {n,1,N1}];
DiscretePlot[f[x], {x,0,400},PlotRange->All];

Using that solution as a guide, I tried this:
 

with(DynamicSystems):
lambda := 15.4;
f[N1_] = sum[exp(-lambda)*lambda^n/n!, {0, N1, n}];
DiscretePlot[f[x], {0, 50, x}, PlotRange -> All];

In the attached image, you can see the output from my attempt.

I am completely flummoxed. Any suggestions would be greatly appreciated.