eulerpi$\_$for

Write a function eulerpi$\_$for.m that takes in a positive integer $n$ and returns an approximate value of pi using the following formula:

$\frac{{}^2}{6}$~~$1+\frac{1}{4}+\frac{1}{9}+\frac{1}{16}+$dots$+\frac{1}{n^2}$

If the input value $n$ is not provided, use $n=1000.$

Solve this problem without using a while loop.

Do not use vectorized code. Do not use sum$()$ function

Do not use while loops.