Курсовая работа: Моделирование непрерывно-стохастической модели на ЭВМ

fprintf(stream, "%f | %f | %f \n",t,y,z);

t1=t1+h/2;

y1=y+((h/2)*z1);

z1=z1+((h/2)*(e*pow((n0/h),1/2)-b1*z-b2*fabs(z)*z-c1*y-c3*pow(y,3)));

mas1[i]=y1;

i=i+1;

}

/* close the file */

// fclose(stream);

y_rez=0;

y_max=mas1[0];

for (i=0; i<=119;i++)

{

y_rez=y_rez+(pow((mas1[i]-mas[i]),2));

if (mas1[i]>y_max)

y_max=mas1[i];

}

eps=pow(y_rez/n,0.5)/fabs(y_max);

printf("%f epsilon for our operations\n",eps);

fprintf(stream, "%f epsilon for our operations\n",eps);

fclose(stream);

return 0;

}