Items tagged with academic academic Tagged Items Feed

I typed (x-a)(x-a) into Maple 18 (the only thing in a new document) and got x(x-a)-a 

I'm wondering why it is evaulated like this?

Is this because of a setting?

BTW, I evaluated at a=4 first. 

How do you use algcurves and Puiseux series in Maple?

Hi,

I am trying to plot only the real part of the values obtained from the MapleSim model.  Probably there is a Modelica function to get only the real part of the values, but I am unable to find it. Could you please have a look at the code and let me know what I should do to solve this problem?

Best

Onder

 

The code is: (I want to plot the real part of R2)

 

model IKTest

parameter Boolean switched = false;
parameter Modelica.SIunits.Length L1 = 0.100;
parameter Modelica.SIunits.Length L2 = 0.102;
parameter Modelica.SIunits.Length L3 = 0.5;
//parameter Real c; //Complex number
//parameter c;

Modelica.Blocks.Interfaces.RealInput u[3]
annotation (Placement(transformation(extent={{-140,-20},{-100,20}},rotation=0)));
Modelica.Blocks.Interfaces.RealOutput y[3]
annotation (Placement(transformation(extent={{90,-10},{110,10}},rotation=0)));

Real Tx = u[1];
Real Ty = u[2];
Real Tz = u[3];
Real R1 = y[1];
Real R2 = y[2];
Real R3 = y[3];

equation

R1=atan(Ty/Tz);
R2=(acos((L1^2+L2^2-(sqrt(Tx^2+Ty^2+Tz^2))^2)/(2*L1*L2)));

//2*Modelica.Math.asin(1.0)
//R3=R2.re;

annotation (uses(Modelica(version = "3.1")));
end IKTest;

I'm taking my first steps with maple and pdsolve, trying to run the example in the maplesoft support page:

http://www.maplesoft.com/support/help/Maple/view.aspx?path=examples/pdsolve_boundaryconditions

which reads

>
> restart; with(PDEtools);
> U := diff_table(u(x, t));
>

and I get a solution that is different from the web page, and when i run

Im using maple 13. Any tips about what's wrong?

 

regards

Hi there

I'm working with maple 13 (2d input) and I would like to hide the input in a worksheet, as follows:

view > show/hide contents > (unchecking) input

however i want to insert comments as text that will remain in my worksheet. Using '#' doesn't help, coz they get hidden too.

Couldn't sort out how to do this, any tips?

cheers

Hi there

I'm taking my first steps with maple and my question is not exactly briliiant

Putting in a simple way, I'm trying to solve sth that looks like this

> li := a*[cos(7*x+8*x*y+9*y)*cos(-7*x-8*x*y-9*y)];
             a [cos(7 x + 8 x y + 9 y) cos(-7 x - 8 x y - 9 y)]

and i want to replace

with

getting (if i do it by hand)

> algsubs(7*x+8*x*y+9*y = 20*y, li);
                    a [cos(20 y) cos(-7 x - 8 x y - 9 y)]

my question, is there an elegant way to do it? sth that would look like this?

 

 

thank you for any tips

Hi There

I'm getting started with maple and facing a doubt about using subs, somehow it does not seem to work, this is my expression:

 

> restart;
> theta := omega*t-k*x;
                                omega t - k x
> phi[1] := -(1/2)*H*c*cosh(k*(z+h))*sin(theta)/sinh(k*h);
                   H c cosh(k (z + h)) sin(-omega t + k x)
                   ---------------------------------------
                                 2 sinh(k h)              
> diff(phi[1], t, t);
                                                             2
                H c cosh(k (z + h)) sin(-omega t + k x) omega
              - ----------------------------------------------
                                 2 sinh(k h)                  
> td := simplify(subs(z = 0, %));
                                                          2
                   H c cosh(k h) sin(-omega t + k x) omega
                 - ----------------------------------------
                                 2 sinh(k h)               
> simplify(subs(cosh(k*h)/sinh(k*h) = 1/tanh(k*h), %));
                                                          2
                   H c cosh(k h) sin(-omega t + k x) omega
                 - ----------------------------------------
                                 2 sinh(k h)               
I cannot subs cosh/sinh = tanh. I would like to know why.

Any tips? Thank you

The Interactive Embedded Components in Physics are of great importance today and will be even more in the future. Hereand leave a small tutorial of Embedded Components in Physics applied to physics. I hope that somehow you motives to continue the development of science.

 

  Interactive_Embedded_Components_in_Physics.mw      (in spanish)                 

 

Atte.

Lenin Araujo Castillo

Physics Pure

Computer Science

 

Hi there

I am stuck with the expression

 

> e := (3+2*sinh(x)^2)/(sinh(x)^2*tanh(x));
                                            2
                               3 + 2 sinh(x)  
                              ----------------
                                     2        
                              sinh(x)  tanh(x)

its simplified way, which i cannot find, is


> f := 3*coth(x)^3-coth(x);
                                     3          
                            3 coth(x)  - coth(x)
> simplify(e-f);
                                      0

I would appreciate immensely if someone could give me some pointers on how to solve it.

 

best regards

 

nvc

 

hi maple-friends !

who can help me??

 

QUESTION 1:

i want to extract (export) the data of let me say the function y=x^2 ...

the result shpuld be:

without the variables !  (x=)   (y=)

1         1

2         4

3         9

and so on...

that means just a matrix...of numbers.

 

QUESTION 2:

and waht i have to do if i need more values, including non-integers as x (example: x=1.1, 1.2,1.2 ... )

sorry, im a maple beginner and a lot of the advices in the books are not very helpful !!!   maybe only written

for experts ??   so by the way, does anybody know about a REALLL good BEGINNER-book for maple?

thanx a lot !   :-)

 

Let x:=2*[x_1,x_2,...,x_n];   y:=3*[x_1,x_2,...,x_n]; If I input x+y, it will turn out  5*[2x_1,2x_2,...,2x_n]; But I want to get  5*[x_1,x_2,...,x_n]. What should I do? Thanks!!

Hi all

I have a mathematical problem and I asked it in various sites but the answers till yet are not correct.

Assume that we have:

T[m]:=t->t^m:
b[n,m]:=unapply(piecewise(t>=(n-1)*tj/N and t<n*tj/N, T[m](N*t-(n-1)*tj), 0), t):

where n,N,tj are known constants. furthermore assume that we want to comute the following integral:

for following approximations:

I have written the following code but it seems to be incorrect:

V1:=Vector([seq(seq(b[n,m](t),m=0..1),n=1..3)]);
V:=evalf(V1.Transpose(V1));

the original program is :

taaylor.mws

I will be so grateful if any one can help me to solve it by maple

Mahmood   Dadkhah

Ph.D Candidate

Applied Mathematics Department

Basically I have to find omega by solving determinant of the following matrix:

M := Matrix(4, 4, {(1, 1) = 0, (1, 2) = -1, (1, 3) = 0, (1, 4) = 1, (2, 1) = -EI*beta^3, (2, 2) = -m*omega^2, (2, 3) = EI*beta^3, (2, 4) = -m*omega^2, (3, 1) = EI*sin(147*beta)*beta+k_r*cos(147*beta)+I*c_r*omega*cos(147*beta), (3, 2) = EI*cos(147*beta)*beta-k_r*sin(147*beta)-I*c_r*omega*sin(147*beta), (3, 3) = -EI*sinh(147*beta)*beta+k_r*cosh(147*beta)+I*c_r*omega*cosh(147*beta), (3, 4) = -EI*cosh(147*beta)*beta+k_r*sinh(147*beta)+I*c_r*omega*sinh(147*beta), (4, 1) = -EI*cos(147*beta)*beta^3+sin(147*beta)*k_h, (4, 2) = EI*sin(147*beta)*beta^3+cos(147*beta)*k_h, (4, 3) = EI*cosh(147*beta)*beta^3+sinh(147*beta)*k_h, (4, 4) = EI*sinh(147*beta)*beta^3+cosh(147*beta)*k_h}):

The remaining values are:

beta=((5000/10^12)*(omega^2))^(1/4),k_r=3.33*10^10,k_h=1.62*10^9,c_r=3.14*10^9,m=350000,L=147,EI=10^12:

What is the proper way to deal with this problem numerically. Or maybe it is even possible to get a reasonable analytical expression?

hello, I put in the code:

> with(plots), with(ColorTools), with(LinearAlgebra), with(RandomTools), with(ExcelTools);
> A := `<|>`(`<,>`(1, 2, 0, 2, 3, 4, 3, 4, 7, 9, 5, 3, 4, 6, 7, 8, 3), `<,>`(0, 4, 7, 2, 2, 2, 4, 5, 6.5, 7, 5, 3, 2, 5, 9, 0, 1), `<,>`(1, 5, 2, 0, 4, 1, 2, 3, 4.3, 7, 8, 5, 3, 2, 9, 6, 4)); J := convert(Import("testB1.xlsx", "Cartesian", "E2:G18"), matrix);
Error, (in convert/matrix) expecting array, rtable or list
> B := matrix([[1], [.2], [.1], [.8], [.5], [.6], [.8], [.764], [.234], [0.4e-1], [.89], [.36], [.687], [.627], [.689], [.328], [.139]]); H := convert(Import("test.xlsx", "Cartesian", "D2:D18"), matrix);
Error, (in convert/matrix) expecting array, rtable or list
> C := [seq(Color([H[i, 1], 0, 0]), i = 1 .. RowDimension(A))];

 

And i get the error message everytime I try import m data list and I dont understand why. Any help would be appreciated. Thank you

 

Jennifer

Can someone please help me to find first 5 coplex valued roots of the followind expression:

Frequency_Equation := -(2.220315292*10^21*I)*csgn(omega)*omega^4*sin(1.236117731*(omega^2)^(1/4))*cosh(1.236117731*(omega^2)^(1/4))-2.354665581*10^22*csgn(omega)*omega^3*sin(1.236117731*(omega^2)^(1/4))*cosh(1.236117731*(omega^2)^(1/4))-2.354665581*10^22*sinh(1.236117731*(omega^2)^(1/4))*csgn(omega)*omega^3*cos(1.236117731*(omega^2)^(1/4))-6.415296703*10^25*(omega^2)^(3/4)*cos(1.236117731*(omega^2)^(1/4))*cosh(1.236117731*(omega^2)^(1/4))-2.177438228*10^22*cosh(1.236117731*(omega^2)^(1/4))*(omega^2)^(3/4)*omega^2*cos(1.236117731*(omega^2)^(1/4))-5.946035576*10^21*omega^2*(omega^2)^(3/4)-6.415296702*10^25*(omega^2)^(3/4)-(3.560760000*10^24*I)*sinh(1.236117731*(omega^2)^(1/4))*omega^3*cos(1.236117731*(omega^2)^(1/4))+(3.560760000*10^24*I)*cosh(1.236117731*(omega^2)^(1/4))*omega^3*sin(1.236117731*(omega^2)^(1/4))-(2.613877239*10^21*I)*cosh(1.236117731*(omega^2)^(1/4))*(omega^2)^(3/4)*omega^3*cos(1.236117731*(omega^2)^(1/4))-(6.049258752*10^24*I)*(omega^2)^(3/4)*omega-1.907153069*10^25*sinh(1.236117731*(omega^2)^(1/4))*omega^2*sin(1.236117731*(omega^2)^(1/4))*(omega^2)^(1/4)-(2.220315292*10^21*I)*sinh(1.236117731*(omega^2)^(1/4))*csgn(omega)*omega^4*cos(1.236117731*(omega^2)^(1/4))-(6.049258753*10^24*I)*(omega^2)^(3/4)*cos(1.236117731*(omega^2)^(1/4))*cosh(1.236117731*(omega^2)^(1/4))*omega-3.500000000*10^21*cosh(1.236117731*(omega^2)^(1/4))*omega^4*sin(1.236117731*(omega^2)^(1/4))+3.500000000*10^21*sinh(1.236117731*(omega^2)^(1/4))*omega^4*cos(1.236117731*(omega^2)^(1/4))-2.156220000*10^25*sinh(1.236117731*(omega^2)^(1/4))*omega^2*cos(1.236117731*(omega^2)^(1/4))+2.156220000*10^25*cosh(1.236117731*(omega^2)^(1/4))*omega^2*sin(1.236117731*(omega^2)^(1/4));

1 2 3 4 5 6 7 Last Page 1 of 23