Dear all

I tried to compute an integral using trapezoidal rule using two differents methods but I get two different results
trap.mw

Thank you for your help 

How to solve and plot the integral equation?
b.mw
 

I am trying to input X-data, Y-data from spreadsheet-section of Numbers (Apple App) to do reression analysis in Maple. I do not understand the export or import tools very well just yet, asis obvious from numerous wrong results.

Since I want to do regression of the X-Y data trends, what format would you recommend using for the data once imported into Maple; two sequences, two-vectors, multi-row column matrix ???

If someone has worked with Numbers-Apple-App (Not Excel) into Maple to do this and can provide simple example with steps/tools used. It would be appreciated. Thanks

Dear all

I would like to sketch the solution of first ode, but by hand we can not get the solution. I used maple, to  solve the ode but unfortunattely there is no solution obtained so I can sketech the behavior of the solution. 

Is there any way to sketch the solution behavior .

sketch_solution.mw

thank you 

Hi. I have solved numerical an ODE and want to use the solution in the integrand of a subsequent integral. I'm being frustrated by my lack of knowledge and the internet isn't helping! I hope someone can put me right.....here's the relevant part of code.

dsysA := {diff(r(s), s) = sqrt((g*r(s) + h)^2 - 1 + 1/r(s)), r(0) = rinit};
dsnA := dsolve(dsysA, numeric, abserr = 10^(-6));

dsnA(0.1);
               [s = 0.1, r(s) = 1.51303353884965]

rA := s -> rhs(dsnA(s)[2]):

rA(0.2);
                        1.53235325247661

tA := s -> tinit + int((g*rA(p) + h)/(1 - 1/rA(p)), p = 0 .. s):

tA(0.01);
Error, (in rA) invalid input: rhs received dsnA(p)[2], which is not valid for its 1st argument, expr
 

Hi, I want to minimize a*x + log2(3+4.2^(-x)+2^(-2x)) , x>0 & 0<=a<=1 where a is some constant. Is there any command in maple which help me to minimize this. For command minimize, I am unable to minimize the function symbolically. I need to minimize it in terms of a.

Hi I am trying to compute the series approximation for the difference between these two divergent integrals

M := int(cosh(p)^2/sqrt(cosh(p)^2 - (1 + x)^2*sinh(p)^2), p = 0 .. 1/2*ln(1/x));

N := int((1 + x)*sinh(p)^2/sqrt(cosh(p)^2 - (1 + x)^2*sinh(p)^2), p = 0 .. 1/2*ln(1/x))

where x,p are positive and x is approaching zero. I would like to get a series expansion of M-N, but I am not quite if this is possible.

Dear Users!
Hope everyone fine here. I have a vector V and want a new vector after transformation like

V := Vector(5, {(1) = 2, (2) = 3, (3) = -1, (4) = 9, (5) = 7});

after square transformation of V I got a new vector W like bellow

W := Vector(5, {(1) = 4, (2) = 9, (3) = 1, (4) = 81, (5) = 49});

and similarly sin transformation of V give the following matrix as

U := Vector(5, {(1) = sin(2), (2) = sin(3), (3) = sin(-1), (4) = sin(9), (5) = sin(7)});
Please help me to get vectors W and U.

Dear all
Using cubic spline interpolation to interpolate a given function I get  
the follwoing error
Warning, solving for expressions other than names or functions is not recommended.
Error, (in solve) list of unknowns contains repeated elements: [1, 0]

and  why the  spline funciton obtained by this code is different to the one given by Maple's built in package.

inter.mw

Thank you for any help 

I want to remove the term XY
f:=(x,y)->2*x²+xy+y²+4x-y-2=0;
eq := simplify(subs(x = X-9/7, y = Y+8/7, f(x, y)));
theta:=Pi/8;
ex := simplify(subs(X = cos(theta)*X-sin(theta)*Y, Y = sin(theta)*X+cos(theta)*Y, eq)); evalf(%);
How to obtain (3+sqrt(2))/2*X²+(3-sqrt(2))/2=36 ? Thank you.

Hello,

I have two ordiary differential equations to be solved numerically. One of the equations can be solved independently (no need to solve the 2nd differential equation). However, in order to solve the 2nd differential equation, I need the numerical solution from the 1st one. 

sol1:=dsolve({ode1,x(0)=x__0, D(x)(0)=x__1},numeric, x(t), parameters=[a,b,c])

sol1(parameters=[4,5,1])

Now, I need to solve the 2nd differential equation which will give me y(t), but the differential equation depends on the function of x(t). What is the syntax of dsolve which incorporates the solution of another differential equation?

Hello, 

How one would define a function of a function-- a functional?

Given

f(g(r)) := int(g(r), r)

I would like to get 

int(h(r+eps),r)

when I call it like

f(h(r+eps))

I could achieve this behaviour by substitions, but is there a nicer and shorter way of doing this?

Thanks in advance.

I am trying to solve this equation, I tried
 

restart; 
ListTools[Categorize]; 
L := []; 
for a to 20 do
 for b to 20 do 
for c to 20 do 
for d to 20 do 
for e to 20 do 
for f to 20 do 
for g to 20 do 
for h to 20 do 
if sqrt(a+b*sqrt(c+d*sqrt(e +f*sqrt(g)))) = h then L := [op(L), [a, b, c, d, e, f, g, h]] end if
 end do end do end do end do end do end do end do end do; nops(L)

where g is not a square of a integer number . I do not get any solutions for a long time. How to reduce timing to solve this equation?

 

In this code, there are 8 vectors, each vector has 5 points except vectors 7 and 8 with 4 points. I made their  (7 and 8) respective fifth point a string. However, this error statement pops up "Error, invalid input: log[10] expects its 1st argument, x, to be of type algebraic, but received 2.51E-10".

Can someone help me with the correction?

Thank you all and kind regards

restart;

B:=<<601,1201,2401,4801,9601>|<2.87E-19,7.94E-21,7.94E-23,1.03E-24,5.91E-26>|
    <2001,4001,8001,16001,32001>|<3.30E-2,5.37E-4,8.47E-6,1.33E-7,2.08E-9>|
    <183,365,728,1476,2910>|<4.58E0,7.54E-8,2.20E-11,4.95E-14,9.15E-15>|
    <1621,3020,6166,12022,"2222">|<2.95E-3,8.51E-6,3.39E-8,2.51E-10,"2.51E-10">>:

for i from 1 to 5 do
   B[i, 2] := log[10](B[i, 2]):                      
   B[i, 4] := log[10](B[i, 4]):         
   B[i, 6] := log[10](B[i, 6]):         
   B[i, 8] := log[10](B[i, 8]):

   B[i, 1] := log[10](B[i, 1]):                      
   B[i, 3] := log[10](B[i, 3]):         
   B[i, 5] := log[10](B[i, 5]):         
   B[i, 7] := log[10](B[i, 7]):

end do:  # computing the log of the max-error
B: # This is the table of values we'll plot.

T:=plot([B[..,[1, 2]],B[1..1,[1, 2]], B[.., [3, 4]],B[1..1,[3, 4]], 
     B[..,[5, 6]],B[1..1,[5, 6]],B[.., [7, 8]],B[1..1,[7, 8]]], 
    legend = ["","BFFM","", "BNM","", "BHM","", "ARKN"],
    #title="Efficiency Curve for Example 5",
    style = ["pointline","point","pointline","point","pointline","point","pointline","point"], 
    symbolsize = 15,axes = framed, 
    symbol = [box,box, circle,circle,solidbox, solidbox,solidcircle, solidcircle],
    color=[ red, red, gold,gold, blue, blue,black, black], 
    axis = [gridlines = [colour = green, majorlines = 1,linestyle = dot]], 
    labels = [typeset(log__10(`NFE`)), typeset(log__10(`Max Err`))]);

This question came to me when I answered in this thread  https://www.mapleprimes.com/questions/231603-Solve-Onevariable-Equation

The equation there is quite cumbersome, I extracted a shorter subexpression from it (I converted all floats to exact constants), but the  solve  command hangs when trying to solve it. fsolve  handles the equation easily. We can see that the equation  Eq  is quite simple and easy to solve even by hand. I ask this question in a separate topic, because this seems to be a serious bug in the  solve  command.

restart;
Eq:=1-1/(1+203808*exp(-342569/506*t)*(1/131537))^(131537/203808)=44983/56599;
solve(Eq);

Edit. I noticed that the  isolate  command solves the problem, but of course the question remains open with  solve .