Consider the planar region bounded by the curves y=x^2 + 4x + 6, x=1 and y=2. Rotate about the line x=3 

i want to display a1,a2,a3 on one side and b1,b2,b3 on another sides  that is in an array (one side a1,a2,a3 graph and another side b1,b2,b3 graphs)

 my codes are 

a1:=plot(x^2,x=1..2);a2:=plot(x^2,x=3..4);a3:=plot(x^2,x=4..5);
b1:=plot(x^3,x=1..2);b2:=plot(x^4,x=1..1);b3:=plot(x^5,x=1..2);

A:=([a1,a2,a3])

B:=([b1,b2,b3])

with(plots):display([A,B]);

Hi 

I want to plot:

plot(int(x^2*e^(PI*x/10^12 M), x = 0 .. infinity))

where M is 10^8 but i want the units of the x axis to be x/M, is this possible if so please can i have some help.

Hi all,

There is a state space system as follows:

X_dot = A.X + B.U

Y = C. X

The transfer function of this system is

H(s) = C. (sI-A)^-1.B                 (1)

H(s), that is calculated in Maple by using the formula (1), is very complex (please see the attached file).

Would you please show me how to find the transfer function of the system? 

In addition, how to plot its bode plot?

I'm trying to plot this equation with initial conditions. I'm following the exact same code that my professor used, but it isn't working. 

restart; with(plots): with(DEtools):

DE := y(t)^2*(4-y(t)^2);

p1 := DEplot( DE, y(t), t=0..1/2, {[y(0)=1],[y(0)=3]}, arrows=none, thickness=2, linecolor=[blue,green]):

Under the last line of code I get the error "Error, (in DEtools/DEplot/CheckInitial) the 'number' option...

Possibly a quick one,

I have two sets of precedures, both return the right values however one returns the values alone e.g.

x := 123123.2312323

The other wrongly returns

x :=  Return(123123.2312323)

I have literally pages of code, any suggestions why this would be the case?

Thank you!!!

Can we find exact solution of PDE with help of maple?

If it is yes the please let me how to write commands for parabolic pde 

u_t=k u_xx with boundary conditions u(0,t)=u(l,t)=0 and initial condition u(x,0)=e^-x

I have a question about the "applyrule" function. I have an expression that looks like:

exprsn1:=sin(z)(Acos(z) + Bcos(x)cos(2z) + Ccos(x));

I would like to express cos(2z)sin(z) in the form 1/2(sin(3z)-sin(z)). However, when I use the "applyrule" as:

exprsn2 := applyrule(sin(z)cos(2z)=1/2(sin(3z)-sin(z), exprsn1);

it returns the same expression. Is there something else I should be doing before I use applyrule?

I have a model with 3 parameter blocks containing a lot of constants. I would like to know how can i print those lists to papper, .pdf, word or something like that.

Thanks

Hi Guys, I hope you can help!

I am stuck trying to plot a function that has a conditional based on 'x' the plotted variable.

e.g.

> pltx := proc (num1::integer, num2::integer, x) 
if 1 < x then 
return 0 
elif x < 0 then 
return 0 
else 
return num1*x^num2 
end if 
end proc;

> d := plot(pltx(2, 5, x), x = 0.1 .. 0.4);

Error, (in pltx) cannot determine if this expression is true or false: 1 < x

Regards,

I am trying to solve a problem in theoretical applied mechanics in which I need to find the solution for 5 unknown but I have 6 equations. None of the equations is redundant and I need the solution of the variable to satisfy each of the equations.

These are the equations 

24-52*L^4*m[0]*Omega^2*(R+1/2)/EI[0]-168*beta[1]-52*F*L^2/EI[0]+104*beta[2]+(120+72*beta[1]+104*L^4*m[0]*Omega^2*R/EI[0]-504*beta[2]+312*beta[3]+84*L^4*m[0]*Omega^2*(R+1/2)/EI[0...

Hello,

I have created a test document to check an equation and I would like to feed a range of k values (ranging from 0 to say 512) in increments of 1) to the main function, and then plot the result as a graph of Zuncomp (y-axis) vs k(x-axis).

I can verify the values of the function by typing say Zuncomp(0) or Zuncomp(10) etc but how do you run this so that you do Zuncomp(k) to read all values of k and generate the output to plot?

Using Maple 13 (Windows 7...

I am new to maplesim. I want to model the human hand and constrain the fingers to a certain range of angle. Then apply a controller to analyse its controllability and also see the dynamic equation formulation

restart;

assume(alph>0);

assume(alph,real);

f_exp:=exp(-alph*r^2);

ff_deriv:=simplify(r^2*f_exp^2);

ff:=simplify(int(ff_deriv,r=0..infinity));

It seems to understand what i'm trying to integrate but when I try and find ff it comes up with an awful expression full of erf functions and with a lim r--> infinity at the front. However on wolfram alpha I get the answer that I want: sqrt(pi/2)/(8*alpha^(3/2)).

I am trying to solve two differential equations numericaly: