Написание программы для вычисления числа пи

N=1000;
n=1:N;
pi=sqrt(...);

function out=thedoctor(pi)
N=1000;n=1:N;
S=0;for T=1\((2.*n-1).^2.*(2.*n+1).^2);S,S=S+T,end
S=(pi^2-8)/16
end

N=1000;
n=1:N;
PI=sqrt(16*sum(1./(4*n.^2-1).^4)+8)

function S=summa(N)
n=1:N;
S=0;for T=1\((2.*n-1).^2.*(2.*n+1).^2);S=S+T;end 
PI=sqrt((16/-8)*S)
 
end

function PHI=summa(N)
s =0; %первоначальная инициализация
for n =1:N; %рассчитываем правую часть
s = s +(1)/(((2*n - 1)^2)*((2*n +1)^2)) %избыток скобок для наглядности
end; %закрываем цикл
s = s*16 % избавляемся от знаминателя
s = s +8 %в переменной s у нас пиквадрат
s = sqrt(s) %получаем число пи.
PHI = s; % присваиваем входному аргументу
end %конец.

N=1000;
n=1:N;
PI=sqrt(16*sum(1./(4*n.^2-1).^2)+8)