@Mark44
I tried looking for the DO loop but haven't been able to find what it yet. What I did instead was modify the code this way to somehow control NSTEP:
11 FORMAT('C(1)=',D20.10,' XMIN(1)=',D20.10,' XMIN(5)=',D20.10)
DO 111, NSTEP = 1,5
C Path call for bifurcation points...
I have attached 4 files which contains NSTEP and ISTEP occurrences. It would be great if anyone can have a look and spot anything relevant to help me write an output code to record the required arrays :)
Been struggling a tad too much :(
The file d02paft.txt has only istep occurrences. The...
Aah, I see what I missed there. My bad. However, this is not what I am trying to achieve :(. I had initially coded the program with IMPLICIT A-Z which is what the reference guide asked me to do. However, it was very hard for me to spot the errors using that. Then I tried using IMPLICIT NONE...
If you look at the output here: http://pastebin.com/kuQA22vS
You'll notice that the values of arrays C(1), XMIN(1),XMIN(5),XMIN(9),XMAX(1),XMAX(5),XMAX(9),AMP(1),AMP(5),AMP(9) are different for every ISTEP. I need to record all these 10 values for N number of ISTEP's.
So suppose I run the...
Hi, thank you for your assistance guys.
@Mark44
What I believe is that, this code
DO I = 1, N
WRITE(4, 11) C(I), XMIN(I), XMAX(I)
END
..will give me values of...
C(1), XMIN(1),XMAX(1) as the first line of output.
C(2), XMIN(2), XMAX(2) as the second line of output
C(3), XMIN(3), XMAX(3)...
Dug through the library and found the code below governing XMIN and XMAX. The changes in C are due to multiple calculations from 84 files in the library. I am unsure on what in particular causes the output above as I am following the guide given to me. If there is some way to figure it out by...
@Mark44
I believe that you are trying to find the mention in the code regarding XMIN, XMAX and AMP. There isn't. The code is calling the library to give that output. I searched the library for the mention of xmin, xmax and amp and it points to a file called libpath.a
I have attached the image...
@Mark44
The code below
JOB = 2
JBIF = 0
NSTEP = 20
CALL PATH (FCN,JAC,N,IJAC,X,TPER,NP,INP,TP,TQ,LAMDA,
* JCASE,JOB,JBIF,ILAM,IER,W,IW)
Gives the output (mentioned only the first 3 entries as PF wouldn't allow me to add all 20 entries:
PATH...
Thank you very much for looking into my problem.
Based on the suggestions, I modified the code such that the output is written in a new file called xmin.dat
So, I modified the code as
OPEN(1,file='propeq.dat')
OPEN(2,file='propeq.dat1')
OPEN(3,file='integrate.csv')...
Hi,
I am having a bit of trouble trying to record the output from my code and any help is highly appreciated.
Attached is the code that needs to be looked into (changed extension to .txt from .f as PF wouldn't allow it) and a screenshot of output that I am trying to record.
What I am trying...
Hi there,
I am trying to achieve a particular set of graphs based on duffing equation graphs which you can find here:
It shows a graph which is Amplitude (X(1)) vs. Omega graph. With the code that I have deduced and the formula used, I am able to see that X(1) is varying at every step...
So, I had initially kept the values of X as 0.0D0 for all my 14 variables and it wouldn't populate any results. This was mainly done based on the guide that I was referring to.
On contacting my tutor, he placed the following values instead?
X(1) = 0.0D-01
X(2) = 0.0D-02...
Oh okay.. I thought that was a way to link the library. My bad, sorry about that. I'll figure that you. Thank you very much for all the help you have provided. I highly appreciate it.