fork download 
  1. #include <math.h>
  2. #include <stdio.h>
  3. #include <string.h>
  4.  
  5. #define P  0.47047
  6. #define A  0.3480242
  7. #define B -0.0958798
  8. #define C  0.7478556
  9.  
  10. double errfunc(double z)
  11. {
  12. 	double	erf, t, sign;
  13.  
  14. 	if (z < 0.0) {
  15. 		z = -z;
  16. 		sign = -1.0;
  17. 	} else {
  18. 		sign = 1.0;
  19. 	}
  20. 	t = 1.0 / (1.0 + P * z);
  21. 	erf = 1.0 - (A * t + B * t * t + C * t * t * t) * exp(-(z * z));
  22. 	return erf * sign;
  23. }
  24.  
  25. int main()
  26. {
  27. 	double	x, z, f;
  28. 	char	buf[50];
  29. 	int	y;
  30.  
  31. 	memset(buf, '*', 50);
  32. 	for (x = -3.0; x <= 3.1; x += 0.2) {
  33. 		z = x / pow(M_PI, 0.5);
  34. 		f = 1.0 / 2.0 * (1.0 + errfunc(z));
  35. 		y = f * 50.0;
  36. 		printf("% 3.1f %5.3f %.*s\n", x, f, y, buf);
  37. 	}
  38. 	return 0;
  39. }
  40.  
Success #stdin #stdout 0.02s 1676KB
comments (?)
stdin 
Standard input is empty
stdout 
-3.0 0.008 
-2.8 0.013 
-2.6 0.019 
-2.4 0.028 *
-2.2 0.040 *
-2.0 0.055 **
-1.8 0.075 ***
-1.6 0.101 *****
-1.4 0.132 ******
-1.2 0.169 ********
-1.0 0.212 **********
-0.8 0.262 *************
-0.6 0.316 ***************
-0.4 0.375 ******************
-0.2 0.437 *********************
 0.0 0.500 *************************
 0.2 0.563 ****************************
 0.4 0.625 *******************************
 0.6 0.684 **********************************
 0.8 0.738 ************************************
 1.0 0.788 ***************************************
 1.2 0.831 *****************************************
 1.4 0.868 *******************************************
 1.6 0.899 ********************************************
 1.8 0.925 **********************************************
 2.0 0.945 ***********************************************
 2.2 0.960 ************************************************
 2.4 0.972 ************************************************
 2.6 0.981 *************************************************
 2.8 0.987 *************************************************
 3.0 0.992 *************************************************
https://ideone.com/mm87U
language:
C (gcc 8.3)
created:
12 years ago
visibility:
public

Share or Embed source code

Discover > Sphere Engine API

The brand new service which powers Ideone!

Discover > IDE Widget

Widget for compiling and running the source code in a web browser!