One Hat Cyber Team

Your IP : 216.73.216.148

Server IP : 194.44.31.54

Server : Linux zen.imath.kiev.ua 4.18.0-553.77.1.el8_10.x86_64 #1 SMP Fri Oct 3 14:30:23 UTC 2025 x86_64

Server Software : Apache/2.4.37 (Rocky Linux) OpenSSL/1.1.1k

PHP Version : 5.6.40

Buat File | Buat Folder

Dir : ~/home/vikul/

Edit File: diag_cont.cpp

#include <stdlib.h>

#include <iostream.h>

#include <fstream.h>

#include <string.h>

#include <math.h>

#include <stdio.h>

main ()

{

double x, y, z, K, Diffr;

double a=0.2, b=0.2, c=5.7, Tau=11.758626;

double Time, TimePrev;

double x0=4.85320967,  y0=3.83371588,  z0=0.39558092, x0l, y0l, z0l;

double StepInt=0.00001, StepK=0.001;

long NumPerSkip=100, NumPerInt=500;

long LongNumSkip=100, LongNumInt=200  ;

double *ArrXIni, *ArrYIni, *ArrZIni;

double *ArrX, *ArrY, *ArrZ, *ArrX1, *ArrY1, *ArrZ1,  *ArrTime,*TmpAdr;

double *ArrR1,*ArrR2,*ArrR3;

double *ArrL1,*ArrL2,*ArrL3;

unsigned long NumStepPer;

long i,j;

char Str[1000], fname[255], StrNum[10];

FILE *FileOut, *FileTrjR, *FileOutL, *FileTrjL, *FileOutTimeL, *FileProd ;

NumStepPer = (unsigned long)(Tau / StepInt)+1;

/*   cout<<"NumStepPer=" <<NumStepPer<<endl; 
*/
  ArrXIni =(double*)malloc((NumStepPer+2)*sizeof(double));

ArrYIni =(double*)malloc((NumStepPer+2)*sizeof(double));

ArrZIni =(double*)malloc((NumStepPer+2)*sizeof(double));

if(!ArrXIni || !ArrYIni || !ArrZIni) {

cout << "Cannot allocate memory for StateIni arrays. Exit. \n";   exit(1);

}

for(i = 0; i <= NumPerSkip; i++) {

if( i == NumPerSkip) {ArrXIni[0] = x0; ArrYIni[0] = y0; ArrZIni[0] = z0;}

for(j = 1; j <= NumStepPer; j++) {

x=x0+StepInt*(-y0-z0);

y=y0+StepInt*(x0+a*y0);

z=z0+StepInt*(b+z0*(x0-c));

x0=x;  y0=y;  z0=z;

if( i == NumPerSkip) {ArrXIni[j] = x0; ArrYIni[j] = y0; ArrZIni[j] = z0;}

}

ArrX =(double*)malloc((NumStepPer+2)*sizeof(double));

ArrY =(double*)malloc((NumStepPer+2)*sizeof(double));

ArrZ =(double*)malloc((NumStepPer+2)*sizeof(double));

if(!ArrX || !ArrY || !ArrZ || !ArrTime) {

cout << "Cannot allocate memory for State arrays. Exit. \n";   exit(1);

}

FileOut=fopen("d_c_2_0.001.dat", "wt");

if(!FileOut) {

cout <<"Cann't open file.Exit"<<endl<<flush;exit(1);

}

for (K=0;K<=0.3;K=K+StepK){

/*  cout<<"K="<<K<<endl;*/

Diffr=0;

for(i=0;i<=NumPerInt;i++){

ArrX[0]=ArrXIni[NumStepPer]; ArrY[0]=ArrYIni[NumStepPer]; ArrZ[0]=ArrZIni[NumStepPer];

for(j = 1; j <= NumStepPer; j++) {

x0 = ArrX[j-1]; y0 = ArrY[j-1]; z0 = ArrZ[j-1];

ArrX[j]= x0+StepInt*((-y0-z0)+K*(ArrXIni[j-1]-x0));

ArrY[j]=y0+StepInt*((x0+a*y0)+K*(ArrYIni[j-1] - y0));

ArrZ[j]=z0+StepInt*((b+z0*(x0-c))+K*(ArrZIni[j-1]-z0));

if( i == NumPerInt) Diffr += fabs(ArrXIni[j-1] - x0)+fabs(ArrYIni[j-1] - y0)+fabs(ArrZIni[j-1] - z0);

}

TmpAdr = ArrX; ArrX = ArrXIni; ArrXIni = TmpAdr;

TmpAdr = ArrY; ArrY = ArrYIni; ArrYIni = TmpAdr;

TmpAdr = ArrZ; ArrZ = ArrZIni; ArrZIni = TmpAdr;

}

fprintf(FileOut,"%15.8lf  %15.8lf\n", K , Diffr/NumStepPer);

fflush(FileOut);

}

fclose( FileOut);

/*  fclose(FileOutR);

fclose( FileTrjL);

fclose(FileOutL);

fclose(FileOutTimeL);

fclose(FileProd); */

free(ArrXIni);

free(ArrYIni);

free(ArrZIni);

free(ArrX);

free(ArrY);

free(ArrZ);

/*  free(ArrR1);

free(ArrR2);

free(ArrR3);

free(ArrL1);

free(ArrL2);

free(ArrL3);           */

}