Поиск:

Ответ в темуСоздание новой темы Создание опроса
> Профиль Жуковского 
:(
    Опции темы
Merhaba
Дата 24.11.2013, 19:33 (ссылка) | (нет голосов) Загрузка ... Загрузка ... Быстрая цитата Цитата


Шустрый
*


Профиль
Группа: Участник
Сообщений: 93
Регистрация: 23.4.2011

Репутация: нет
Всего: нет



Добрый вечер! Помогите Пожалуйста исправить код. Рисует профиль Жуковского в левом верхнем углу экрана, а надо в центре.
Код

logical(4)::d
integer(4)::status
real(8)::Xmin, Xmax, Ymin, Ymax, p, q 
real(4)::beta1
type (windowconfig) wc
type (wxycoord) wxy
 
 
d=GETWINDOWCONFIG(wc); p=wc.numxpixels; q=wc.numypixels
window.Xmin=-0.5d0*line_str.length*L_ref;     Xmin=window.Xmin
window.Xmax=+0.5d0*line_str.length*L_ref;     Xmax=window.Xmax
window.Ymin=-0.5d0*line_str.length*L_ref*q/p; Ymin=window.Ymin
window.Ymax=+0.5d0*line_str.length*L_ref*q/p; Ymax=window.Ymax
 
status = SETWINDOW(.TRUE., Xmin, Ymin, Xmax, Ymax)
 
    call clearscreen($gclearscreen)
    status=setcolor(15)
    status=rectangle_w($gfillinterior,Xmin,Ymax,Xmax,Ymin)
    status=setcolor(ellips%color)
 
        pi=4.*atan(1.)
        c=100 
        a=110 
        beta1=0.02
        step=2*pi/200
        aa=a*sin(beta1)
        bb=a*cos(beta1)-c
        xm=sqrt(aa*aa+bb*bb)
        delta=atan(aa/bb)
        
 
do o=0.,2*pi,step
 
        x=xm*cos(delta)+a*cos(o)
        y=xm*sin(delta)+a*sin(o)
        r=x*x+y*y
        if (r.ne.0.) then
        xi=x*(1.+c*c/r)
        et=y*(1.-c*c/r)
 
        status=setpixel(xi,et)
 
        end if
 
        end do 


    
PM MAIL   Вверх
  
Ответ в темуСоздание новой темы Создание опроса
1 Пользователей читают эту тему (1 Гостей и 0 Скрытых Пользователей)
0 Пользователей:
« Предыдущая тема | Fortran | Следующая тема »


 




[ Время генерации скрипта: 0.0875 ]   [ Использовано запросов: 21 ]   [ GZIP включён ]


Реклама на сайте     Информационное спонсорство

 
По вопросам размещения рекламы пишите на vladimir(sobaka)vingrad.ru
Отказ от ответственности     Powered by Invision Power Board(R) 1.3 © 2003  IPS, Inc.