Iza

30 Reputation

3 Badges

3 years, 224 days

MaplePrimes Activity


These are questions asked by Iza

Hello,

I will explain my problem on a simple function and some data.

mg:=78.54*7.85*1e-9*1000*9.81;l:=10000;l0:=12000;h:=0;H := 28.399;
eta:=mg*l/(2*H);
zeta:=arcsinh(mg*h/(2*H*sinh(eta)))-eta;
z:=H/mg*(cosh(mg/H*x+zeta)-cosh(zeta));
plot(z,x=0..l);

Above is a definition of an initial plot. Now I would like to draw a new plot in which every x coordinate is shifted horizontally by a function:

deltax:=H*arcsinh(x)-x;

As you see I can't use translate or scale so I've tried to use a transform function but I can't do it the right way.

I would be grateful for any ideas.

 

Greetings,

Iza

Hello,

I am solving a simple system of algebraic equations in a "for" loop. I would like to write all the values of solutions to a matrix, but I only get the symbols. Here is my code:

A:=1:E:=0.1:L0:=<5,5,5,5>:m0g:=0.05:u:=<-5/4,-1/4,-3/4,-1/4>:v:=<-1/2,-3/2,1/2,1/2>:w:=<0,1,0,-1>:
Force:=Matrix(m,2);
for i from 1 to m do
 r1||i:=H||i*L0(i)/(A*E)+H||i/m0g*(arcsinh(V_A||i/H||i)+arcsinh((m0g*L0(i)-V_A||i)/H||i))-sqrt(u(i)^2+v(i)^2);
 r2||i:=V_A||i*L0(i)/(A*E)*(m0g*L0(i)/(2*V_A||i)-1)+H||i/m0g*(sqrt(1+((m0g*L0(i)-V_A||i)/H||i)^2)-sqrt(1+(V_A||i/H||i)^2))-abs(w(i));
 rozw||i:=fsolve({r1||i, r2||i}, {H||i,V_A||i},H||i=0..1);Force(i,1):=H||i;Force(i,2):=V_A||i;subs(rozw||i,Force(i,..));
end do:
evalf(Force);

Thanks in advance,

Iza

Hello,

I've tried to find the solution to my problem, but none of my attempts was succesful.

I have a function which is one-to-one in a particular domain which I am interested in. I would like to get the inverse function of it only in this domain. Here is my function and plot for xp=0..10000:

x := xp-> (-1)*720.5668720*sinh(0.2043018094e-3*xp-0.8532729286)+84952.59423+4.014460003*10^5*arcsinh(0.1219272144e-1*sinh(0.2043018094e-3*xp-0.8532729286)-0.2032498888)

I would appreciate any help,

Iza

Hello,

I use a map() command to get the values of function for each element of vector. Here is an example of a simple task:

restart:

A:=[1,2,8];

f1:=x->2*x;
f1_table_A:=Vector([map(f1,A)]);

I would like to do the same with piecewise function but it doesn't work. Here is an example:

restart:
with(plots):
f:=2*x:
g:=x^2:

h:=x->piecewise( 0<=x and x<= 5, f, 5<x and x <= 10, g ) ;
A:=[1,2,8]:
h_table_A:=map(h,A):

Is it possible to work with piecewise function? Maybe someone has an idea how to do it in different way?

Kind regards,

Iza

Hallo,

I am rather a beginner in Maple, so I will be grateful for understanding.

I try to solve the equation with one unknown. To create this equation I firstly had to solve the set of four equations so it is very complicated. Maple is able to create the plot of it but command fsolve is not sufficient to solve it. Here is my procedure:

restart:
z1:=H/mg*cosh(mg/H*(x+C11))-C21:
z2:=H/mg*cosh(mg/H*(x+C12))-C22:
tg_alpha1:=diff(z1,x):
tg_alpha2:=diff(z2,x):
r1:=H*subs(x=xP,tg_alpha2)-H*subs(x=xP,tg_alpha1)-P:
r2:=subs(x=0,z1):
r3:=subs(x=l,z2)-h:
r4:=subs(x=xP,z1)-subs(x=xP,z2):
mg:=2:l:=20:h:=5:xP:=l/5:P:=10:
C := solve({r1, r2, r3, r4}, {C11, C12, C21, C22}):
s1 := int(sqrt(1+(diff(subs(C, z1), x))^2), x):
s2 := int(sqrt(1+(diff(subs(C, z2), x))^2), x):
s := subs(x = xP, s1)-subs(x = 0, s1)+s2-subs(x = xP, s2):
eq1 := l0-subs(x = l, s):
l0 := 40:
plot(eq1, H=10 .. 11, colour = red);
H := fsolve(eq1 = 0, H = 10.7 .. 10.8);

I would be grateful for any help and indicating my mistakes.

Thanks

Iza

Page 1 of 1