260/CH9/EX9.6/9_6.sce


1
2
3
4
5
6
7
8
9
10
11
12
13
14

//Eg-9.6
//pg-400

clc
clear

x=0.6;
p=.3275911;
t=1/(1+p*x);

erfx=1-(.254829592*t-.284496736*t^2+1.421413741*t^3-1.453152027*t^4+1.061405429*t^5)*exp(-x^2);

disp("erf(.6)")
disp(erfx)