Unanswered Questions

This page lists MaplePrimes questions that have not yet received an answer

Hi 

Sometimes this happen to me.

why ?  and how can I fix it?

Please recover the content of the attached file if it is possible

Thanks

Check.mw

hi.i have a problem with gain Eigenvectors

please help me

thanks

lorenz_1.mw

restart; with(LinearAlgebra[Modular]); A := Matrix(3, 3, [[-a, a, 0], [0, 0, -sqrt(l.b.c.(h+k))/(h+k)], [(2.*h.sqrt(l.b.c.(h+k)))/(l.(h+k)), (2.*k.sqrt(l.b.c.(h+k)))/(l.(h+k)), -c]])

A := Matrix(3, 3, {(1, 1) = -a, (1, 2) = a, (1, 3) = 0, (2, 1) = 0, (2, 2) = 0, (2, 3) = -(`.`(l, b, c, h+k))^(1/2)/(h+k), (3, 1) = 2.*(h.((`.`(l, b, c, h+k))^(1/2)))/(l.(h+k)), (3, 2) = 2.*(k.((`.`(l, b, c, h+k))^(1/2)))/(l.(h+k)), (3, 3) = -c})

(1)

Sol := LinearAlgebra[Eigenvectors](A)

Error, (in evala/preproc3) floats not handled yet

 

``


Download lorenz_1.mw

DEAR SIR,

PLEASE HELP ME WITH THAT QUESTION

I have a gensym routine (gensym = generate symbol) which appends Chinese characters to some base symbols to create 'new' symbols.  This works well on Windows, Mac and tty maple.  But I get all square boxes in the linux GUI.

I know the fonts exist as, on the same machine, it works in tty.  But I have no idea how to tell the GUI to go and use those fonts.  Nor why it isn't.

Amusingly, if I try to paste something from the GUI as MapleMath here on Primes, in the 'paste' box, it shows just fine, but then it says it is invalid Maple.

For example, from the GUI I see

but the same thing in TTY is

(x三0不 + x三1下)^(-M)

and I just did cut-and-paste of the same thing.  By this I mean that I took that output with boxes above, put it into a terminal window, re-copied that, pasted it in here, and voila!

So clearly the issue is with the GUI, and only on linux (Ubuntu 16.04).

I watched a webinar training:   Maple Training for Engineers, Researchers, and Scientists.

 

The file demonstrated topics and used collapsible regions with a colored band across the document and a + / - sign on the left to open/close the section or region.  Different than the document gray down arrow.

 

I have searched all the documentation and can not find how to do this.   I have done this in other programs like MathCAD to close and lock a region from viewing or editing of specific formulas.

 

Any help always appreciated.

 

Regards,

Bill

I'm currently working on building a Grid Layout for a project, and I'm having trouble coding in the RunWindow and GetFile elements into buttons under the grid layout. I've gone through the overviews and examples for them, but had no luck. I'm using Maple 2016.1 for OS X.

Additionally, the structure of the code is slightly different as to how many of the example worksheets structure their Grid Layout code, since the code originated from a Maplet Builder file. I.e. in the example worksheets they would follow as:

maplet := Maplet('onstartup' = 'Action1', 'reference' = 'Maplet1',
         BoxLayout('background' = "#D6D3CE", 'border' = 'false', 'halign' = 'center', 'inset' = '5', 'reference' = 'BoxLayout1', 'valign' = 'center', 'vertical' = 'false', 'visible' = 'true',
                       BoxColumn( BoxCell('hscroll' = 'never', 'value' = 'Button1', 'vscroll' = 'never'),
         GridLayout('background' = "#D6D3CE", 'border' = 'false','halign'='center','inset'='5', 'reference' = 'GridLayout1', 'valign' = 'center', 'visible' = 'true',
                   GridRow('valign' = 'top', GridCell('height' = '1', 'hscroll' = 'never', 'value' = 'BoxLayout1', 'vscroll' = 'never', 'width' = '1' ))),
         Window('layout'= 'GridLayout1', 'reference' = 'W1', 'resizable' = 'true', 'title' = "Maplet"),
          Action('reference' = 'Action1', RunWindow('window'= 'W1'))

However the structure for the code I am working with has action at the very start of the code, follwed by the the code for the buttons then layouts/window.  E.g. (the code has been shortened)

with (Maplets[Elements]):
maplet :=
Maplet('onstartup'='Action1','reference'='Maplet1',
Action('reference'='clickButton1'),
Action('reference'='clickButton9',
Evaluate('function'='plot3d(x^2*cos(y),x = -1 .. 1,y = -2*Pi .. 2*Pi)','target'='Plotter1','waitforresult'='true')),
Action('reference'='clickButton11'),
Action('reference'='clickButton12'),
Action('reference'='clickButton10'),
Button('background'="#D6D3CE",'caption'="Insert Molecular Geometry",'enabled'='true','foreground'="#000000",'onclick'='clickButton1','reference'='Button1','visible'='true'),

....

BoxLayout('background'="#D6D3CE",'border'='false','halign'='center','inset'='5','reference'='BoxLayout1','valign'='center','vertical'='false','visible'='true',
BoxColumn(
BoxCell('hscroll'='never','value'='Button1','vscroll'='never'),
BoxCell('hscroll'='never','value'='BoxLayout2','vscroll'='never'),
BoxCell('hscroll'='never','value'='BoxLayout3','vscroll'='never'),
BoxCell('hscroll'='never','value'='BoxLayout9','vscroll'='never'),
BoxCell('hscroll'='never','value'='BoxLayout14','vscroll'='never')),
BoxColumn(
BoxCell('hscroll'='never','value'='Label3','vscroll'='never'),
BoxCell('hscroll'='never','value'='Plotter1','vscroll'='never'),
BoxCell('hscroll'='never','value'='Slider1','vscroll'='never'))),
GridLayout('background'="#D6D3CE",'border'='false','halign'='center','inset'='5','reference'='GridLayout1','valign'='center','visible'='true',
GridRow('valign'='top',
GridCell('height'='1','hscroll'='never','value'='BoxLayout1','vscroll'='never','width'='1'))),
Window('layout'='GridLayout1','reference'='Window1','resizable'='true','title'="Maplet"),
Action('reference'='Action1',
RunWindow('window'='Window1'))):

Maplets[Display](maplet);

 

If anyone would be able to provide an example of code or some guidance I could follow that would be greatly appreciated! 

MAPLE 15 on Windows 10

I have a file on which I can no longer open. This happened after my computer made me long on again after being gone for a little while. Maple was not running when I got back on. This is not unusal and I just clinck on the file and restart.

This time this did not work.

It comes up with a message asking me to select an input mode. Whatever I select even 'Plain text' it just hangs.

I'd send you an image, but I don't see how to do it. The button provide only allows for something with  url and the snip I made is on my desktop.

I need to get thiis worksheet back. I've put a lot of effort into it.

-Traruh

After running Maple in a shell file, I come up with this error that I do not understand on my Mac,

gap_long := 0.117647058823529 Pi

gap_lat := 0.0588235294117647 Pi

lat_begin := 0.441176470588235 Pi

long_begin := -Pi

lat_begin_0 := 0.441176470588235 Pi

long_begin_0 := -Pi

long_max := 0.882352941176471 Pi

lat_max := -0.441176470588235 Pi

33

Warning, `parameter` is implicitly declared local to procedure `set_par_eff`

distance eff distance_eff
im in has not
im in has not
im in has not
im in has not
im in has not
im in has not
im in has not
im in has not
im in has not
Im in has par
Im in has par
Error, invalid input: eval expects its 2nd argument, eqns, to be of type
{integer, equation, set(equation)}, but received par_eff_post
13

32

31

17

hou := 0

mini := 0

seci := 0

memory used=4.0MB, alloc=32.3MB, time=0.23



If needed, I can attach more files if my question is still a bit too cryptic. Please let me know asap as this is urgent. Thank you so much,
-Z

hi.please help me for solve this equation

i encounter with error''

Error, (in StringTools:-IsPrefix) second argument must be a string''

equations which be solved attached as pdf file

thanks

Kernel4.mw

root.pdf


 

restart

with(LinearAlgebra):

Typesetting:-Settings(functionassign=false):

NULL

Constants

 

landa := 0.404e11; -1; mu := 0.27e11; -1; alpha := 0.23e-4; -1; rho := 2707; -1; k := 204; -1; c := 903; -1; nu := .3; -1; E := 0.70e11; -1; T0 := 293; -1; omega := 0.1e-1

0.1e-1

(1.1.1)

beta := alpha*(3*landa+2*mu):

NULL

varpi := 0.1e-1; -1; No := 15

15

(1.1.2)

 

 

Eq[1] := besselj(0, xi*b)*(eval(diff(bessely(0, xi*r), r), r = a))-(eval(diff(besselj(0, xi*r), r), r = a))*bessely(0, xi*b):

 

wf1 := unapply(Eq[1], xi):

1

 

1.794010904

 

1

 

2

 

1.794010904

 

1

 

3

 

4.802060761

 

2

 

4

 

4.802060761

 

2

 

5

 

4.802060761

 

2

 

6

 

7.908961712

 

3

 

7

 

7.908961712

 

3

 

8

 

7.908961712

 

3

 

9

 

11.03509457

 

4

 

10

 

11.03509457

 

4

 

11

 

11.03509457

 

4

 

12

 

11.03509457

 

4

 

13

 

14.16798650

 

5

 

14

 

14.16798650

 

5

 

15

 

14.16798650

 

5

 

16

 

17.30400975

 

6

(1.2.1)

Eq[2] := MTM:-besselj(1, eta*b)*(eval(diff(MTM:-bessely(1, eta*r), r), r = a))-(eval(diff(MTM:-besselj(1, eta*r), r), r = a))*MTM:-bessely(1, eta*b):

wf2 := unapply(Eq[2], eta):

1

 

1.958510605

 

1

 

2

 

1.958510605

 

1

 

3

 

4.857021628

 

2

 

4

 

4.857021628

 

2

 

5

 

4.857021628

 

2

 

6

 

7.941288451

 

3

 

7

 

7.941288451

 

3

 

8

 

7.941288451

 

3

 

9

 

11.05802155

 

4

 

10

 

11.05802155

 

4

 

11

 

11.05802155

 

4

 

12

 

11.05802155

 

4

 

13

 

14.18576207

 

5

 

14

 

14.18576207

 

5

 

15

 

14.18576207

 

5

 

16

 

17.31852918

 

6

(1.2.2)

 

for m to MM do K0[m] := proc (r, m) options operator, arrow; BesselJ(0, xi[m]*r)*BesselY(0, xi[m]*b)-BesselJ(0, xi[m]*b)*BesselY(0, xi[m]*r) end proc; KK0[m] := proc (r, m) options operator, arrow; diff(K0[m](r, m), r) end proc; K1[n] := proc (r, n) options operator, arrow; BesselJ(1, eta__n*r)*BesselY(1, eta__n*b)-BesselJ(1, eta__n*b)*BesselY(1, eta__n*r) end proc; KK1[n] := proc (r, n) options operator, arrow; diff(K1[n](r, n), r) end proc end do

proc (r, m) options operator, arrow; BesselJ(0, xi[m]*r)*BesselY(0, xi[m]*b)-BesselJ(0, xi[m]*b)*BesselY(0, xi[m]*r) end proc

 

proc (r, m) options operator, arrow; MTM:-diff(K0[m](r, m), r) end proc

 

proc (r, n) options operator, arrow; BesselJ(1, eta__n*r)*BesselY(1, eta__n*b)-BesselJ(1, eta__n*b)*BesselY(1, eta__n*r) end proc

 

proc (r, n) options operator, arrow; MTM:-diff(K1[n](r, n), r) end proc

 

proc (r, m) options operator, arrow; BesselJ(0, xi[m]*r)*BesselY(0, xi[m]*b)-BesselJ(0, xi[m]*b)*BesselY(0, xi[m]*r) end proc

 

proc (r, m) options operator, arrow; MTM:-diff(K0[m](r, m), r) end proc

 

proc (r, n) options operator, arrow; BesselJ(1, eta__n*r)*BesselY(1, eta__n*b)-BesselJ(1, eta__n*b)*BesselY(1, eta__n*r) end proc

 

proc (r, n) options operator, arrow; MTM:-diff(K1[n](r, n), r) end proc

 

proc (r, m) options operator, arrow; BesselJ(0, xi[m]*r)*BesselY(0, xi[m]*b)-BesselJ(0, xi[m]*b)*BesselY(0, xi[m]*r) end proc

 

proc (r, m) options operator, arrow; MTM:-diff(K0[m](r, m), r) end proc

 

proc (r, n) options operator, arrow; BesselJ(1, eta__n*r)*BesselY(1, eta__n*b)-BesselJ(1, eta__n*b)*BesselY(1, eta__n*r) end proc

 

proc (r, n) options operator, arrow; MTM:-diff(K1[n](r, n), r) end proc

 

proc (r, m) options operator, arrow; BesselJ(0, xi[m]*r)*BesselY(0, xi[m]*b)-BesselJ(0, xi[m]*b)*BesselY(0, xi[m]*r) end proc

 

proc (r, m) options operator, arrow; MTM:-diff(K0[m](r, m), r) end proc

 

proc (r, n) options operator, arrow; BesselJ(1, eta__n*r)*BesselY(1, eta__n*b)-BesselJ(1, eta__n*b)*BesselY(1, eta__n*r) end proc

 

proc (r, n) options operator, arrow; MTM:-diff(K1[n](r, n), r) end proc

 

proc (r, m) options operator, arrow; BesselJ(0, xi[m]*r)*BesselY(0, xi[m]*b)-BesselJ(0, xi[m]*b)*BesselY(0, xi[m]*r) end proc

 

proc (r, m) options operator, arrow; MTM:-diff(K0[m](r, m), r) end proc

 

proc (r, n) options operator, arrow; BesselJ(1, eta__n*r)*BesselY(1, eta__n*b)-BesselJ(1, eta__n*b)*BesselY(1, eta__n*r) end proc

 

proc (r, n) options operator, arrow; MTM:-diff(K1[n](r, n), r) end proc

 

proc (r, m) options operator, arrow; BesselJ(0, xi[m]*r)*BesselY(0, xi[m]*b)-BesselJ(0, xi[m]*b)*BesselY(0, xi[m]*r) end proc

 

proc (r, m) options operator, arrow; MTM:-diff(K0[m](r, m), r) end proc

 

proc (r, n) options operator, arrow; BesselJ(1, eta__n*r)*BesselY(1, eta__n*b)-BesselJ(1, eta__n*b)*BesselY(1, eta__n*r) end proc

 

proc (r, n) options operator, arrow; MTM:-diff(K1[n](r, n), r) end proc

(1.2.3)

U1 := -(int(r*K0[m]*(diff(K1[n], r)+K1[n]/r), r = a .. b))/(int(r*K0[m]^2, r = a .. b)); -1; U2 := -(int(r*K1[n]*(diff(K0[m], r)), r = a .. b))/(int(r*K1[n]^2, r = a .. b)); -1; U3 := (int(r^2*omega^2*K1[n], r = a .. b))/(int(r*K1[n]^2, r = a .. b))

0.1555555555e-3/K1[0]

(1.2.4)

m := 0; -1; for m to MM do M__m := int(r*K1[m](r, m)^2, r = a .. b); bb__m := 1/M__m end do

int(r*K1[1](r, 1)^2, r = 1 .. 2)

 

1/(int(r*K1[1](r, 1)^2, r = 1 .. 2))

 

int(r*K1[2](r, 2)^2, r = 1 .. 2)

 

1/(int(r*K1[2](r, 2)^2, r = 1 .. 2))

 

int(r*K1[3](r, 3)^2, r = 1 .. 2)

 

1/(int(r*K1[3](r, 3)^2, r = 1 .. 2))

 

int(r*K1[4](r, 4)^2, r = 1 .. 2)

 

1/(int(r*K1[4](r, 4)^2, r = 1 .. 2))

 

int(r*K1[5](r, 5)^2, r = 1 .. 2)

 

1/(int(r*K1[5](r, 5)^2, r = 1 .. 2))

 

int(r*K1[6](r, 6)^2, r = 1 .. 2)

 

1/(int(r*K1[6](r, 6)^2, r = 1 .. 2))

(1.2.5)

MM; 1; n; 1; m; 1; U1; 1; U2; 1; U3; 1; xi

6

 

0

 

7

 

-(2/3)*K1[0]/K0[7]

 

0

 

0.1555555555e-3/K1[0]

 

xi

(1.2.6)

for m to MM do for n to MM do dsys := {diff(S[m][n](t), t, t, t)+xi^2*[m]*(diff(S[m][n](t), t, t))+(-U1*U2+`η__η__n__`^2)*(diff(S[m][n](t), t))+xi[m]^2*`η__η__n__`^2*S[m][n](t) = -(2*U2*bb[m]/(Pi*xi[m])*(-BesselJ(0, xi[m]*b)/BesselJ(1, xi[m]*a)))*q+xi^2*[m]*U3} end do end do; sol := dsolve(dsys)

Error, (in StringTools:-IsPrefix) second argument must be a string

 

 

NULL

for m to MM do for n to MM do dsys2 := {diff(Q__mn(t), t, t, t)+xi[m]^2*(diff(Q__mn(t), t, t))+(-U1*U2+eta__n^2)*(diff(Q__mn(t), t))+xi[m]^2*eta__n^2*Q__mn(t) = -2*BesselJ(0, xi[m]*b)*U1*U2*b__m*(1-exp(-xi[m]^2*t))/(BesselJ(1, xi[m]*a)*Pi*xi[m]^3)} end do end do;

sol2 := dsolve(dsys2)

Error, (in dsolve) invalid input: `PDEtools/sdsolve` expects its 1st argument, SYS, to be of type Or(set({`<>`, `=`, algebraic}), list({`<>`, `=`, algebraic}), `casesplit/ans`(list, list)), but received [{Q__mn(t)*pochhammer(1-n, n)+(1497143767/5000000)*(diff(Q__mn(t), [`$`(t, t)]))+eta__n^2*(diff(Q__mn(t), t))+(1497143767/5000000)*eta__n^2*Q__mn(t) = 0}]

 

``

NULL

NULL

 

Download Kernel4.mw

I'm tryinv to find the contents of specific bins in a maple histogram. I don't find this elementary 'function' anywhere in the help.

How do I do it? Surely it's possible w/o taking a ruler to the screen...

 

 

Hi, I want to ask. the maple program that i have done have something wrong somewhere.

for an example changes basis

 from e_{1}e_{1}=e_{1} , e_{1}e_{2}=e_{2} maple program reading :A1 : (1,1,1)=1,B1: (1,2,2)=1

to e_{2}e_{2}=e_{2} , e_{2}e_{1}=e_{1} maple program reading :A2: (2,2,2)=1,B2: (2,1,1)=1

these changes basis above are 2 operation, left product and right product

A1 and A2 are left product, while B1 and B2 are right product,

i need to make A1 isomorphic to A2, and B1 isomorphic to B2.

by using maple program, i should get identity in matrix form 2x2

[0 1] but i get [0         1]

[1 0],            [C_{21} 0],

For isomorphism, the determinant should not be zero

here's are my maple program:

>isom := proc (A1, A2, B1, B2, n)

local i, j, k, s, r, eqns, t, TEST, BChange, sols, m, S1, S2, C;

C := matrix(n, n);

BChange := matrix(n, n);

TEST := 0; eqns := {};

for i to n do for j to n do for m to n do

S1 := sum(A1[i, j, k]*C[k, m], k = 1 .. n); S2 := sum(C[i, r]*(sum(A2[r, s, m]*C[j, s], s = 1 .. n)), r = 1 .. n);

eqns := `union`(eqns, {S1 = S2})

end do end do end do;

for i to n do for j to n do for m to n do

S1 := sum(B1[i, j, k]*C[k, m], k = 1 .. n); S2 := sum(C[i, r]*(sum(C[j, s]*B2[r, s, m], s = 1 .. n)), r = 1 .. n);

eqns := `union`(eqns, {S1 = S2})

end do end do end do;

sols := [solve(eqns)];

t := nops(sols);

for i to t do for j to n do for k to n do

BChange[k, j] := subs(sols[i], C[k, j])

end do end do;

if simplify(linalg:-det(BChange)) <> 0 then print("BChange", BChange);

print("s1", S1); print("s2", S2); print("The det is", simplify(linalg:-det(BChange)));

TEST := 1 end if end do;

if TEST = 0 then print("These two algebras are not isomorphic")

end if end proc

input maple program:

> DENDA1 := array(sparse, 1 .. 2, 1 .. 2, 1 .. 2, [(1, 1, 1) = 1]);
> DENDB1 := array(sparse, 1 .. 2, 1 .. 2, 1 .. 2, [(1, 2, 2) = 1]);
> DENDA2 := array(sparse, 1 .. 2, 1 .. 2, 1 .. 2, [(2, 2, 2) = 1]);
> DENDB2 := array(sparse, 1 .. 2, 1 .. 2, 1 .. 2, [(2, 1, 1) = 1]);
> isom(DENDA1, DENDA2, DENDB1, DENDB2, 2);

Couldn't any question with title beginning with http be removed automatically? There has been quite a few containing nothing but spam.

Dear All,

I am going to solve the following systems of ODEs but get the error: Newton iteration is not converging.
Could you please share your idea with me. In the case of AA=-0.2,0,0.2,0.4,...; I could get the solution.
Thank you in advance.


restart;
with(plots);
Pr := 2; Le := 2; nn := 2; Nb := .1; Nt := .1; QQ := .1; SS := .1; BB := .1; CC := .1; Ec := .1; MM := .2;AA:=-0.4;

Eq1 := diff(f(eta), `$`(eta, 3))+f(eta).(diff(f(eta), `$`(eta, 2)))-2.*nn/(nn+1).((diff(f(eta), eta))^2)-MM.(diff(f(eta), eta)) = 0; Eq2 := 1/Pr.(diff(theta(eta), `$`(eta, 2)))+f(eta).(diff(theta(eta), eta))-4.*nn/(nn+1).(diff(f(eta), eta)).theta(eta)+Nb.(diff(theta(eta), eta)).(diff(h(eta), eta))+Nt.((diff(theta(eta), eta))^2)+Ec.((diff(f(eta), `$`(eta, 2)))^2)-QQ.theta(eta) = 0;
Eq3 := diff(h(eta), `$`(eta, 2))+Le.f(eta).(diff(h(eta), eta))+Nt/Nb.(diff(theta(eta), `$`(eta, 2))) = 0;

bcs := f(0) = SS, (D(f))(0) = 1+AA.((D@@2)(f))(0), theta(0) = 1+BB.(D(theta))(0), phi(0) = 1+CC.(D(phi))(0), (D(f))(etainf) = 0, theta(etainf) = 0, phi(etainf) = 0

Error, (in dsolve/numeric/ComputeSolution) Newton iteration is not converging

how to translate python code which use scipy, numpy to maple code

 

 

import numpy as np
from scipy.sparse.linalg import svds
from functools import partial


def emsvd(Y, k=None, tol=1E-3, maxiter=None):
    """
    Approximate SVD on data with missing values via expectation-maximization

    Inputs:
    -----------
    Y:          (nobs, ndim) data matrix, missing values denoted by NaN/Inf
    k:          number of singular values/vectors to find (default: k=ndim)
    tol:        convergence tolerance on change in trace norm
    maxiter:    maximum number of EM steps to perform (default: no limit)

    Returns:
    -----------
    Y_hat:      (nobs, ndim) reconstructed data matrix
    mu_hat:     (ndim,) estimated column means for reconstructed data
    U, s, Vt:   singular values and vectors (see np.linalg.svd and 
                scipy.sparse.linalg.svds for details)
    """

    if k is None:
        svdmethod = partial(np.linalg.svd, full_matrices=False)
    else:
        svdmethod = partial(svds, k=k)
    if maxiter is None:
        maxiter = np.inf

    # initialize the missing values to their respective column means
    mu_hat = np.nanmean(Y, axis=0, keepdims=1)
    valid = np.isfinite(Y)
    Y_hat = np.where(valid, Y, mu_hat)

    halt = False
    ii = 1
    v_prev = 0

    while not halt:

        # SVD on filled-in data
        U, s, Vt = svdmethod(Y_hat - mu_hat)

        # impute missing values
        Y_hat[~valid] = (U.dot(np.diag(s)).dot(Vt) + mu_hat)[~valid]

        # update bias parameter
        mu_hat = Y_hat.mean(axis=0, keepdims=1)

        # test convergence using relative change in trace norm
        v = s.sum()
        if ii >= maxiter or ((v - v_prev) / v_prev) < tol:
            halt = True
        ii += 1
        v_prev = v

    return Y_hat, mu_hat, U, s, Vt

     It is known that ODE boundary value problem is similar to the problem of solving systems of nonlinear equations. Equations are the boundary conditions, and the variables are the values of the initial data.
For example:

y '' = f (x, y, y '), 0 <= x <= 1,

y (0) = Y0, y (1) = Y1;

Where y (1) = Y1 is the equation, and Z0 is variable, (y '(0) = Z0).

     solve () and fsolve () are not directly suitable for such tasks. Directly should work the package of optimization in relation to a system of nonlinear equations. (Perhaps it has already been implemented in Maple.)
Personally, I am very small and unprofessional know Maple and cannot do it. Maybe there is someone who would be interested, and it will try to implement this approach to solving ODE boundary value problems?  

First 186 187 188 189 190 191 192 Last Page 188 of 362