MaplePrimes Questions

Search Questions:

Latest Questions Latest Questions Feed

my int is in below.

my_int.mw

my int has hyperbolic and trigonal parts thus when they are multiplied , Maple is not able to integrate!

 

Now, How can I solve this integral?

thank you for help

my int is:

 

``

restart

eq1 := m*(diff(diff(w(x, t), t), t))+diff(diff(EIy*(diff(diff(w(x, t), x), x))+(EIz-EIy)*(diff(diff(w(x, t), x), x))*(theta(x)+phi(x, t))^2, x), x)-Pz = 0

l := 16; m := .75; EIy := 0.2e5; EIz := 0.400e7; GJ := 0.1e5; mj := .1; Pz := 5000; Mz := 0; theta := proc (x) options operator, arrow; 0 end proc

w := proc (x, t) options operator, arrow; q[1](t)*(cosh(1.8751*x/l)-cos(1.8751*x/l)+(-1)*.734096*(sinh(1.8751*x/l)-sin(1.8751*x/l)))+q[2](t)*(cosh(4.6941*x/l)-cos(4.6941*x/l)+(-1)*1.01847*(sinh(4.6941*x/l)-sin(4.6941*x/l)))+q[3](t)*(cosh(7.8548*x/l)-cos(7.8548*x/l)+(-1)*.999224*(sinh(7.8548*x/l)-sin(7.8548*x/l))) end proc

phi := proc (x, t) options operator, arrow; q[4](t)*sqrt(2)*sin(1.5708*x/l)+q[5](t)*sqrt(2)*sin(4.7124*x/l)+q[6](t)*sqrt(2)*sin(7.8540*x/l) end proc

NULL

f[1] := int(lhs(eq1)*(cosh(1.8751*x/l)-cos(1.8751*x/l)-.734096*(sinh(1.8751*x/l)-sin(1.8751*x/l))), x = 0 .. l)

``

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

restart;
eq:=diff(y(x),x) - y(x)^2 - x*y(x) - x + 1;
r:=dsolve(eq,y(x));
latex(r);

Gives an error "(in fprintf) string expected for string format". The problem seems when the math contains "erf"? Both examples below fail with result that contains this.

Maple 18, on windows 7.

Another example:

restart;
eq:=diff(y(x),x) + a*y(x)*(y(x)-x) - 1;
r:=dsolve(eq,y(x));
latex(r);

Another example, which seems to fail for different reason

restart;
eq:=diff(y(x),x) + a*x*y(x)^3+b*y(x)^2;
r:=dsolve(eq,y(x));
latex(r);

I've been using a maple package (written by another author) on earlier versions of Maple.  Recently I have been trying to rewrite it so it's a bit more Maple 16+ friendly (since I've read that Maple is not forward-compatible). However, I have been experiencing similar strange behavior in both the original and my current version: 

After saving the module to my personal library, I load with(MyPackage). Most functions work fine (all but a handful out of dozens).  However, some functions won't be recognized (e.g., Coproduct(...) will return itself Coproduct(...)).  

HOWEVER, when I execute the code in the module in its own Maple worksheet, all functions are recognized and work perfectly.  I look at the difference in the code between the functions that work upon loading the package (using with(MyPackage) command in a separate sheet) and those that don't, and there is no difference.

I'm assuming I'll just have to cut and paste the hundreds of lines of code for my package into each sheet I'd like to use it in, since I haven't been able to find others experience this problem in the knowledge base.  I'm hoping someone will be able to help out, however.

hi all.

I have wrore the following program for optimization with bernstein and block pulse hybrid functions.

the program have some errors which i can't understand.

Bernestien1.mws

restart:

alias(C=binomial):
with(LinearAlgebra):
macro(LA= LinearAlgebra):


HybrFunc:=proc(N, M,  tj)               # N=Number of subintervals,  M=Number of functions in subintervals
 
local B, n, m;

global b;

for n from 1 to N do
for m from 0 to M-1 do

B := (i,m,t) -> C(m,i)*(1-t)^(m-i)*t^i:

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


Array(1..N, 0..M-1, (n,m)->b[n,m](t)):

#convert(%,vector);
end proc:

HybrFunc(3, 3, 1);




                                       # End Of Definition
 
g2(t):=t;            #*exp(t-1):                      # Any other function can be replaced here
    

g1(t):=add(add(c[n,m]*b[n,m](t), m=0..2), n=1..3);
Optimization[Minimize](sqrt(int((g2(t)-g1(t))^2, t=0.. 1)));
assign(op(%[2]));
plot([g2(t),g1(t)], t=0..1, 0..5, color=[blue,red],thickness=[1,3],discont, scaling=constrained);

Array(1 .. 3, 0 .. 2, {(1, 0) = piecewise(0 <= t and t < 1/3, (1-3*t)^2, 0), (1, 1) = piecewise(0 <= t and t < 1/3, (6*(1-3*t))*t, 0), (1, 2) = piecewise(0 <= t and t < 1/3, 9*t^2, 0), (2, 0) = piecewise(1/3 <= t and t < 2/3, (2-3*t)^2, 0), (2, 1) = piecewise(1/3 <= t and t < 2/3, (2*(2-3*t))*(3*t-1), 0), (2, 2) = piecewise(1/3 <= t and t < 2/3, (3*t-1)^2, 0), (3, 0) = piecewise(2/3 <= t and t < 1, (3-3*t)^2, 0), (3, 1) = piecewise(2/3 <= t and t < 1, (2*(3-3*t))*(3*t-2), 0), (3, 2) = piecewise(2/3 <= t and t < 1, (3*t-2)^2, 0)}, datatype = anything, storage = rectangular, order = Fortran_order)

g2(t) := t

"g1(t):=c[1,0] ({[[(1-3 t)^2,0<=t and t<1/3],[0,otherwise]])+c[1,1] ({[[6 (1-3 t) t,0<=t and t<1/3],[0,otherwise]])+c[1,2] ({[[9 t^2,0<=t and t<1/3],[0,otherwise]])+c[2,0] ({[[(2-3 t)^2,1/3<=t and t<2/3],[0,otherwise]])+c[2,1] ({[[2 (2-3 t) (3 t-1),1/3<=t and t<2/3],[0,otherwise]])+c[2,2] ({[[(3 t-1)^2,1/3<=t and t<2/3],[0,otherwise]])+c[3,0] ({[[(3-3 t)^2,2/3<=t and t<1],[0,otherwise]])+c[3,1] ({[[2 (3-3 t) (3 t-2),2/3<=t and t<1],[0,otherwise]])+c[3,2] ({[[(3 t-2)^2,2/3<=t and t<1],[0,otherwise]])"

Error, (in Optimization:-NLPSolve) complex value encountered

Error, invalid left hand side in assignment

(1)



Download Bernestien1.mws

 I'll be so grateful if any one can help me.

Mahmood   Dadkhah

Ph.D Candidate

Applied Mathematics Department

can i compute this integral so that the answer does not have any complex part ? i just want the integral returns me real function or number .
 i also use evalf(int ... ) and evalf(Int ... ) and two different answers were given ! what is wrong ?
tnx for help in advance .

restart:

int(log(sec(x)+ tan(x)),x);

-I*ln(exp(I*x))*ln((-I*(exp(I*x))^2+I+2*exp(I*x))/((exp(I*x))^2+1))+(1/2)*Pi*ln(-I*(exp(I*x)-I))+I*dilog(-I*exp(I*x))+I*dilog(-I*(I+exp(I*x)))+I*ln(exp(I*x))*ln(-I*(I+exp(I*x)))

(1)

evalf(Int(log(sec(x)+ tan(x)),x=0..10));

3.493589070+16.61151993*I

(2)

evalf(int(log(sec(x)+ tan(x)),x=0..10));

7.035851984+23.35343547*I

(3)

 

 

Download integral.mw

Solve Integral...

July 12 2014 vhha1972 10
1 4

Hi there,

I am trying to do a Fourier transform of the function

f(t) = 3*cos(2000*Pi*t)

I use Fourier transform definition to solve the following integral in Maple:

Int(3*cos(2000*Pi*t)*exp(-(I*2)*Pi*f*t), t = -infinity .. infinity)

I get the answer as "undefined".

The correct answer is:

(3/2)*delta(f+1000)+(3/2)*delta(f-1000)

Thanks!

Vhha

 

Hi,


Is it possible to solve linear matrix inequality in Maple?
For example, using Matlab and yalmip we can easily
solve Lyapunov equation
  
    A'P + P * A <0, P> 0

Is it possible to solve this equation in the Maple?


Regards,

 

hello

we have an exam next week and I want to know

how I can write (fordo) in maple for numerical integration

in different methods such as trapezoid , newton cotes and so on.

thanks

Is it possible to solve piecewise differential equations directly instead of separating the pieces and solving them separately.

like for example if i have a two dimensional function f(t,x) whose dynamics is as follows:

dynamics:= piecewise((t,x) in D1, pde1, pde2); where D1 is some region in (t,x)-plane

now is it possible to solve this system with one pde call numerically?

pde(dynamics, boundary conditions, numeric); doesnot work

This integral is defined as

Where

 

Now, How can I solve this integral?

tnx for your help.

 

 

GroupTheory...

July 12 2014 Kurt Ewald 40
1 1

G with Order 12 acts on the Set 1.2...12.

How can I get an isomorphic group which acts on the Set 1.2...4  ?

Kurt Ewald

 

 

With Maplets, it is possible to create popup windows to tell or ask the user something, but only when it is appropriate. That is, it is not there all the time. A particular usage would be to give the user a facilitating question when they have incorrect entries.

Is that possible in some way with Möbius apps also? Or does one have to always keep a blank TextArea available, just a blank space if invisible perhaps but using up space in the presentation nonetheless, to give the student the information?

i have read data from a text file. i want to convert this data to array,how should i do that ?

for example :
 i have readed these files and i want convert them to array which for the first one,has 3 rows and 4 columns.
what should i do ?

restart:

FileTools[Text][ReadFile]("NLIST.txt");

"       1     0.00000000000       0.00000000000       0.00000000000    
       2     50.0000000000       0.00000000000       0.00000000000    
       3     25.0000000000       0.00000000000       0.00000000000
        

"

(1)

op(FileTools[Text][ReadFile]("ELIST.txt"));

"       1   1   1   1   0   1      1     3     0
       2   1   1   1   0   1      3     2     0

      
      "

(2)

 

NULL


Download convert_to_array.mw   ELIST.txt  NLIST.txt

4 5 6 7 8 9 10 Last Page 6 of 1031