computing nonlinear weights w1,w2,w3
c1 = 1.0e-1/((eps+s1)^2)
c2 = 6.0e-1/((eps+s2)^2)
```

```
c3 = 3.0e-1/((eps+s3)^2)
w1 = c1/(c1+c2+c3)
w2 = c2/(c1+c2+c3)
w3 = c3/(c1+c2+c3)