Not a member of Pastebin yet?
Sign Up,
it unlocks many cool features!
- program Biqrt;
- {$APPTYPE CONSOLE}
- uses
- SysUtils,Math;
- var
- a,b,c,D,x,y,f1,f2,modul1, modul2,z1,z2: double;
- begin
- Readln(a,b,c);
- Writeln(a:2:2,'x^4+',b:2:2,'x^2+',c:2:2,'=0');
- D:= sqr(b) - 4*a*c;
- if D >= 0 then
- begin
- y:= 0;
- x:= (-b + sqrt(D))/(2*a);
- f1:= (pi*(0.5*sign(x)-0.5));
- modul2:= sqrt(sqrt(sqr(x)+sqr(y)));
- x:= (-b - sqrt(D))/(2*a);
- f2:= (pi*(0.5*sign(x)-0.5));
- modul1:= sqrt(sqrt(sqr(x)+sqr(y)));
- end;
- if D < 0 then
- begin
- x:= -b / (2*a);
- y:= sqrt(-D)/(2*a);
- f1:= sign(y) * (pi*(0.5*(sign(x))-0.5) + arctan(y/x));
- f2:= sign(-y) * (pi*(0.5*(sign(x))-0.5) + arctan(-y/x));
- modul1:= sqrt(sqrt(sqr(x)+sqr(y)));
- modul2 := modul1;
- end;
- z1:= modul1*cos((f2+2*pi*0)/2); z2:= modul1*sin((f2+2*pi*0)/2);
- writeln(z1:2:2,'+(',z2:2:2,'i)');
- z1:= modul1*cos((f2+2*pi*1)/2); z2:= modul1*sin((f2+2*pi*1)/2);
- writeln(z1:2:2,'+(',z2:2:2,'i)');
- z1:= modul2*cos((f1+2*pi*0)/2); z2:= modul2*sin((f1+2*pi*0)/2);
- writeln(z1:2:2,'+(',z2:2:2,'i)');
- z1:= modul2*cos((f1+2*pi*1)/2); z2:= modul2*sin((f1+2*pi*1)/2);
- writeln(z1:2:2,'+(',z2:2:2,'i)');
- readln;
- end.
Advertisement
Add Comment
Please, Sign In to add comment