Dear all

I tried to compute an integral using trapezoidal rule using two differents methods but I get two different results
trap.mw

Thank you for your help 

Dear all

I would like to sketch the solution of first ode, but by hand we can not get the solution. I used maple, to  solve the ode but unfortunattely there is no solution obtained so I can sketech the behavior of the solution. 

Is there any way to sketch the solution behavior .

sketch_solution.mw

thank you 

In using maple 17, for a for loop, the alias command is valid outside the for loop, but entering the for loop, maple verifies that alias is not recognized.The code is as follows:
with(PDEtools):
for k from 1 to 4 do
alias(u[k]=u[k](x,y,t));
alias(a[k]=a[k](x,y,t));
od:
print(has(a[1],x));
sol:=[]:
for i from 1 to 4 do
sol:=[sol[],a[i]=rhs(isolate(a[i]+diff(a[i],x)+i=0,diff(a[i],x)))];
print(has(a[i],x));
end do;
The output is as follows: thequestionofalias
true
[a[1] = a[1] + 1]
false
[a[1] = a[1] + 1, a[2] = a[2] + 2]
false
[a[1] = a[1] + 1, a[2] = a[2] + 2, a[3] = a[3] + 3]
false
[a[1] = a[1] + 1, a[2] = a[2] + 2, a[3] = a[3] + 3, a[4] = a[4] + 4]
false
Why is it outside the loop? A [1] ... A [4] has the variable x, but either a [ I ] or a [1] has no variable X.Send the question to the Forum and see which expert can answer it. See how you can make the alias command work inside a for loop! I'd appreciate it if you could answer that!

How to solve this integral equation?

rho:=1/(2);

mu:=1/(2);

T:=1/(4)*(t-x)^(1/(2));

E[rho,mu](T):=Sum((T^(k))/(GAMMA(k*rho+ mu)),k=0..5);
                       
eq1 := h(t) = (8/3)*(int((t-x)^(1/2)*E[rho, mu](T)*(1/(sqrt(Pi)*x^(1/2))-(3/32)*x^2+(3/16)*exp(-x)-1), t = -1 .. 1));       

 

Download 1.mw

Hi!

There is some procedure/function to convert C code into maple code? For instance, the above text file

C_code.txt

contains the C code of this paper  http://www.dcs.bbk.ac.uk/~jkl/pubs/JL1_00a.pdf to compute the image of a point in [0,1] under the N-dimensional Hilbert curve and vice versa (i.e., given a point P in the N-dimensional cube [0,1]^{N} find the point t in [0,1] such that t is mapped into P under the Hilbert curve)

Unfortunately, I don't know anything about the C programming language.

Many thanks in advance for your comments.

When i trying to solve the BVP the following error will found 

Error, (in fproc) unable to store 'HFloat(1.0141653815612295)+HFloat(0.0020355870386861838)*I' when datatype=float[8]
 

I have attached my worksheet.

Please anyone help me to get solution to this problem.

Thank you so much

Download Main.mw

I'm trying to solve the problem which is published in the research article. The paper is attached. 

This work done with the help of shooting technique in maple 18

1.Is it possible to solve in maple 17

2. How to obtain the triple soln

3. How to implement Shooting technique

Please anyone help me to get this 

Download Triple.mw

Hi

I want to write this system and solve it by picard iteration method

but I don't know how to write a system of ODE in maple

x'=cos (x)    x(0)=1

y'= sin (z)    y(0)=-1

z'= zy           z(0)=2

 

In the below paper 

ODEs are solved by using Matlab Software

Is it possible to solve this in Matlab software?  [I think you mean Maple!--Carl Love]

I have attached the Article

j.physa.2019.123959.pdf

Hi 
I need to represent the contourplot of sum a function as

contourplot(sum(BesselJ(n,r),n=0..5),r=-10..10).
can anyone know what the problem is ?

Thanks

I m Stucked with this step

Download Pblm2.mw

When I try to solve this problem above error displays

Here is the file
 

Download Main.mw

When i try to solve the problem i got error like this 

Please help me get the solution

Download 1.mw
 

Hi

cant get the answer, any comments plz

f := (x, y) ->(1/2)*a*(sinh(y-x^2)+tanh(x-y^3));

evalf(int(int(f(x, y), x = -6 .. 5), y = -5 .. 5)):

One forum had a topic related to such a platform. ( Edited: the video is no longer available.)

The manufacturer calls the three-degrees platform, that is, having three degrees of freedom. Three cranks rotate, and the platform is connected to them by connecting rods through ball joints. The movable beam (rocker arm) has torsion springs.  I counted 4 degrees of freedom, because when all three cranks are locked, the platform remains mobile, which is camouflaged by the springs of the rocker arm. Actually, the topic on the forum arose due to problems with the work of this platform. Neither the designers nor those who operate the platform take into account this additional fourth, so-called parasitic degree of freedom. Obviously, if we will to move the rocker with the locked  cranks , the platform will move.
Based on this parasitic movement and a similar platform design, a very simple device is proposed that has one degree of freedom and is, in fact, a spatial linkage mechanism. We remove 3 cranks, keep the connecting rods, convert the rocker arm into a crank and get such movements that will not be worse (will not yield) to the movements of the platform with 6 degrees of freedom. And by changing the length of the crank, the plane of its rotation, etc., we can create simple structures with the required design trajectories of movement and one degree of freedom.
Two examples (two pictures for each example). The crank rotates in the vertical plane (side view and top view)
PLAT_1.mw


and the crank rotates in the horizontal plane (side view and top view).

The program consists of three parts. 1 choice of starting position, 2 calculation of the trajectory, 3 design of the picture.  Similar to the programm  in this topic.