Advanced Search

Axel Vogt

5821 Reputation

20 Badges

20 years, 221 days
Munich, Bavaria, Germany
Contact Axel Vogt

MaplePrimes Activity


These are replies submitted by Axel Vogt

Cengel2_1.mw...

Axel Vogt 5821
December 18 2022
0 0

This is what I get opening Cengel2.mw using Maple 2021.2 (which asks for repairing it and warns me about auto execution) and save it as Cengel2_1.mw

Cengel2_1.mw

like this?...

Axel Vogt 5821
December 17 2022
0 0
> 

A:=<1,2;3,4>: B:=<6,7;8,9>:
> 

'A.B': '%'=%;

A.B = (Matrix(2, 2, {(1, 1) = 22, (1, 2) = 25, (2, 1) = 50, (2, 2) = 57}))

 (1)
> 

 

Download MP_235397.mw

thx for correcting me...

Axel Vogt 5821
December 11 2022
0 0

thx for correcting me

"break"...

Axel Vogt 5821
December 11 2022
0 0

Look up the help: the command is just "break", not "break something".

Maple 2021 Standard on Win 7...

Axel Vogt 5821
December 08 2022
0 0

I can work with it. But I am aware that the Win version is too old. And, Christopher, it may depend on specific configurations of your machine as well, especially concerning GUI / Java. Personally I would not try it ... for example you can not be sure that it will help for the next Maple version.

and without output?...

Axel Vogt 5821
December 04 2022
0 0

And what happens if you suppress the output (using ":" instead of ";") ?

here you go ......

Axel Vogt 5821
November 19 2022
0 0

Here it is enough that it is divergent in x=0 if you "know" about Int(1/ln(x), ...)
 

> 

# https://www.mapleprimes.com/questions/235198--How-To-Consider-The-Convergence-Of
> 

restart: interface(version);

`Standard Worksheet Interface, Maple 2021.2, Windows 7, November 23 2021 Build ID 1576349`

 (1)
> 

g:= x -> (x+2)/(2*x^3+1);
plot(g(x), x=1 ..2);

proc (x) options operator, arrow; (x+2)/(2*x^3+1) end proc

  

 
> 

g(2)/ln(x) <= 'g(x)/ln(x)'; is(%) assuming 1<x,x<=2;
[g(2)/ln(x) , g(x)/ln(x)]: plot(%, x=1..2, color=[red, blue]);

(4/17)/ln(x) <= g(x)/ln(x)

  

true

  

 
> 

 


 

Download MP_235198.mw

It was possible in former versions, for ...

Axel Vogt 5821
October 31 2022

It was possible in former versions, for example 2017, no?

file attached...

Axel Vogt 5821
October 12 2022
0 0

@golnaz Sorry for being too thin-lipped, the following should make it clear

NULL

> 

# https://www.mapleprimes.com/questions/234977-Why-Fsolve-Doesnt-Work
restart:
> 

N := 60: L := 10: a := 2*10^(-12):
> 

eqn1:=ns=1-((((N/A)+((A*L)^(L/(1-L))))^((1-L)/L))/(A*L))*(2+(a*exp((((N/A)+((A*L)^(L/(1-L)))))^(1/L)))/((A*L*(((N/A)+((A*L)^(L/(1-L))))^((L-1)/L)))+((a*exp((((N/A)+((A*L)^(L/(1-L)))))^(1/L)))))):
> 

#rhs(eqn1): plot(%, A=0..10); MultiSeries:-series(%%, A=0): evalf(%);
> 

eval(eqn1, A= 1/B^9): Eqn1:=simplify(%) assuming positive:
[seq(fsolve((Eqn1)), ns=[0.9603,0.9647,0.9691])]:
map(x -> 1/x^9, %);

[0.930068277335412e-9, 0.305414431895731e-8, 0.117266280509612e-7]

 (1)
> 

#plot(rhs(Eqn1));
> 

 

> 

eqn2:=r=16*(((A*L*(((N/A)+((A*L)^(L/(1-L))))^((L-1)/L)))+((a*exp((((N/A)+((A*L)^(L/(1-L)))))^(1/L))))))/(((A*L*(((N/A)+((A*L)^(L/(1-L))))^((L-1)/L))))^2):
> 

#rhs(eqn2): plot(%, A=0..10); MultiSeries:-series(%%, A=0): evalf(%);
> 

eval(eqn2, A= 1/B^9): Eqn2:=simplify(%) assuming positive:
[seq(fsolve((Eqn2)), r=[0.2,0.1])]:
map(x -> 1/x^9, %);

[0.994026002208275e-7, 0.105660856235243e-3]

 (2)
> 

#plot(rhs(Eqn2));
> 

 

NULL

 

Download MP_234977.mw

"smooth ?"...

Axel Vogt 5821
October 06 2022
0 0

I do not know (and what curves you have in mind), you have to study your integrand, essentially that question is like asking for the continuity of the following second term (t ~ your afa, F ~ your integrand depending on parameters)


 

> 

Int(Int(F(x, y, t), x=0..1), y=0 .. t): diff(%, t);

Int(Int(diff(F(x, y, t), t), x = 0 .. 1), y = 0 .. t)+Int(F(x, t, t), x = 0 .. 1)

 (1)
> 

 

 

another plot...

Axel Vogt 5821
October 06 2022
0 0

'sin(x)*T2'; #log10('%');
plot3d(%, x=Pi/2 - 0.05 .. Pi/2.0, y=-0.05 ..0.05, axes=boxed, numpoints=100^2);

'sin(x)*T2'; #log10('%');
plot3d(%, x=evalf(Pi/2) - 0.05 .. evalf(Pi/2), y=-0.05 ..0.05, view= 10.0^8 .. 10.0^10, axes=boxed, numpoints=100^2);

You can use the mouse for the plot to see, the ugly cases are roughly on a parabola (and that does not depend on adding some small imaginary part to gama).

A more careful study *may* give: outside that integration works and estimating the ugly region may contribute only a small portion (using vv's suggestion?).

But perhaps it would make sense you de-compose the task to the steps from where you built it

H(0) = 1/2...

Axel Vogt 5821
October 02 2022
0 0

@vv  yes, you are right

where does it stem from?...

Axel Vogt 5821
October 01 2022
0 0

May be it helps a bit if you can say how one achieves the integral (some volume?)

just "execute" the file...

Axel Vogt 5821
September 26 2022
0 0

@fhn2022 

Just execute the file, I attached it.

Either step by step  or  using the !!! symbol in the menu bar  or  through the menu execute / worksheet.

Re Question on ' ' and " "...

Axel Vogt 5821
September 13 2022
0 0

@Dkunb Here I just use it for a pretty presentation, you may look up the help on quotes and Unevaluated Expressions, 'expr'

For your question on N see acer's post  first image

5 6 7 8 9 10 11 Last Page 7 of 207