Здравейте!
Искам да помоля някой по опитен от мен (ниво начинаещ) да ми обясни как се съставя програма, която да намира и извежда корените на кубично уравнение (уравнение, в което най-висшата степен е трета), при въведени коефициенти пред степените. Успях да съставя код на програма, която да прави същото, но с квадратно уравнение. Знам как се решава кубично уравнение и мога да го реша на лист хартия, но изпитвам затруднения при реализацията на кода. Мисля че частта, при която се използва схема на Хорнер може да се изрази с цикъл while, ала незнам как да намеря стойностите на променливата, която да дели коефициентите във схемата на Хорнер (в цикъла). Тази променлива трябва да е делителите на свободния член върху делителите на коефициента пред висшата степен. Принципно знам как да намеря делители на число с цикъл for, обаче този цикъл тук е безполезен, понеже делителите може да са отрицателни числа, тоест да са от реален тип, а доколкото знам с цикъл for може да се инициализират само променливи от целочислен тип. Та въпросът ми е как да намеря тези стойности (делителите на  свободния член върху делителите на коефициента пред висшата степен) и да ги присвоя към променлива, която да дели след това коефициентите в цикъл while. Оттам мисля, че мога да се спрява.
Моля някой да ми помогне.

А това е кода, който съставих за решаване на квадратно уравнение, ако това може да помогне по някакъв начин:

#include <iostream>// reshavane na kvadratno uravnenie
#include <cmath>
using namespace std;
int main ()
{
	double a, b, c, x, D, x1, x2;
	cin>>a>>b>>c;
	a*pow(x,2)+b*x+c==0;
	D=(b*b)-(4*a*c);
	if (a==0) 
	{
	if (b==0) 
	{
	if (c==0) cout<<"Vsqko x e reshenie"<<endl;
	else cout<<"Nqma reshenie"<<endl;
    }
	else 
	{x=(-(c/b));
	cout<<x<<endl;}
}
else 
{
	if (D==0) 
	{x1=x2=(-(b/(2*a))); 
	cout<<x1<<endl;}
	else
	{
		if (D<0) cout<<"Nqma reshenie"<<endl;
		else 
		{x1=(-b+sqrt(D))/(2*a); 
		x2=(-b-sqrt(D))/(2*a); 
		cout<<x1<<"  "<<x2<<endl;}
		
	}
}
system ("pause");
return 0;
}