Nm= p1.
p2 ...pm + 1, for m more than or equals 1.

So N1 = p1 + 1 = 2 + 1 = 3, N2 = p1 p2 + 1 = 2
3 + 1 = 7, etc.

We prove that Nm is not divisible by any of p1, p2, . . . , pm, so that Nm is either a prime or it is divisible by a prime larger than pm.

(c) Use Maple to find out which of these numbers Nm, for m = 1, 2, . . . , 15, is actually prime.

Use Maple to compare pm with the smallest prime number that divides Nm, for m =1, 2, . . . , 15.

So I am brand new to using Maple and am trying to solve a system of ODE's. When I try and use dsolve however I reciever the error "function expected" and I am not sure why. The code I have is as follows

sys_ode := diff(A(t), t) = a*(A(t)+B(t))-b1*A(t)*(C(t)+D(t))-g1*A(t)*B(t), diff(B(t), t) = -b2*B(t)*(C(t)+D(t))+g1*A(t)*B(t), diff(C(t), t) = e*(b1*A(t)+b2*B(t))*C(t)-g2*b2*C(t)*B(t)-g3*C(t)*D(t), diff(D(t), t) = g2*b2*C(t)*D(t)-m*D(t)

dsolve([sys_ode])

I'm sure it is something simple, but I just can't seem to figure out where the issue is. Any help is much appreciated.

I'm trying to substitute one Differential equation into another differential equation.

eq1:=d*n(t)/dt = (rho(t)-beta)*n(t)/Lambda+lambda*C(t)+q

eq2:=diff(eq1, t)

resulting in -> eq2 := d*(diff(n(t), t))/dt = (diff(rho(t), t))*n(t)/Lambda+(rho(t)-beta)*(diff(n(t), t))/Lambda+lambda*(diff(C(t), t))

then I'm given that (diff(C(t), t)) is given by another equation: 

eq3:=d*C(t)/dt = beta*n(t)/Lambda-lambda*C(t)

At this point I'm trying to substitute equation 3 into equation 2 for diff(C(t),t)

eq4 := subs(diff(C(t), t) = rhs(eq2), eq5)

however no matter what way's I try this I get an error: 

Error, (in simpl/reloprod) invalid terms in product: (d*(diff(n(t), t))/dt = (diff(rho(t), t))*n(t)/Lambda+(rho(t)-beta)*(diff(n(t), t))/Lambda+lambda*(diff(C(t), t)))^-1

I then tried to map it but again i got an error specifically about the first parameter: 

Error, invalid operator parameter name

eq5:=map((d/dt C(t))->beta/Lambda*n(t)-lambda*C(t),eq2)

I'm just wondering if what I am trying to do is even possible in Maple? 

If anyone can help I would greatly appreciate it! 

The system :

Hello ..

I need codes bifurcation and cobweb for logistic map by maple ..

can you help me ?

Thanks ...

I am trying to recreate journal work for validating using another computer program so I am trying to use maple to solve the ODE, based on further research I found using laplace might be the best but I am having some trouble.

eq8:=d*(n(t)+C(t))/drho = -rho(t)/(l*alpha*K_c)

given the initial conditions of:

ICs:= n(0) = n_0, rho(0) = rho_0, C(0) = (beta-rho_0)*n_0/(l*lambda)

therefore: 

equation9 := dsolve({equation8, ICs}, {C(t), n(t)}, method = laplace)

Following this process I get the error: 

Error, (in dsolve) invalid initial condition

According to the journal work the solution I am looking for is: 

C(t)=-n(t)+(rho_0^2+rho(t)^2)/(2*l*alpha*K_c)+((Beta+l*lambda-rho_0)*n_0)/(l*lambda)

is there something that I'm doing wrong or missing? 

Any help would be greatly Appreciated! 

I want to simulate Inelastic collision

There has 2 ball which I can change  Quality and Radius.

one ball  move to another stirless ball with diferent angle

I am trying to define some functions by recursion relations but I get this error:

Error, (in type/realcons) too many levels of recursion

Please see the attached file for details:

File.mw

I`ve written this program:



It determines the power of the Ithfactor, only one problem. When I try to compute the factor of 512 I get 2_2 (subscript) which is extremly stupid. I want to know how I can change my program so that it can give the power of every ith factor. Thanks!

Dear all,

I woul line integrate e gradient moltiply to a vector. I will try to explain better.

I have define this function:

and I have a vector a

v=(0,1)

I would like to applu a dot product between the gradient of phi an the vector and integrate the results. 

I have already try in many way but without succeed. Someone could please help me?

Thanks 

Hi all

I'm having a hard time, making Maple plot a pretty huge expression in my project.

I have solved a differential equation with initial conditions with method=laplace. The differential equation contains a fourier serie equation, so the more accurate i want the equation to be, the larger the differential equation will be.

Maple solves the equation just fine, and i can plot the solution with 2-4 fourier parts, but when i go higher as i need, the graph ends up empty?

with 20 parts i get the following equation: 

0.*sin(52.88*t)+0.*cos(74.03*t)-0.*sin(74.03*t)-0.*cos(52.88*t)+0.*cos(200.95*t)-0.*sin(200.95*t)+0.*cos(158.65*t)-5.55*10^(-8)*sin(105.76*t)-0.*sin(116.34*t)+0.*cos(31.73*t)-.45*sin(10.58*t)+1.02*cos(10.58*t)+0.*sin(95.19*t)+0.*cos(116.34*t)+0.*sin(179.80*t)-0.*cos(179.80*t)+0.*sin(137.49*t)-0.*sin(31.73*t)-0.*cos(95.19*t)+5.53*10^(-993)*(-1.13*10^992*cos(10.61*t)+8.14*10^991*sin(10.61*t))*exp(-0.7e-1*t)+4.23*10^(-7)*cos(211.53*t)-6.69*10^(-7)*cos(63.46*t)-6.11*10^(-7)*cos(105.76*t)+5.79*10^(-7)*cos(126.92*t)+6.67*10^(-8)*sin(42.31*t)-5.88*10^(-8)*sin(148.07*t)+5.88*10^(-8)*sin(211.53*t)+7.09*10^(-7)*cos(42.31*t)+5.45*10^(-8)*sin(84.61*t)+6.40*10^(-7)*cos(84.61*t)+5.72*10^(-8)*sin(126.92*t)-9.01*10^(-7)*cos(21.15*t)+5.97*10^(-8)*sin(169.22*t)+5.06*10^(-7)*cos(169.22*t)-5.98*10^(-8)*sin(190.38*t)-4.65*10^(-7)*cos(190.38*t)-5.44*10^(-7)*cos(148.07*t)-1.33*10^(-7)*sin(21.15*t)-5.61*10^(-8)*sin(63.46*t)-0.*cos(137.49*t)-0.*sin(158.65*t)

if i plot that expression, the graph ends up empty?

I did also try to solve the equation numerical to plot it with odeplot, but when i try to solve it without the laplace method i get this error message:
"Error, (in dsolve) found the following equations not depending on the unknows of the input system:"

The differential equation is:

ode:=diff(Theta(t), t, t)+2*Zeta*omega[balanceue]*(diff(Theta(t), t))+omega[balanceue]^2*Theta(t) = M[p]/m[balanceue]

and the initial conditions:

ICS := Theta(0) = (1/8)*Pi, (D(Theta))(0) = 0;

when i do:

dsolve({ICS, ode}, Theta(t), method = laplace) it solves just fine.

but when i try with:

dsolve({ICS, ode}, Theta(t))

or

dsolve({ICS, ode}, Theta(t),numeric)

I get the message: 

Error, (in dsolve) found the following equations not depending on the unknowns of the input system: {Theta(0) = (1/8)*Pi, (D(Theta))(0) = 0}

It doesnt seem logical at all, is it a bug? Or can anybody help me with this problem?

Regards

Nicolai

I have denotation like A[0], A[1], A[2], A[3]... But one package doesn't allow to use indexed variables.

I'd like to change denotation. For example, to A0, A1, A2, A3, but I don't know how to do it automatically...

I am trying to perform the following manipulation (This is a minimum working example).

a < b  < c;
(1)*2;

Error, invalid terms in product: a < b and b < c

 Can anyone tell if it is possible to manipulate inequalities exactly as it is the case with equations?

I have two Reissner Nordstrom black holes that are near extreme. How do I show they move? 

When i copy expression and past it in word, i can change the size of the picture whitout loosing the detials.

How can i export the expression to a file, such that when i will open it in word i could change the size without loosing details? much thnks :)