Доброго времени суток, уважаемые эксперты!
Требуется реализовать метод простой итерации для решения нелинейной системы состоящей из семи уравнений. Вот такой:

x[0]=x[1]^3+x[2]-x[3]+2x[5]+x[6]-13
x[1]=10x[0]-x[2]+x[3]+100x[4]^2-x[5]-x[6]+6
x[2]=x[0]^2+x[1]+x[3]+x[4]+x[5]+x[6]-2
x[3]=x[0]-x[1]-x[2]+x[4]^2+x[5]-x[6]+2
x[4]=x[0]^2+2x[1]-2x[2]-x[3]+x[5]+x[6]-5
x[5]=4x[0]^2+2x[1]+2x[2]-4x[3]+5x[4]-x[6]-8
x[6]=5x[0]+5x[1]-3x[2]+4x[3]-x[4]^2+x[5]-3

корни системы уравнений:
x[0]=0
x[1]=2
x[2]=1
x[3]=-1
x[4]=0
x[5]=1
x[6]=1

И есть еще другие корни...

Начальное приближение задаётся таким образом:
x[0]:=-0.75+(random()/2);
x[1]:=1.75+(random()/2);
x[2]:=0.75+(random()/2);
x[3]:=-1.25+(random()/2);
x[4]:=-0.25+(random()/2);
x[5]:=0.75+(random()/2);
x[6]:=0.75+(random()/2);

Я пробовал реализовать метод, но он у меня почему-то не сходится... Кто может, напишите алгоритм пожалуйста в какой-нибудь из этих сред: дельфи, maple
алгоритм нужен сегодня.