Итерация С Фиксированной Точкой. Поиск корней.

---

Я пытаюсь написать программу для поиска корней с помощью метода итерации с фиксированной точкой, и каждый раз, когда я запускаю ее, я получаю ноль. ввод p0=1, Tol=.01

Не мог бы кто-нибудь помочь? Я попытался следовать алгоритму в книге, но я все еще новичок в программировании и не очень хорошо читаю их. Заранее спасибо.


Алгоритм:

Входное начальное приближение p0; допуск TOL; максимальное число итераций N0.
Вывод приближенного решения P или сообщения об аварии.

Шаг 1 Установите i=1.
Шаг 2, пока i <= N0, выполните шаги 3-6.
Шаг 3 Установите p=g(p0). (Вычислить число Пи.)
Шаг 4 Если |p-p0|<tol mode="hold"> OUTPUT (p); (процедура прошла успешно.)
ОСТАНОВИТЬ.
Шаг 5 Установите i=i+1.
Шаг 6 Установите p0=p. (Обновление p0.)
Шаг 7 вывод ('метод не удался после N0 итераций, N0=', N0);
(Процедура оказалась неудачной.)
ОСТАНОВИТЬ.

#include <stdio.h>
#include <math.h>

double f(double x)
{
return x*x*x*x-3*x*x-3;  //change equation for each problem
}

int main()
{
    double p, p0, Tol;
    int i=1;
    int No;

    printf("Enter approximate p: ");
    scanf ("%lf", &p0);

    printf("Desired Tolerance: ");
    scanf ("%lf", &Tol);

    printf("Maximum Iterations: ");
    scanf ("%d", &No);


    while (i<=No)
    {
        p = p0;

        if((fabs(p-p0))<Tol)
        {
            printf("%lf", &p);
            break;
        }
        printf("Iteration %d: Current value = %lf\n", i, p);

        i++;  //i=i+1
        p0=p;

        if (i>No)
        {
        printf("Method Failed after %d", No);
        printf(" iterations");
        }

Вот обновленный вариант:

#include <stdio.h>
#include <math.h>

double f(double x)
{
return x*x*x*x-3*x*x-3;  //change equation for each problem
}

double g(double x)
{
return pow(3*x*x+3,.25);
}

int main()
{
    double p, p0, Tol;
    int i=1;
    int No;

    printf("Enter approximate p: ");
    scanf ("%lf", &p0);

    printf("Desired Tolerance: ");
    scanf ("%lf", &Tol);

    printf("Maximum Iterations: ");
    scanf ("%d", &No);


    while (i<=No)
    {
        p = g(p0);

        if((fabs(p-p0))<Tol)
        {
            //printf("%lf", &p);
            break;
        }
        printf("Iteration %d: Current value = %lf\n", i, p);

        i++;  //i=i+1
        p0=p;

        if (i>No)
        {
        printf("Method Failed after %d", No);
        printf(" iterations");
        }

    }

}

---