Табулирование функции, полярные координаты

[Su-34]

Табулирование функции, полярные координаты

Помогите пжлста. Мне надо затабулировать функцию р=sin(f/2), заданную полярными координатами с помощью Visual Fortran.
Что я не так сделал, что нужно исправить чтобы в екселе получился график 2). Пока он мне рисует только правую часть 1)

program lab1
real a,b,h,x,y
integer i,N
open(10,file='lab1.txt')
a=-7
b=1
N=31
h=(b-a)/(N-1)
do i=1,N
f=a+h*i
r=sin(f/2)
x=r*cos(f)
y=r*sin(f)
print*,x,y
write(10,*)x,y
end do
end

 

[Vladimir_S]

a=-7 b=1

А почему Вы берёте такие странные пределы? Ведь справа указано: от 0 до 4π, т.е от 0 до где-то 12.57. Такие и надо брать.
Переписал Вашу программу с VF на Паскаль, практически ничего в ней не меняя:

var
 a,b,h,x,y,r,f:real;
 i,N:integer;
 ff:TEXT;
begin
 Assign(ff,'D:\lab1.dat');
 ReWrite(ff);
 a:=0;
 b:=Pi*4;
 N:=100;
 h:=(b-a)/N;
 for i:=0 to N do
  begin
   f:=a+h*i;
   r:=sin(f/2);
   x:=r*cos(f);
   y:=r*sin(f);
   writeln(ff,x:8:5,'  ',y:8:5);
  end;
 Close(ff);
end.

запустил и - пожалуйста (построено в пакете ORIGIN):

Так что всё дело только в правильной установке пределов.

 

[Su-34]

Можно еще вопрос? То же самое надо сделать с р=(phi)^2. Ведь должна получиться парабола? Тогда надо писать:
program lab1
real a,b,h,x,y
integer i,N
open(10,file='lab1.txt')
a=-1
b=1
N=31
h=(b-a)/(N-1)
do i=1,N
f=a+h*i
r=f**2
x=r*cos(f)
y=r*sin(f)
print*,x,y
write(10,*)x,y
end do
end

Что не так?:tehnari_ru_117:

 

[Vladimir_S]

Ведь должна получиться парабола?

Это еще в честь чего? В малой окрестности нуля должна быть функция вида у=х**(3/2). А у Вас что получается?