for loop counter

for loop counter - Messages

#1 Posted: 10/11/2017 2:17:44 PM
timgow

timgow

0 likes in 27 posts.

Group: User

Hi,
I wish to loop through the YR formula, in pic below, for values of z being 0,0.2,0.4,0.6,0.8,1.0,1.2,1.4,1.589.
My progress so far is to produce an array evenly spaced z values and if the 1.0001 factor is not applied to hc_full then the last element of the z array is 1.362, rather than 1.589.
It is obvious to me that I have some missing skill here, please throw me some hints.
graph.JPG
Regards, Tim
#2 Posted: 10/11/2017 3:25:59 PM
Jean Giraud

Jean Giraud

983 likes in 6866 posts.

Group: User

Just define your 'x' vector of values, include in small loop.

AtWill.sm (6 KiB) downloaded 60 time(s).
#3 Posted: 10/11/2017 3:49:55 PM
timgow

timgow

0 likes in 27 posts.

Group: User

Hi Jean,
Thank you for the example file, I will study this for better understanding.
I have longer arrays to create also so wish to automate; if I make z as in the pic or with equal spacing, is there a means of appending the final (1.589) value to the end of the array?
Tim
#4 Posted: 10/11/2017 5:42:16 PM
CBG

CBG

77 likes in 312 posts.

Group: User

I do not see any problem.

hc_full.png
#5 Posted: 10/11/2017 10:22:09 PM
Jean Giraud

Jean Giraud

983 likes in 6866 posts.

Group: User

Wrote

I have longer arrays to create also so wish to automate;



If you have longer array, equally spaced => as easy.
Naturally, if not equally spaced, must be manual.

AtWill.sm (13 KiB) downloaded 39 time(s).
#6 Posted: 10/11/2017 10:28:28 PM
Jean Giraud

Jean Giraud

983 likes in 6866 posts.

Group: User

... if they are not equally spaced, but if there is a bit
of logic in the spacing, we can select from a list that
can be made/refactored as you wish at will from a little
counter loop.
#7 Posted: 10/12/2017 5:28:21 AM
timgow

timgow

0 likes in 27 posts.

Group: User

Jean,
Thank you for the counter loop example, this will stretch my programming comprehension.

CBG,
Thank you for typing the work from scratch, and showing no issue with the 1.589 value.
I redefined hc_full immediately above the z range call and then the output changed to match your results.

The hc_full is defined with an if statement.
Capture1a.JPG

The 1.589 value is absent from the z range.
Capture1.JPG

The hc_full is redefined directly above the z range call (to override the value produced by the if statement) and the 1.589 value is now present in the z range.
Capture2.JPG

Have I contravened a working principle in some way?

Tim
#8 Posted: 10/12/2017 5:38:09 AM
timgow

timgow

0 likes in 27 posts.

Group: User

Redefining hc_full using hc_full:=hc_full also fixes the z range array.
Tim
#9 Posted: 10/12/2017 1:05:33 PM
Jean Giraud

Jean Giraud

983 likes in 6866 posts.

Group: User

Wrote

Have I contravened a working principle in some way?



Maybe not ? work sheet would help and no help.
Smath 6179 does not take index on the RHS.

hopper.sm (187 KiB) downloaded 34 time(s).
#10 Posted: 10/12/2017 5:45:07 PM
Martin Kraska

Martin Kraska

1222 likes in 2150 posts.

Group: Moderator

for(3)-loops with initial, second and maximum value aren't robust if floating point values are specified. The internal representation of such numbers may differ from what is displayed due to round-off.

Due to this it is not good programming style to test the equality of two floating point numbers without specifying a tolerance. If you want to be safe, you either introduce such a tolerance or you generate the list of values based on a list of integers.

It should be possible to modify for(3) such that the tolerance is generated as a fraction of the difference between first and second value.

2017-10-12 22_36_42-SMath Studio - [Blatt4_].png
loop.sm (4 KiB) downloaded 36 time(s).
Martin Kraska Pre-configured portable distribution of SMath Studio: https://en.smath.info/wiki/SMath%20with%20Plugins.ashx
#11 Posted: 10/12/2017 10:12:01 PM
Jean Giraud

Jean Giraud

983 likes in 6866 posts.

Group: User

True Martin: that damned floating point often creates a very small imaginary.
#12 Posted: 10/12/2017 11:22:25 PM
Jean Giraud

Jean Giraud

983 likes in 6866 posts.

Group: User

In the context of this topic, there was an example in Mathcad 8.
About solving an RK DE relative to "Orbit" by Steven Finch.
  • New Posts New Posts
  • No New Posts No New Posts