Mariusz Iwaniuk

1381 Reputation

14 Badges

7 years, 115 days

Social Networks and Content at Maplesoft.com

MaplePrimes Activity


These are answers submitted by Mariusz Iwaniuk

If we use MMA to solve series:

inttrans:-fourier(sech(x), x, -k);

#Pi*sech(k*Pi/2)

 

L := 1;
g := 9.81;
f := 4*sqrt(L/g)*EllipticF(Pi/2, sin(theta0/2));
plot([Re(f), Im(f)], theta0 = 1 .. 20);

With Maple 2022.2 and the option  adaptive I have:

f := proc (x) options operator, arrow; arcsin(2*x/(1+x^2)) end proc

plot(f(x), x, adaptive); plot(f(x), x, view = [-10 .. 10, -10 .. 10], adaptive); plot(f(x), x = -10 .. 10, adaptive); plot(f(x), x, 'adaptive' = true, 'numpoints' = 20)

 

 

 

 

help("Plot Adaptive Options#Execue this command")

Download ASin.mw

value(sol);
simplify(%);

I executed on Maple 2022.2 ,give me:

#                  1      /               /
theta__1(y) = ---------- |27 (y + sigma) |
                      12 \               \
              20 sigma                    
  5                     9   /  10                 
- - lambda (k - 1) sigma  + |- -- lambda (k + 1) y
  9                         \  27                 

     /  149     149\  4   /49      49 \  3   /  29     29\  2
   + |- --- n + ---| k  + |--- n - ---| k  + |- -- n + --| k 
     \  432     432/      \108     108/      \  72     72/   

     /49      49 \     10            149     149\      8   11   
   + |--- n - ---| k - -- Q lambda - --- n + ---| sigma  - -- (k
     \108     108/     9             432     432/          6    

        /10  2          5          / 2   14      \          
   - 1) |-- y  lambda + -- (k + 1) |k  - -- k + 1| (n - 1) y
        \33             36         \     11      /          

               / 2   2       \\      7   ///  1     1\  4
   + Q (n - 1) |k  - -- k + 1|| sigma  + |||- - n + -| k 
               \     11      //          \\\  8     8/   

     /1     1\  3   /  3     3\  2   /1     1\     1  
   + |- n - -| k  + |- - n + -| k  + |- n - -| k - - n
     \2     2/      \  4     4/      \2     2/     8  

     10            1\  2
   - -- Q lambda + -| y 
     9             8/   

     145                   / 2   50      \  
   - --- (k + 1) Q (n - 1) |k  - -- k + 1| y
     108                   \     29      /  

     23  2 / 2   26      \        \      6                     /
   - -- Q  |k  - -- k + 1| (n - 1)| sigma  - 4 (k - 1) (n - 1) |
     6     \     23      /        /                            \
  1                 2  3   1          2  2   5  2              3\ 
- -- (k + 1) (k - 1)  y  + - Q (k - 1)  y  + - Q  (k + 1) y + Q | 
  32                       4                 8                  / 

       5             /  1         4  4   3                  2  3
  sigma  - 2 (n - 1) |- -- (k - 1)  y  - - Q (k + 1) (k - 1)  y 
                     \  32               8                      

     3  2        2  2   5  3              4\      4               
   + - Q  (k - 1)  y  + - Q  (k + 1) y + Q | sigma  - 4 (k - 1) Q 
     2                  6                  /                      

          /  1        2  2   3                2\  2      3
  (n - 1) |- - (k - 1)  y  - - Q (k + 1) y + Q | y  sigma 
          \  8               8                 /          

        2          2 /  3        2  2   1                2\ 
   - 2 Q  (n - 1) y  |- - (k - 1)  y  - - Q (k + 1) y + Q | 
                     \  4               2                 / 

       2      3  4                          4  4        \   
  sigma  + 2 Q  y  (n - 1) (k - 1) sigma + Q  y  (n - 1)| (y
                                                        /   

              \
   - sigma) Br|
              /

Analytic_result_Help_ver2.mw

Using:

restart;
kernelopts(version);
#`Maple 2022.1, X86 64 WINDOWS, May 26 2022, Build ID 1619613`

A := factor(sum(sum(sum(sum(sum(sum(sum(sum(((x - i1)^2 + (y - i2)^2)*((x - j1)^2 + (y - j2)^2)*((x - k1)^2 + (y - k2)^2)*((x - l1)^2 + (y - l2)^2), i1 = 1 .. N), j1 = 1 .. N), k1 = 1 .. N), l1 = 1 .. N), i2 = 1 .. N), j2 = 1 .. N), k2 = 1 .. N), l2 = 1 .. N));

N := 8:
B := CodeTools:-Usage(A);

#memory used=1.77KiB, alloc change=0 bytes, cpu time=0ns, real time=0ns, gc time=0ns
#B := (5764801*(3*x^2 + 3*y^2 - 24*x - 24*y + 120)^4)/81

 

plot([Re((-2)^x), Im((-2)^x)], x = -2 .. 2, legend = ["Real", "Imaginary"]);

 

kernelopts(version)
#`Maple 2022.1, X86 64 WINDOWS, May 26 2022, Build ID 1619613`

(int(sin(x)*exp(-2*I*Pi*f*x)/x, x = -infinity .. infinity) assuming (0 < f));

#signum(0, -2*Pi*f + 1, 0)*Pi/4 + Pi/2 - signum(0, 2*Pi*f - 1, 0)*Pi/4

convert(%, piecewise);

#piecewise(f < 1/(2*Pi), Pi, f = 1/(2*Pi), Pi/2, 1/(2*Pi) < f, 0)

 

Order := 50;
ode := diff(x(u), u, u) + 1/4*exp(-u^2)*x(u) = 0;
sol := convert(dsolve([ode, x(0) = 1, D(x)(0) = 0], x(u), 'series'), polynom);
plot([rhs(sol), diff(rhs(sol), u)], u = 0 .. 2, view = [0 .. 2, -1/2 .. 1])

 

Looks like a bug in integrate.

int(x^n*exp(x), x = 0 .. 2, method = _RETURNVERBOSE); #Can't solve.

I have only a workround.

See attached file:2_ver1.mw

As a infinite sum:

Int(sin(sqrt(-x^2 + 1)), x) = Sum(2^(-1/2 - k)*sqrt(Pi)*x^(1 + 2*k)*BesselY(-1/2 + k, 1)/((1 + 2*k)*GAMMA(1 + k)), k = 0 .. infinity) + C

for: -1<=x<=1

Check:

int(sin(sqrt(-x^2 + 1)), x = 0 .. 1, numeric);
evalf(add(eval(2^(-1/2 - k)*sqrt(Pi)*x^(1 + 2*k)*BesselY(-1/2 + k, 1)/((1 + 2*k)*GAMMA(1 + k)), x = 1), k = 0 .. 2000));

 


 

restart

interface(rtablesize = 100)

eqns := {x(0) = 1, x(t+1) = (1+10^(-6)*(0.4e-2-0.6e-2)/(0.1e-3))*x(t)+0.7e-1*10^(-6)*y(t), y(0) = 857.1428571, y(t+1) = 0.6e-2*10^(-6)*x(t)/(0.1e-3)+(1-0.7e-1*10^(-6))*y(t)}

{x(0) = 1, x(t+1) = .9999800000*x(t)+0.7000000000e-7*y(t), y(0) = 857.1428571, y(t+1) = 0.6000000000e-4*x(t)+.9999999300*y(t)}

(1)

SOL := evalf(rsolve(convert(eqns, rational), {x, y}))

{x(t) = 2.959071538*1.000000138^t-1.959071539*.9999797913^t, y(t) = 851.3060776*1.000000138^t+5.836778775*.9999797913^t}

(2)

seq([rhs(SOL[1]), rhs(SOL[2])], t = 0 .. 20)

[.999999999, 857.1428564], [1.000039997, 857.1428559], [1.000079996, 857.1428555], [1.000119993, 857.1428549], [1.000159988, 857.1428545], [1.000199984, 857.1428540], [1.000239979, 857.1428536], [1.000279972, 857.1428531], [1.000319966, 857.1428526], [1.000359958, 857.1428522], [1.000399950, 857.1428517], [1.000439940, 857.1428513], [1.000479930, 857.1428509], [1.000519919, 857.1428504], [1.000559907, 857.1428499], [1.000599894, 857.1428495], [1.000639882, 857.1428491], [1.000679867, 857.1428487], [1.000719852, 857.1428482], [1.000759837, 857.1428478], [1.000799820, 857.1428474]

(3)

````

rtable([%], subtype = Vector[column])

Vector[column](%id = 36893490192963129204)

(4)

NULL


 

Download SystemRecursive_new2.mw

.

f := (x, y) -> piecewise(x <> 0 and y <> 0, (y^2 + x^2)^x, x = 0 and y = 0, 1);
DX := diff(f(x, y), x);
DY := diff(f(x, y), y);
[eval(DX, [x = 1, y = 1]), eval(DY, [x = 1, y = 1])]

#at x=1 and y=1;
#[2*ln(2) + 2, 2]

 


 

NULL

Vector(3, {(1) = (l+s(t))*cos(beta(t))+xo+s(t)-x(t), (2) = sin(beta(t))*(l+s(t))*sin(alpha(t))-y(t), (3) = -sin(beta(t))*(l+s(t))*cos(alpha(t))-z(t)})

Vector[column](%id = 36893489901319436636)

(1)

"solve(?[2],[beta(t)])[][]"

beta(t) = arcsin(y(t)/(sin(alpha(t))*(l+s(t))))

(2)

"subs(beta(t) = arcsin(y(t)/sin(alpha(t))/(l+s(t))),?[3])"

-sin(arcsin(y(t)/(sin(alpha(t))*(l+s(t)))))*(l+s(t))*cos(alpha(t))-z(t)

(3)

solve(-sin(arcsin(y(t)/(sin(alpha(t))*(l+s(t)))))*(l+s(t))*cos(alpha(t))-z(t), [alpha(t)])

[[alpha(t) = -arctan(y(t)/z(t))]]

(4)

" solve(simplify(subs(beta(t) = arcsin(y(t)/sin(alpha(t))/(l+s(t))),?[3])),[alpha(t)])"

[[alpha(t) = -arctan(y(t)/z(t))]]

(5)

NULL

NULL


 

Download Arctan_ver2.mw

 

restart;
with(LinearAlgebra);
w := (2*Pi)/14;
v := Vector([1, sin(w*t), cos(w*t)]);
simplify(eval(sum(v . (Transpose(v)), t = k .. k + n), n = 13));

 

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