alvgaona

A = [
  -1 -1 1;
  1 3 3;
  -1 -1 5
];

[QQ,RR] = qr(A);

[m, n] = size(A);

alvgaona

x = 2;
 
ax = abs(x);

if (ax < 3.75)
  y = x/3.75;
  y = y*y;
  ans = 1.0+y*(3.5156229+y*(3.0899424+y*(1.2067492 + y*(0.2659732+y*(0.360768e-1+y*0.45813e-2)))));
else
  y=3.75/ax;