mehdi jafari

749 Reputation

13 Badges

12 years, 108 days

MaplePrimes Activity


These are questions asked by mehdi jafari

hi everyone , i need your point of view in my question,any help would be appreciated in advance .

we have a discrete function named g(t) and a continous function f(t) in in convolution integral just like this :
int(f(t-x)*g(x),x=0..t) ; 
we have just g(x) in some special points int the interval (0..t) , thus i need to convert this integral to a series.
how should i do this ? can anyone help or any idea ? i need at first a mathematical solution or idea about how to do this and then, how to do this in a software ?
tnx again.

i want to define a loop for like this : 

for i from 0 to 100 by 0.2 do; ... ; od ;

i want when this loop finishes,the the loop turns to exactly vice versa respect,and continue , i mean it continue like this :
for i from 100 to 0 by -0.2 ...  , 

can anybody help me to find an algorithm for this ? tnx in advance


``

restart:

l := 16; m := .75; EIy := 0.2e5; EIz := 0.400e7; GJ := 0.1e5; mj := .1; Pz := 5000;

16

 

.75

 

0.2e5

 

0.400e7

 

0.1e5

 

.1

 

5000

(1)

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;

.75*(diff(diff(w(x, t), t), t))+0.2e5*(diff(diff(diff(diff(w(x, t), x), x), x), x))+0.398e7*(diff(diff(diff(diff(w(x, t), x), x), x), x))*(theta(x)+phi(x, t))^2+0.1592e8*(diff(diff(diff(w(x, t), x), x), x))*(theta(x)+phi(x, t))*(diff(theta(x), x)+diff(phi(x, t), x))+0.796e7*(diff(diff(w(x, t), x), x))*(diff(theta(x), x)+diff(phi(x, t), x))^2+0.796e7*(diff(diff(w(x, t), x), x))*(theta(x)+phi(x, t))*(diff(diff(theta(x), x), x)+diff(diff(phi(x, t), x), x))-5000 = 0

(2)

w := unapply(q[1](t)*cosh(1.8751*x/l),x, t);

proc (x, t) options operator, arrow; q[1](t)*cosh(.1171937500*x) end proc

(3)

theta := proc (x) options operator, arrow; 0 end proc;

proc (x) options operator, arrow; 0 end proc

(4)

#theta(x);eq1;

phi := unapply(q[4](t)*sin(1.5708*x/l),x, t);

proc (x, t) options operator, arrow; q[4](t)*sin(0.9817500000e-1*x) end proc

(5)

eq2:=unapply(lhs(eq1),x,t);

proc (x, t) options operator, arrow; .75*(diff(diff(q[1](t), t), t))*cosh(.1171937500*x)+3.772661154*q[1](t)*cosh(.1171937500*x)-302.9569534*q[1](t)*cosh(.1171937500*x)*q[4](t)^2*sin(0.9817500000e-1*x)^2+2515.691179*q[1](t)*sinh(.1171937500*x)*q[4](t)^2*sin(0.9817500000e-1*x)*cos(0.9817500000e-1*x)+1053.716523*q[1](t)*cosh(.1171937500*x)*q[4](t)^2*cos(0.9817500000e-1*x)^2-5000 end proc

(6)

f[1] := int(expand(eq2(x,t)*cosh(1.8751*x/l)), x = 0 .. l);

int(.75*(diff(diff(q[1](t), t), t))*cosh(.1171937500*x)^2+3.772661154*q[1](t)*cosh(.1171937500*x)^2-302.9569534*q[1](t)*cosh(.1171937500*x)^2*q[4](t)^2*sin(0.9817500000e-1*x)^2+2515.691179*cosh(.1171937500*x)*q[1](t)*sinh(.1171937500*x)*q[4](t)^2*sin(0.9817500000e-1*x)*cos(0.9817500000e-1*x)+1053.716523*q[1](t)*cosh(.1171937500*x)^2*q[4](t)^2*cos(0.9817500000e-1*x)^2-5000*cosh(.1171937500*x), x = 0 .. 16)

(7)

 

 

``


Download 1.mw

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

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

First 16 17 18 19 20 21 22 Page 18 of 25