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...

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...

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:

I recently get a new macbook pro with retina display. Unfortunately Maple 16 seems to be not retina ready. Searching the net I found a tool called retinizer. Using it Maple 16 seems to be displaying a little better, but still the outputs of calculations looks blured. however, whenever certain texts are selected, the whole Maple window display becomes very sharp, so I guess Maple is somehow retina ready but this functionality is disabled by default. Checking the "Get Info" window...

Hi 

In trying to solve:

 a := diff(u(t), t$2)+(p^2-I+t^2)*u(t) = 0;

I get the following solutions,              

sol1 := u(t) = _C1*WhittakerM(-1/4-(1/4*I)*p^2, 1/4, I*t^2)/sqrt(t)+_C2*WhittakerW(-1/4-(1/4*I)*p^2, 1/4, I*t^2)/sqrt(t)        

It should be possible to expand these into parabolic cylinder functions and but im not sure how, i would appreciate any help.

thanks

Hi I have three differential equations: 

u := diff(P(t), t) = -7*10^(-8)*P(t)*t/(P(t)*t+R(t))^(1/2),

diff(R(t), t) = 7*10^(-8)*t^2*P(t)/(P(t)*t+R(t))^(1/2)+600*(Z(t)^2-10^5*t^3*(1/(1.15*10^12))^(2/3)*e^(-1.15*10^12))/(t*(P(t)*t+R(t))^(1/2)),

diff(Z(t), t) = -4*10^5*(Z(t)^2-10^5*t^3*(1/(1.15*10^12))^(2/3)*e^(-1.15*10^12))/(t^2*(P(t)*t+R(t))^(1/2))

and i want to solve them with initial conditions:

initial := R(0) = 0, Z(0) = 0, P(0) = P;

Hello and thank you in advance,

I am used to Mathematica and am looking to switch, conditionals in Mathematica seem easy, like say

How do you replicate this in an easy to read way without lots of nested if statements?

Hi,

r:=[3,3,3,3,4,4,4,3,3,3,3,3,3,3,3,2,4,3,1];

What's the shortest way to create a frequency table from it?

Ie, tells me elements are 1,2,3,4

the corresponding counts are: 1,1,13 and 4

and then if it can return the number with highest counts: 3

Many thanks!

Casper

Hellow, can any body give the suggestion regarding multiple plot in array

i have 

V:=(2*x+3*gamma)*y^2

d0:=subs(x=0.1,gamma1=0,V):
d1:=subs(x=0.1,gamma1=2,V):
d2:=subs(x=0.1,gamma1=4,V):

plot([d0,d1, d2], y=0..0.1,  legend = ["gamma = 0", "gamma= 2","gamma= 4"], titlefont = ["ROMAN", 15], labels = ["y ", "v"], labeldirections = ["horizontal", "vertical"], labelfont = ["HELVETICA", 10], linestyle = [solid, longdash,dashdot],color = [black, red,green], axesfont = ["HELVETICA", "ROMAN", 8], legendstyle = [font = ["HELVETICA", 9], location = bottom],axes=boxed);

and 

e0:=subs(x=0.5,gamma1=0,V):
e1:=subs(x=0.5,gamma1=2,V):
e2:=subs(x=0.5,gamma1=4,V):

plot([e0,e1, e2], y=0..0.1,  legend = ["gamma = 0", "gamma= 2","gamma= 4"], titlefont = ["ROMAN", 15], labels = ["y ", "v"], labeldirections = ["horizontal", "vertical"], labelfont = ["HELVETICA", 10], linestyle = [solid, longdash,dashdot],color = [black, red,green], axesfont = ["HELVETICA", "ROMAN", 8], legendstyle = [font = ["HELVETICA", 9], location = bottom],axes=boxed);

I want to display  this side by side  i,e set {do,d1,d2} and {eo,e1,e2} 

 thanks in advance

 exa.doc

Hi, I am trying to plot a solution curve on a vector field with an initial condition of y(0)=0 and I keep getting error messages. This is what I have so far:

with(DEtools); dfieldplot(diff(y(x), x) = x^2+y(x)^2-1, y(x), x = 0 .. 5, y = -1 .. 5, arrows = line, title = 'Slope*Field');

Thank you. 

Greetings,

I have a finite Fourier series of the form, Asin(x)cos(z) + Bsin(2x)cos(z) + Csin(3x)cos(2z) +... I would like to be able to extract the coefficients, A,B,C etc., that correspond to a specific mode. I have tried using the coeff and coeffs functions. They work for a one-dimensional Fourier series (i.e. if S= Asin(x) + Bsin(2x) + Csin(3x); is my series, then coeff(S,sin(2x) returns B). I cannot however get this to work for the 2-D case. Any suggestions?