Вопрос по Compaq Visual Fortran

[Fencer]

Как задать такую разрывную функцию? Я дал параметр c,d. Но график не рисует...

program lab1
real a,b,c,d,h,x,y
integer i,N
open(10,file='lab1.txt')
a=-5
b=-13./10.
c=23./10.
d=5
N=31
h=(b-a)/(N-1)
do i=1,N
x=a+(i-1)*h
y=2-(3./4.)*(sqrt(x**2-x-3.))
print*,x,y
write(10,*)x,y
end do
end

[Selyd]

У тебя под корнем функция, которая на интервале [-1;+2] (приблизительно,
посчитай точно) не существует. Поэтому надо дать два цикла -
один цикл для левой ветки, а второй - для правой ветки кривой.
Дырка, кажется, получится в автомате. Удачи.

[Fencer]

как прописать правильно два цикла?

[Selyd]

как прописать правильно два цикла?

Найди корни подкоренного выражения.
Левая ветвь
Меньший корень - правый конец для левого интервала. Левый конец левого интервала выбирай сам.
Левая ветвь на отрезке [Левый конец; Меньший корень]
Назначь N для левой ветви и - цикл.
Правая ветвь
Больший корень - левый конец интервала для правого интервала.
Правый конец правого интервала выбирай сам.
Назначь N для правой ветви и - цикл.

Назначь N для правой ветви и - свой цикл.

[Fencer]

Как назначить N

[Selyd]

Если у тебя отрезок [a;b] и шаг h, тогда N=(b-a)/h.
Как ты делал сам, только теперь два раза. Это ты выбираешь, назначаешь.
Выбирая левую и правую границы, и шаг будеш менять детализацию рисунка.
В разумных пределах конечно, чтобы влезло в поле екрана и было наглядно.
Дыра от меньшего до большего корня подкоренного выражения
стабильна. Тут ты ничего сделать не можешь.
Посмотрел твой рисунок ещё раз. Ну там же всё правильно!
Есть левая ветвь, обрыв в -1. Есть правая, обрыв в +2.
И дыра на месте.

[Fencer]

Какая-то ошибка в sqrt...что там не так?

[Selyd]

Ты двигаешься слева направо. Аргумент
становится больше меньшего корня и ты
под корнем получаешь отрицательное число.
Речь и идёт о том, чтобы пропустить отрезок
[Xmin; Xmax] за счёт двух отдельных циклов
для левой и правой ветвей графика.

[Fencer]

При таком раскладе в екселе рисует только правую часть графика..!

program lab1
real a,b,h,x,y
integer i,N
open(10,file='lab1.txt')
a=-5
b=-7./5.
N=31
h=(b-a)/(N-1)
do i=1,N
x=a+(i-1)*h
y=2-(3./4.)*(sqrt(x**2-x-3))
print*,x,y
END DO
a=12./5.
b=5
N=31
h=(b-a)/(N-1)
do i=1,N
x=a+(i-1)*h
y=2-(3./4.)*(sqrt(x**2-x-3.))
print*,x,y
write(10,*)x,y
end do
end

Что делать??

[wital1984]

При таком раскладе в екселе рисует только правую часть графика..!
....

Честно говоря, в Visual Fortran я не в зуб ногой, но почему у вас нет в первом цикле write(10,*)x,y ?, может, в этом дело.
и почему в первом цикле нет точки после 3 при нахождении y: y=2-(3./4.)*(sqrt(x**2-x-3))
А так вроде все логично верно, если рисует правую часть, то и левую должен рисовать

[ScoRpiON]

всё сделал.

[ScoRpiON]

Вообщем вопрос в следующем. Есть функция y=tan(x/2-1). Период беру от -1.13 до 5.13. Но в екселе рисует так:
Как сделать чтобы была плавная линия, взять ограничения по y? Помогите пожалуйста

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

[Selyd]

Возьми большое количество точек, малый шаг.
Тогда все пиксели будут засвечены и не будет пропусков.