Reply to thread

Message

<blockquote data-quote="noze2000" data-source="post: 45538">I wrote this small program for C on a windows machine. It comlies and runs perfect.yet it will not compile in Xcode for some reason. There are a few mistakes but most of them comes from xcode not recognizing #include &lt;fstream&gt;. Can anyone help me clean up the code a little since I cannot find any documentation on it on the web<img src="data:image/gif;base64,R0lGODlhAQABAIAAAAAAAP///yH5BAEAAAAALAAAAAABAAEAAAIBRAA7" class="smilie smilie--sprite smilie--sprite3" alt=":(" title="Frown&nbsp; &nbsp; :(" loading="lazy" data-shortname=":(" /> ?regardsnoze2000[code]/*&nbsp; &nbsp; bits of the code has been adapted from&nbsp; &nbsp; [url]http://www.cs.caltech.edu/~arvo/code/Rand.C[/url]*//*-------------------------------------------------------------------------*&nbsp;* M E T H O D 2&nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp;&nbsp; *&nbsp;*&nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp;&nbsp; *&nbsp;* From &quot;The Multiple Prime Random Number Generator,&quot; by Alexander Haas,&nbsp;&nbsp; *&nbsp;* ACM Transactions on Mathematical Software, Vol. 13, No. 4, December&nbsp; &nbsp;&nbsp; *&nbsp;* 1987, pp. 368-381.&nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; *&nbsp;*&nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp;&nbsp; *&nbsp;*-------------------------------------------------------------------------*/#include &lt;fstream&gt;#include &lt;stdlib.h&gt;#include &lt;stdio.h&gt;#include &lt;math.h&gt;#include &lt;time.h&gt;#define pi 3.14159265358979&nbsp; &nbsp; #ifndef ABS#define ABS( x ) ((x) &gt; 0 ? (x) : -(x))#endifint time_steps_per_periode;double cvar,perioder,tid,gnst, s0, rr, X, dt, tilf_tal, tilf_tal_1, tilf_tal_2,stdaf,geometrisk_pris, vol,Pris_estimat,var_a ,tvar_a,var_g,tvar_g ,covarians, varians_estimat, pris_a,pris_g,pris_1=0.0,pris_2=0.0;int m, i, j ,r , brugt_normale_tal, nye_tal,u,u0,u1,u2,u3,u1000;clock_t tid1,tid2,temptid;void Seed(){&nbsp; &nbsp; &nbsp; &nbsp; m=971,i=11113,j=104322;&nbsp; &nbsp; r=4181;}void initialize() {&nbsp; &nbsp; &nbsp; &nbsp; tid1=clock();&nbsp; &nbsp; Seed();&nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; //&nbsp; &nbsp; sÊtter begyndelse for tilfÊldighedrÊkke&nbsp; &nbsp; }&nbsp; &nbsp; double random()&nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; // lav eet tilfÊldigt tal&nbsp; &nbsp; {&nbsp; &nbsp; &nbsp; &nbsp; if( (m += 7&nbsp; &nbsp; ) &gt;=&nbsp;&nbsp; 9973 ) m -=&nbsp; 9871;&nbsp; &nbsp; if( (i += 1907 ) &gt;=&nbsp; 99991 ) i -= 89989;&nbsp; &nbsp; if( (j += 73939) &gt;= 224729 ) j -= 96233;&nbsp; &nbsp; r = ((r * m + i + j) % 100000) / 10;&nbsp; &nbsp; return r * 1.00010001E-4;&nbsp; &nbsp; }void polar()&nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; // brug Polar metoden til at lave 2 normalfordelte tal{&nbsp; &nbsp; double v1,v2,u1,u2,n=0.0;&nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; // sÂ der laves 2 tal til at starte med&nbsp; &nbsp;&nbsp; do {&nbsp; &nbsp; &nbsp; &nbsp; u1=random();&nbsp; &nbsp; &nbsp; &nbsp; u2=random();&nbsp; &nbsp; &nbsp; &nbsp; v1=2*u1-1;&nbsp; &nbsp; &nbsp; &nbsp; v2=2*u2-1;&nbsp; &nbsp; &nbsp; &nbsp; n=v1*v1+v2*v2;&nbsp; &nbsp; }while (n&gt;1);&nbsp; &nbsp; &nbsp; &nbsp; tilf_tal_1 = v1 * sqrt(-2 * log(n) / n);&nbsp; &nbsp; &nbsp; &nbsp; tilf_tal_2 = v2 * sqrt(-2 * log(n) / n);}void taeller_af_tilf_tal ()&nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; //&nbsp; &nbsp; kontrollerer at de tilfÊldige tal ikke benyttes flere gange{&nbsp; &nbsp; switch(nye_tal){&nbsp; &nbsp; case 0 : {polar(); tilf_tal = tilf_tal_1 ;}&nbsp; &nbsp; &nbsp; &nbsp; break;&nbsp; &nbsp; case 1 : {tilf_tal = tilf_tal_2 ; nye_tal=0;}&nbsp; &nbsp; &nbsp; &nbsp; break;&nbsp; &nbsp; }}&nbsp; &nbsp; void tid_tilbage(int x,int y)&nbsp; &nbsp; {&nbsp; &nbsp; &nbsp; &nbsp; int u1000=(int)y/500,u0=(int)y/100,u1=(int)y/3,&nbsp; &nbsp; //&nbsp; &nbsp; variable til tidstagning&nbsp; &nbsp; u2=(int)y/2, u3=(int)(2.0/3.0*y);// rutine til estimering af tid tilbage for udregningen&nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; double a,b,c;;&nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; if(x==u1000){&nbsp; &nbsp; temptid = clock () - tid1;&nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; a=(500.0*(temptid)/((double)CLOCKS_PER_SEC));&nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; b=fmod(a , 60); // sekunder tilbage&nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; c=int(a/60);&nbsp; &nbsp; // minutter tilbage&nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; printf(&quot; 0.2%% : Estimeret tid tilbage: %2.0f min %2.0f sek\n\n&quot;,c ,b);}&nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; if(x==u0) {&nbsp; &nbsp; temptid=clock ();&nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; a=(100.0*(temptid)/((double)CLOCKS_PER_SEC));&nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; b=fmod(a , 60); // sekunder tilbage&nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; c=int(a/60);&nbsp; &nbsp; // minutter tilbage&nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; printf(&quot; 1%% : Estimeret tid tilbage: %2.0f min %2.0f sek\n\n&quot;,c ,b);}&nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; if(x==u1) {temptid=clock ();;&nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; a=(2*(temptid)/((double)CLOCKS_PER_SEC));&nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; b=fmod(a , 60); // sekunder tilbage&nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; c=int(a/60);&nbsp; &nbsp; // minutter tilbage&nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; printf(&quot;33%% : Estimeret tid tilbage: %2.0f min %2.0f sek\n\n&quot;,c ,b);}&nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; if(x==u2) {temptid=clock ();&nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; a=((temptid)/((double)CLOCKS_PER_SEC));&nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; b=fmod(a , 60); // sekunder tilbage&nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; c=int(a/60);&nbsp; &nbsp; // minutter tilbage&nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; printf(&quot;50%% : Estimeret tid tilbage: %2.0f min %2.0f sek\n\n&quot;,c ,b);}&nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; if(x==u3) {temptid=clock ();&nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; a=((temptid)/(2.0*(double)CLOCKS_PER_SEC));&nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; b=fmod(a , 60); // sekunder tilbage&nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; c=int(a/60);&nbsp; &nbsp; // minutter tilbage&nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; printf(&quot;66%% : Estimeret tid tilbage: %2.0f min %2.0f sek\n\n&quot;,c ,b);}&nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; }&nbsp; &nbsp; void gem(char *file, double a){&nbsp; &nbsp; using namespace std;&nbsp; &nbsp; ofstream savefile( file, ios::binary | ios::app );&nbsp; &nbsp; savefile &lt;&lt; a &lt;&lt; &quot;\n&quot;;&nbsp; &nbsp; savefile.close();}void gem(char *file, long double a,long double b){&nbsp; using namespace std;&nbsp; &nbsp; &nbsp; ofstream savefile( file, ios::binary | ios::app );&nbsp; savefile &lt;&lt; a &lt;&lt; &quot;\t&quot; &lt;&lt; b &lt;&lt; &quot;\n&quot;;&nbsp; savefile.close();}void init_gem(char *file){&nbsp; &nbsp; using namespace std;&nbsp; &nbsp; ofstream savefile(file,ios::trunc);&nbsp; &nbsp; savefile.close();}&nbsp; &nbsp; void main() {&nbsp; &nbsp; printf(&quot;Monte Carlo simulation med Control variable\n\n&quot;);&nbsp; &nbsp; initialize();&nbsp; &nbsp; &nbsp; &nbsp; double b,c;&nbsp; &nbsp; int files = 100 ;&nbsp; &nbsp; const int tal=100000;&nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; double e[tal];&nbsp; &nbsp; char a[20];&nbsp; &nbsp; for(int w = 0;w&lt;files;w++)&nbsp; &nbsp; {&nbsp; &nbsp; &nbsp; &nbsp; if (w&lt;10) sprintf(a, &quot;normal_%d%d.txt&quot;,0, w+1);&nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; else sprintf(a, &quot;normal_%d.txt&quot;, w+1);&nbsp; &nbsp; &nbsp; &nbsp; char *fil = a;&nbsp; &nbsp; &nbsp; &nbsp; init_gem(fil);&nbsp; &nbsp; &nbsp; &nbsp; printf(&quot;laver 100000 tal...\t&quot;);&nbsp; &nbsp; &nbsp; &nbsp; for(int s=0;s&lt;tal;s++) // lav tal&nbsp; &nbsp; &nbsp; &nbsp; {&nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; taeller_af_tilf_tal();&nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; e[s] = tilf_tal;&nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; }&nbsp; &nbsp; &nbsp; &nbsp; printf(&quot;laver en fil af tallene....\n&quot;);&nbsp; &nbsp; &nbsp; &nbsp; printf(&quot;fil #%d&quot;,w+1);&nbsp; &nbsp; &nbsp; &nbsp; using namespace std;&nbsp; &nbsp; &nbsp; &nbsp; ofstream savefile( fil, ios::binary | ios::app );&nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; for(s=0;s&lt;tal;s++)&nbsp; &nbsp; &nbsp; &nbsp; {savefile &lt;&lt; e[s] &lt;&lt; &quot;\n&quot;;}&nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; savefile.close();&nbsp; &nbsp; }&nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; tid2 = clock();double tid = (tid2 - tid1)/(CLOCKS_PER_SEC / (double) 1000.0);&nbsp; &nbsp; b=fmod(tid , 60); // sekunder tilbage&nbsp; &nbsp; c=int(tid/60);&nbsp; &nbsp; // minutter tilbage&nbsp; &nbsp; printf(&quot;\ntid brugt: %2.0f min %2.5f sek\n\n&quot;,c ,b);&nbsp; &nbsp; printf(&quot;\n%f&quot;,tid);&nbsp; &nbsp; }[/code]</blockquote>

[QUOTE="noze2000, post: 45538"] I wrote this small program for C on a windows machine. It comlies and runs perfect. yet it will not compile in Xcode for some reason. There are a few mistakes but most of them comes from xcode not recognizing #include <fstream>. Can anyone help me clean up the code a little since I cannot find any documentation on it on the web:( ? regards noze2000 [code]/* bits of the code has been adapted from [url]http://www.cs.caltech.edu/~arvo/code/Rand.C[/url] */ /*-------------------------------------------------------------------------* * M E T H O D 2 * * * * From "The Multiple Prime Random Number Generator," by Alexander Haas, * * ACM Transactions on Mathematical Software, Vol. 13, No. 4, December * * 1987, pp. 368-381. * * * *-------------------------------------------------------------------------*/ #include <fstream> #include <stdlib.h> #include <stdio.h> #include <math.h> #include <time.h> #define pi 3.14159265358979 #ifndef ABS #define ABS( x ) ((x) > 0 ? (x) : -(x)) #endif int time_steps_per_periode; double cvar,perioder,tid,gnst, s0, rr, X, dt, tilf_tal, tilf_tal_1, tilf_tal_2,stdaf,geometrisk_pris, vol,Pris_estimat,var_a ,tvar_a,var_g,tvar_g ,covarians, varians_estimat, pris_a,pris_g,pris_1=0.0,pris_2=0.0; int m, i, j ,r , brugt_normale_tal, nye_tal,u,u0,u1,u2,u3,u1000; clock_t tid1,tid2,temptid; void Seed(){ m=971,i=11113,j=104322; r=4181; } void initialize() { tid1=clock(); Seed(); // sÊtter begyndelse for tilfÊldighedrÊkke } double random() // lav eet tilfÊldigt tal { if( (m += 7 ) >= 9973 ) m -= 9871; if( (i += 1907 ) >= 99991 ) i -= 89989; if( (j += 73939) >= 224729 ) j -= 96233; r = ((r * m + i + j) % 100000) / 10; return r * 1.00010001E-4; } void polar() // brug Polar metoden til at lave 2 normalfordelte tal { double v1,v2,u1,u2,n=0.0; // sÂ der laves 2 tal til at starte med do { u1=random(); u2=random(); v1=2*u1-1; v2=2*u2-1; n=v1*v1+v2*v2; }while (n>1); tilf_tal_1 = v1 * sqrt(-2 * log(n) / n); tilf_tal_2 = v2 * sqrt(-2 * log(n) / n); } void taeller_af_tilf_tal () // kontrollerer at de tilfÊldige tal ikke benyttes flere gange { switch(nye_tal){ case 0 : {polar(); tilf_tal = tilf_tal_1 ;} break; case 1 : {tilf_tal = tilf_tal_2 ; nye_tal=0;} break; } } void tid_tilbage(int x,int y) { int u1000=(int)y/500,u0=(int)y/100,u1=(int)y/3, // variable til tidstagning u2=(int)y/2, u3=(int)(2.0/3.0*y);// rutine til estimering af tid tilbage for udregningen double a,b,c;; if(x==u1000){ temptid = clock () - tid1; a=(500.0*(temptid)/((double)CLOCKS_PER_SEC)); b=fmod(a , 60); // sekunder tilbage c=int(a/60); // minutter tilbage printf(" 0.2%% : Estimeret tid tilbage: %2.0f min %2.0f sek\n\n",c ,b);} if(x==u0) { temptid=clock (); a=(100.0*(temptid)/((double)CLOCKS_PER_SEC)); b=fmod(a , 60); // sekunder tilbage c=int(a/60); // minutter tilbage printf(" 1%% : Estimeret tid tilbage: %2.0f min %2.0f sek\n\n",c ,b);} if(x==u1) {temptid=clock ();; a=(2*(temptid)/((double)CLOCKS_PER_SEC)); b=fmod(a , 60); // sekunder tilbage c=int(a/60); // minutter tilbage printf("33%% : Estimeret tid tilbage: %2.0f min %2.0f sek\n\n",c ,b);} if(x==u2) {temptid=clock (); a=((temptid)/((double)CLOCKS_PER_SEC)); b=fmod(a , 60); // sekunder tilbage c=int(a/60); // minutter tilbage printf("50%% : Estimeret tid tilbage: %2.0f min %2.0f sek\n\n",c ,b);} if(x==u3) {temptid=clock (); a=((temptid)/(2.0*(double)CLOCKS_PER_SEC)); b=fmod(a , 60); // sekunder tilbage c=int(a/60); // minutter tilbage printf("66%% : Estimeret tid tilbage: %2.0f min %2.0f sek\n\n",c ,b);} } void gem(char *file, double a) { using namespace std; ofstream savefile( file, ios::binary | ios::app ); savefile << a << "\n"; savefile.close(); } void gem(char *file, long double a,long double b) { using namespace std; ofstream savefile( file, ios::binary | ios::app ); savefile << a << "\t" << b << "\n"; savefile.close(); } void init_gem(char *file){ using namespace std; ofstream savefile(file,ios::trunc); savefile.close(); } void main() { printf("Monte Carlo simulation med Control variable\n\n"); initialize(); double b,c; int files = 100 ; const int tal=100000; double e[tal]; char a[20]; for(int w = 0;w<files;w++) { if (w<10) sprintf(a, "normal_%d%d.txt",0, w+1); else sprintf(a, "normal_%d.txt", w+1); char *fil = a; init_gem(fil); printf("laver 100000 tal...\t"); for(int s=0;s<tal;s++) // lav tal { taeller_af_tilf_tal(); e[s] = tilf_tal; } printf("laver en fil af tallene....\n"); printf("fil #%d",w+1); using namespace std; ofstream savefile( fil, ios::binary | ios::app ); for(s=0;s<tal;s++) {savefile << e[s] << "\n";} savefile.close(); } tid2 = clock(); double tid = (tid2 - tid1)/(CLOCKS_PER_SEC / (double) 1000.0); b=fmod(tid , 60); // sekunder tilbage c=int(tid/60); // minutter tilbage printf("\ntid brugt: %2.0f min %2.5f sek\n\n",c ,b); printf("\n%f",tid); } [/code] [/QUOTE]

Verification

What is Apple's smallest desktop computer called?