220 Reputation

17 years, 34 days
Kharkov
PhD in Numerical Optimization.

(+)...

I tested a couple of example. You are right!

evalf(Int) return the same values each time.

So my problem is solved. Regards.

I figured out.

int / Int...

Regarding symbolic calculation of int, look here:

evalf(int) does numerical integration!

Спасибо...

Thank you, I'll try to insight into the problem.

OK...

I perfectly undestand, that the different values appear due to roundoff errors. But why are the errors different? If an algorithm is deterministic, the errors should be equal.

You are right...

You are right. But in some cases text files can have significantly different lengths of strings. Use of extra spaces in this case will drastically increase the file size and, hence, the speed of processing.

Acceptable solution...

Proposed solution is obvious, but I would like to find more elegant way.......

Wanted solution...

Thanks, this is a solution I was looking for.

My workaround...

Usually I use plots[display] for this purpose, но это не очень удобно.

Thanks ........

Indeed, this works, but unfortunatelly it is difficult to manage the plots by mouse.

More exactly........

This situation in our languge is called "difficulties of translation" (I don't know, if there is a similar idiom in English). In item 1 I wished to say, that read operator does not read any content from *.m file: neither S1, no something else. But for any other extentsion it is possible to save and read to/from a file without problems. I used the term "returns" in wrong way.

Nevertheless, can you suggest any approach for the problem I described above:

Next, I see, that it is necessary to use output = array for correct saving of solution. But I have to calculate integrals on z(x) and z'(x); besides, the second derivatives of solution (z''(x)) are incorporated into right part of ODE, which I solve later. If z(x) and z'(x) are procedures, I can find z''(x) from initial ODE and fulfill all calculations, But if z(x) and z'(x) are presented by arrays, I should write special procedure for numerical integration of such functions (presented by a set of points). It is desirable to avoid this.

Re:...

Thank you for assistance, I hope, finally I will be able to solve a problem.

What do you mean about puzzling 1 and 2?

1. If S1 saved into external file with m extention (filename.m), it is not possible to read it by read operator. In other words, read returns nothing. Any other extention allows to save and read S1 as needed.

2. Assuming with subs indeed is more simple, now I use this workaround.

Clearly?

Next, I see, that it is necessary to use output = array for correct saving of solution. But I have to calculate integrals on z(x) and z'(x); besides, the second derivatives of solution (z''(x)) are incorporated into right part of ODE, which I solve later. If z(x) and z'(x) are procedures, I can find z''(x) from initial ODE and fulfill all calculations, But if z(x) and z'(x) are presented by arrays, I should write special procedure for numerical integration of such functions (presented by a set of points). It is desirable to avoid this.

BTW, I have the same configuration (M 18.02, W 7 x86 Ult).

Strange behavior of Maple ........

I could not reproduce some your steps. Maybe, the version of Maple or platform matter?

1. When I save the solution to *m file (Maple internal formft), read operator reads nothing from the file. Saving into a file with any other extension don't produce any issues (all works).
2. Syntax subs(S1,z(x)) is ideed simpler, I shall use it for further.
3. There is no problem with printing by showstat(`dsolve/numeric/hermite`) (as you metioned above). But assigning `dsolve/numeric/hermite`:=proc (................... doesn't solve the problem. After entering F1(0.12345)  the following error occurs: Error, (in F1) invalid input: diff received .13, which is not valid for its 2nd argument.

Correct this...

The cause is that the upper limit for variable y is not constant.

Alignment...

Since 2D Output is a character style, I suppose, that alignment is defined by the following:

 1 2 3 4 5 Page 3 of 5
﻿