Obliczanie pierwiastków metodą Newtona

Witam, nie mogę dojść do tego czemu to nie działa jak powinno :/

#include <cstdio>
#include <cstdlib>
#include <cmath>
#include <iostream>
using namespace std;

int main() {
double s, a, b;
cin >> s;
a = 0;
b = 1;
  while(abs(a-b) > 0.0001) {
   a = b;
	b = 0,5*(s/a+a);
	
	}
	cout << b;
	}

cokolwiek wpiszę to wyskakuje 0

do
{
   a = b;
   b = a - (a*a - s) / 2*b;            //imho
}while( abs(a - b) > 0.0001 );

Rany, wszystko przez przecinek, niepotrzebne zawracanie głowy, temat do skasowania, dobranoc ;)