How can i solve for G(p,q,r,t) in the following pde

diff(G(p, q, r, t), t) = {lambda*(p-1)+gamma*(r-1)}*G(p, q, r, t)+delta*(1-p)*(diff(G(p, q, r, t), p))+kappa*(N*r-q)*(diff(G(p, q, r, t), q))+mu*(1-r)*(diff(G(p, q, r, t), r))+beta*(q*e^(-`ρτ`)-p*r)*(diff(G(p, q, r, t), r, p))

Dear All

I trying to solve the following in MAPLE, the code works fine in MATLAB but I am not sure why it doesn't give me values in MAPLE

n:=3: alpha_p:=1: p:=2: mu0:=4*Pi*1e-7: Br:=1.12:

A1:=unapply(sin((n*p+1)*alpha_p*Pi/(2*p))/((n*p+1)*alpha_p*Pi/(2*p)),n):

A2:=unapply(sin((n*p-1)*alpha_p*Pi/(2*p))/((n*p-1)*alpha_p*Pi/(2*p)),n):

M1:=unapply((Br/mu0)*alpha_p*(A1(n)+A2(n)),n):

M2:=unapply((Br/mu0)*alpha_p*(A1(n)-A2(n)),n):

How can i solve for X(t) in the following equation

diff(mean(X(t)), t) = lambda-delta*mean(X(t))-beta*mean(X(t)*V(t))

and how can i plot X(t) against t.

Here is what I tried.

I need to get the sum of all the real roots of 3*sin(x^2) = x

So I tried this. 

>solve(3*sin(x^2)-x = 0, x, AllSolutions)

Which gaved me this 3*sin(RootOf(_Z-9*sin(_Z)^2))

Then >sum(3*sin(RootOf(_Z-9*sin(_Z)^2)), _Z = -infinity .. infinity)

Which gives me this answer which I cannot make sense of.  signum(sin(RootOf(_Z-9*sin(_Z)^2)))*infinity

I'm trying to solve for n in a simple sum() function:  Here is the function.

I have a question regarding on how to graph and prove that csc ^2(theta)-cot ^2(theta)=1 our professor provided us with a very poor example, i have to graph the left side of the equation first with a 0 to 2*Pi x-axis and then graph the right side on the same plot. 

I examined some rules to verify that if a certain function for example f(x)=x^2+5 is a 1-1 function or not.

What I did after defining the function was as follows:

[> assume(x>0,s>0);

[> solve(f(x)=f(s),{x,s});

OR

[> solve(f(x)=f(s),{x,s}) assuming x:: positive;

They didn't give me just

s=x

at all. Thanks so much.

Could anyone help with solving this integral?

> evalf((∫)[-infinity]^(infinity)((∫)[-infinity]^(v)( uv)/(sqrt(2 Pi))(1/2+1/2 erf((1 )/2 sqrt(2) u))(1/2-1/2 erf((1 )/2 sqrt(2) v))(e)^(-((u^(2+v^(2))))/(2)) ⅆu )dv);

thanks,

Hi!

I have a piecewise defined function alpha. Its functional directives are stored in a vector A. But when I call up the command

alpha:=(x1,x2,x3)->piecewise(
    x1 < 0 , A[1],

   0< x1<1 A[2], ...):,

the result (x1,x2,x3)->piecewise(.... , A[1], ...) remains symbolic. What can I do that the entries of A[i] appear in the functional directive?

Thanks!

I'm new to using Maple. I'm trying to evaluate an equation where all unknowns are given. So to say use Maple

as a kind of desktop calculator with the possibility to change certain input variables and see what the outcome is.

I'm wondering why the trigonometric functions are not evaluated.

> restart; with(plots); with(LinearAlgebra); with(linalg);
> m1 := 2280; m2 := 2243; m3 := 2220; k1 := 300; k2 := 200; k3 := 100; F0 := 4.1; omega1 := sqrt((k1+k2)/m1);
> Sys := m1*(diff(u1(t), t, t))+(k1+k2)*u1(t)-k2*u2(t) = F0*cos(omega*t);
print(??); # input placeholder
> IC := u1(0) = 0, (D[1](u1))(0) = 0;
> csol := dsolve({IC, Sys}, {u1(t)});
%;
> p1 := subs(csol, u1(t));
> 
> plot(p1, t = 0 .. 10*Pi, 'color' = "red", 'thickness' = 2, labels = ...

I have 5 term which each term includes three function of time (autonomous system ), i want to obtain time derivative of each term as other terms for example if i define each term as a dimensionless variable , e.g.  x , y , z , v , w , i want to know dx/dt = as  x , y , z , v , w.

i differentiate each term but i dont know how to obtain each differentiated term as other terms

Hello,

I need to plot M vs phi from the following ode. For this I am using implicitplot but unfortunately

I am getting empty plot. Please guide me through this issue.

Thanks 

eq1 := 1/(1-phi)^2.5/(1+4117/783*phi)*diff(f(eta),`$`(eta,3))+1/2*f(eta)*diff(f(eta),`$`(eta,2))

+M*(1-diff(f(eta),eta)) = 0;

bc:=f(0)=0,D(f)(0)=beta/(1-phi)^(2.5)*(D@D)(f)(0),D(f)(N)=1;

respar:=dsolve({eq1,f(0)=0,D(f)(0)=beta/(1-phi)^(2.5)*(D@D)(f)(0)},numeric,

I have a program in  matlab,

I want to change below code in maple,

henon_map = inline('[a-x.^2+b*y x]','a','b','x','y');

%a=input('a [0] = ');
b=input('b [0.4] = ');

%if(isempty(a))
% a=0
%end

if(isempty(b))
 b=0.4
end

Npts=200;
Ntrans=100;
Ntot=Ntrans+Npts;

% For the whole window:
fig1=1; fig2=2;xmin=-2;xmax=2; Nmin=1; Nmax=50;
amin = -0.09;
amax = 1.249;

I am getting the reecursive assignment error on the folowing code in the band[i]:= [band[]i[], ... section.

CreateMatrix is defined and I know that works and creates a global Matrix H