Dexp a75: Радар-детектор DEXP A75 | Радар-детекторы | Автотовары | Каталог

Информация на этом веб-сайте предоставляется «как есть, как доступно» без каких-либо гарантий. DeviceSpecifications не несет ответственности за какие-либо упущения, неточности или другие ошибки в публикуемой информации.Все гарантии в отношении этой информации не предоставляются. Воспроизведение любой части этого веб-сайта полностью или частично, в любой форме или на любом носителе без предварительного письменного разрешения запрещено. Торговые марки, марки и логотипы производителей устройств, программного обеспечения, оборудования и т. Д. Являются собственностью соответствующих владельцев.

политика в отношении файлов cookie | Условия эксплуатации
о нас | свяжитесь с нами | новости | отзывы
© 2021 технические характеристики устройств.ком

Вармт! DEXP Ixion ML450 Super Force Sag Høj kvalitet, Exklusiv flip læder phone taske cover tilfældet med Передний слот для слайд-корта

Velkommen til vores Butik.vi er professionalelle factory sælger, all mobiltelefon

sagen er håndlavet .Det vil 100% Perfekt Match for Kamera / Flash / Volumen Huller.

Велкоммен Прямая поставка и отгрузка комбайнерами, работа для хранения!

Varmt! DEXP Ixion ML450 Super Force Sag Høj kvalitet flip læder phone taske cover tilfældet for DEXP Ixion ML450 Super Force med Передний слот для слайд-корта

Forsendelse

1.Kina Post Almindelige Lille Pakke Plus for de små rækkefølge (rækkefølge

ОБРАТНАЯ СВЯЗЬ

1.Vi vil efterlade и 5-STJERNET ПОЛОЖИТЕЛЬНЫЙ ОТЗЫВ til all købere.

2.Din tilfredshed er vores højeste priority, og vores mål er at gøre dig så glad som muligt, samtidig med at de handler med os.

3.Vi er afhængige af din tilfredshed at lykkes.Derfor er din feedback er ekstrem vigtig for os.

Hvis du ikke er tilfreds med dit køb, skal ikke efterlader nogen negativ feed back, før du kontakter os.

4.Vi lover, at vi vil løse проблема для того, чтобы копать så hurtigt som vi kan.

FAQ

F … Produkter kan tilmelde hvad tid?

ru >>> Нормальный тидер 20-60 дней; nogle gange er det, ifølge сказал бэчендлингстид кёбер лэнд; vi køber du fuldføre forsendelse inden for 3 dage (100% bestand)

F … Det er et særligt hylster det?

Ru >>> ja! Hvis du ikke svarer til de huller, vi fuld refundering

F … prisen er en smule dyrt

Ru >>> Распространять все продукты фабрики для продажи по адресу

Ф…Kan ikke finde det ønskede telefonnummer, sag

En >>> kan Du fortælle os; de af os, der er all handmad

F … Billedet ser så kedeligt

En >>> for at sikre ægtheden af ​​produktet; Alle Disse Billeder, Du Har Taget Sælgeren (800 пикселей)

Теги: dexp ixion ml450, afrikanske george, ml450, ixion ml450, dexp glas, sagen dexp ml450, dexp, dexp ml450, brude-shirt, sr125.

5 ядер Originais! DEXP Ixion MS650 Caso da Moda de Luxo Ultra-fino, de Couro de Alta qualidade Exclusiva Caso DEXP Ixion MS650

Bem-vindo à loja

nós somos cometidos para Fornecer os Melhores Produtos e Serviços, o guia da marca:

Усадо DEXP Ixion MS650

Материал Do Produto

De Alta Qualidade Do Couro Flip, микрофибра

Função

1.Material de alta qualidade, ambiental, durável, anti-derrapante

2.Pode ser lavado

3. Судейра, арранхойнс, колизы, tocar Confortável sensação de

4. 100% Совершенное сочетание для камеры / МИКРОФОН / Громкость / Potência Buraco

Cor Do Produto

Королевское шампанское Preto azul, ouro de Rosa-cor-de-Rosa Toque

Handfeel Macio

Embalagem De Varejo

Couro do plutônio de telefone caso da embalagem de revenda

Tempo de entrega

3-5 дней до полуночи

Método de envio

Pacote, HK Post, EMS China Post, DHL, e-pacote de viagem

Serviço

Варехо, Атакадо,

Olá a todos

quente aviso:

1.Se Você comprar mais de 10 computadores, temos atacado link com desconto, acesse a nossa loja para encontrar

2.Se voiceê é um atacadista, o preço pode falar.bem-vindo a contactar-nos.

3.Возможность подключения к телефону в режиме онлайн или по индивидуальному заказу для смартфона, номер телефона

4. Estamos de acordo com o local do furo de seu telefone para soco, Cada caso de telefone é

profissional personalizado, talvez a gente não pode exibir a imagem do local do furo, Mas nós

promessa de que cada caso de telefone é 100% com o seu telefone de Câmera / MTC

/ Volume / Se não for, damos-lhe um reembolso total.Обригадо !!!

nós somos cometidos para Fornecer os Melhores Produtos e Serviços, o guia da marca, 100% de Satisfação

Mais Produtos na nossa, por Favor, vá em:

Этикетки: caso dexp g253, es950, dexp ms650, caso bs155 dexp, samsung s20 edição de fã, ixion, caso dexp, dexp tundra, dexp bl350 filme, caso mi3.

Влияние дексаметазона и троглитазона на чувствительность к инсулину….

Контекст 1

… дексаметазон в высоких дозах (2 мг два раза в день) в течение 3-4 дней оказал сильное влияние на снижение толерантности к глюкозе и чувствительности к инсулину. В присутствии дексаметазона уровни глюкозы и инсулина натощак и через 2 часа после перорального введения глюкозы 75 г были значительно увеличены (таблица 2), как и площади под соответствующими кривыми глюкозы и инсулина (рис. 1A и B). ) по сравнению с исходным исследованием. Из 10 субъектов 4 соответствовали диагностическим критериям диабета, а еще 5 продемонстрировали нарушение толерантности к глюкозе при приеме дексаметазона.Одной из причин такого влияния на толерантность к глюкозе было значительное снижение на 34% максимальной инсулино-стимулированной GDR у лиц, получавших дексаметазон, по сравнению с исходными значениями (рис. …

Контекст 2

… лечились одним троглитазоном (400 мг / день) в течение 4-6 недель, и метаболические эффекты тиазолидиндиона также показаны в таблице 2. Уровни инсулина и глюкозы после перорального введения глюкозы существенно не изменились по сравнению с исходным уровнем после введения одного троглитазона.Тем не менее, троглитазон увеличивал чувствительность, что проявлялось в увеличении GDR на 20% по сравнению с исходным уровнем (P 0,001), как показано на рис. 2. Субъекты набирали в среднем 1,7 кг в течение 4-6 недель при приеме этого лекарства. Других нежелательных эффектов не было, и повышение печеночных трансаминаз не было …

Контекст 3

… исследования были завершены через 4-6 недель приема троглитазона, дексаметазон был повторно введен (2 мг два раза в день в течение 4 дней) в то время как тиазолидиндион был продолжен.Когда данные у субъектов, получавших комбинацию дексаметазона и троглитазона, сравнивали с одним дексаметазоном, было обнаружено, что предварительная обработка троглитазоном резко блокирует пагубные эффекты дексаметазона на толерантность к глюкозе (рис. 1, таблица 2). Кроме того, троглитазон полностью предотвращал индуцированную глюкокортикоидами инсулинорезистентность, поскольку максимальные GDR у субъектов, получавших комбинацию дексаметазона и троглитазона, были аналогичны исходным значениям (рис. 2). Однако средний GDR не увеличивался по сравнению с исходным уровнем у субъектов, принимавших оба препарата, как это наблюдалось с одним троглитазоном.Воздействие на СЖК. Не было обнаружено значительных изменений уровней FFA натощак, уровней триглицеридов или липидных панелей (таблица 2). Однако значительные эффекты наблюдались в подавлении FFAs во время гиперинсулинемико-эугликемического зажима, показателя антилиполитического эффекта инсулина (2 мг два раза в день) за последние 3 дня (дексаметазон плюс троглитазон) (OE). Данные представляют собой средние значения SE для 10 субъектов. на жировой ткани. Подавление FFA инсулином было снижено на 22% под действием дексаметазона и существенно не изменилось при использовании одного троглитазона по сравнению с исходными значениями.Однако троглитазон полностью предотвратил нарушение подавляющей способности FFA дексаметазоном, как показано на рис. 3. Мы не обнаружили корреляции между уровнями FFA натощак и либо GDR в подгруппах лечения, либо изменениями GDR, вызванными дексаметазоном и троглитазоном (данные не показаны. ). Однако, когда данные были проанализированы до и после лечения дексаметазоном, мы обнаружили, что эффект этого препарата на поглощение глюкозы был пропорционален его влиянию на подавляемость FFA (рис. 4A). Напротив, способность троглитазона увеличивать утилизацию глюкозы по сравнению с исходным уровнем не коррелировала с эффектами подавления FFA (рис.4Б). Даже в этом случае улучшение GDR у субъектов, получавших троглитазон плюс дексаметазон, по сравнению с одним дексаметазоном, сильно коррелировало с изменениями в подавляющей способности FFA (рис. 4C). Таким образом, характер изменений подавляемости СЖК был связан с уровнем глюкокортикоид-индуцированной резистентности к инсулину и ее отменой троглитазоном, но не был связан с повышением чувствительности к инсулину в результате действия тиазолидиндиона как такового. Воздействие на лептин. Уровни лептина натощак были увеличены 2.В 2 раза под действием дексаметазона и не изменялись одним троглитазоном по сравнению с исходными значениями (рис. 5). Однако предварительная обработка троглитазоном блокировала удвоение значений лептина под действием дексаметазона и возвращала значения к исходному уровню. Уровни лептина натощак не коррелировали с GDR, уровнем глюкозы или инсулина натощак или 2-часовым OGTT или уровнями FFA (данные не …

Контекст 4

… исследования были завершены через 4-6 недель троглитазона был повторно введен дексаметазон (2 мг б.я бы. в течение 4 дней), в то время как продолжали прием тиазолидиндиона. Когда данные у субъектов, получавших комбинацию дексаметазона и троглитазона, сравнивали с одним дексаметазоном, было обнаружено, что предварительная обработка троглитазоном резко блокирует пагубные эффекты дексаметазона на толерантность к глюкозе (рис. 1, таблица 2). Кроме того, троглитазон полностью предотвращал индуцированную глюкокортикоидами резистентность к инсулину, поскольку максимальные GDR у субъектов, получавших комбинацию дексаметазона и троглитазона, были аналогичны исходным значениям (рис.2). Однако средний GDR не увеличивался по сравнению с исходным уровнем у субъектов, принимавших оба препарата, как это наблюдалось с одним троглитазоном. Воздействие на СЖК. Не было обнаружено значительных изменений уровней FFA натощак, уровней триглицеридов или липидных панелей (таблица 2). Однако значительные эффекты наблюдались в подавлении FFAs во время гиперинсулинемико-эугликемического зажима, показателя антилиполитического эффекта инсулина (2 мг два раза в день) за последние 3 дня (дексаметазон плюс троглитазон) (OE). Данные представляют собой средние значения SE для 10 субъектов.на жировой ткани. Подавление FFA инсулином было снижено на 22% под действием дексаметазона и существенно не изменилось при использовании одного троглитазона по сравнению с исходными значениями. Однако троглитазон полностью предотвратил нарушение подавляющей способности FFA дексаметазоном, как показано на рис. 3. Мы не обнаружили корреляции между уровнями FFA натощак и либо GDR в подгруппах лечения, либо изменениями GDR, вызванными дексаметазоном и троглитазоном (данные не показаны. ). Однако при анализе данных до и после лечения дексаметазоном мы обнаружили, что влияние этого препарата на поглощение глюкозы было пропорционально его влиянию на подавляемость FFA (рис.4А). Напротив, способность троглитазона увеличивать утилизацию глюкозы по сравнению с исходным уровнем не коррелировала с эффектами подавления FFA (рис. 4B). Даже в этом случае улучшение GDR у субъектов, получавших троглитазон плюс дексаметазон, по сравнению с одним дексаметазоном, сильно коррелировало с изменениями в подавляющей способности FFA (рис. 4C). Таким образом, характер изменений подавляемости СЖК был связан с уровнем глюкокортикоид-индуцированной резистентности к инсулину и ее отменой троглитазоном, но не был связан с повышением чувствительности к инсулину в результате действия тиазолидиндиона как такового.Воздействие на лептин. Уровни лептина натощак были увеличены дексаметазоном в 2,2 раза и не изменились одним троглитазоном по сравнению с исходными значениями (рис. 5). Однако предварительная обработка троглитазоном блокировала удвоение значений лептина под действием дексаметазона и возвращала значения к исходному уровню. Уровни лептина натощак не коррелировали с GDR, уровнем глюкозы или инсулина натощак или 2-часовым OGTT или уровнями FFA (данные не …

Контекст 5

… исследования были завершены через 4-6 недель троглитазона был повторно введен дексаметазон (2 мг б.я бы. в течение 4 дней), в то время как продолжали прием тиазолидиндиона. Когда данные у субъектов, получавших комбинацию дексаметазона и троглитазона, сравнивали с одним дексаметазоном, было обнаружено, что предварительная обработка троглитазоном резко блокирует пагубные эффекты дексаметазона на толерантность к глюкозе (рис. 1, таблица 2). Кроме того, троглитазон полностью предотвращал индуцированную глюкокортикоидами резистентность к инсулину, поскольку максимальные GDR у субъектов, получавших комбинацию дексаметазона и троглитазона, были аналогичны исходным значениям (рис.2). Однако средний GDR не увеличивался по сравнению с исходным уровнем у субъектов, принимавших оба препарата, как это наблюдалось с одним троглитазоном. Воздействие на СЖК. Не было обнаружено значительных изменений уровней FFA натощак, уровней триглицеридов или липидных панелей (таблица 2). Однако значительные эффекты наблюдались в подавлении FFAs во время гиперинсулинемико-эугликемического зажима, показателя антилиполитического эффекта инсулина (2 мг два раза в день) за последние 3 дня (дексаметазон плюс троглитазон) (OE). Данные представляют собой средние значения SE для 10 субъектов.на жировой ткани. Подавление FFA инсулином было снижено на 22% под действием дексаметазона и существенно не изменилось при использовании одного троглитазона по сравнению с исходными значениями. Однако троглитазон полностью предотвратил нарушение подавляющей способности FFA дексаметазоном, как показано на рис. 3. Мы не обнаружили корреляции между уровнями FFA натощак и либо GDR в подгруппах лечения, либо изменениями GDR, вызванными дексаметазоном и троглитазоном (данные не показаны. ). Однако при анализе данных до и после лечения дексаметазоном мы обнаружили, что влияние этого препарата на поглощение глюкозы было пропорционально его влиянию на подавляемость FFA (рис.4А). Напротив, способность троглитазона увеличивать утилизацию глюкозы по сравнению с исходным уровнем не коррелировала с эффектами подавления FFA (рис. 4B). Даже в этом случае улучшение GDR у субъектов, получавших троглитазон плюс дексаметазон, по сравнению с одним дексаметазоном, сильно коррелировало с изменениями в подавляющей способности FFA (рис. 4C). Таким образом, характер изменений подавляемости СЖК был связан с уровнем глюкокортикоид-индуцированной резистентности к инсулину и ее отменой троглитазоном, но не был связан с повышением чувствительности к инсулину в результате действия тиазолидиндиона как такового.Воздействие на лептин. Уровни лептина натощак были увеличены дексаметазоном в 2,2 раза и не изменились одним троглитазоном по сравнению с исходными значениями (рис. 5). Однако предварительная обработка троглитазоном блокировала удвоение значений лептина под действием дексаметазона и возвращала значения к исходному уровню. Уровни лептина натощак не коррелировали с GDR, уровнями глюкозы или инсулина натощак или 2-часовой OGTT, или уровнями FFA (данные не …

Контекст 6

… оценивали эффекты глюкокортикоида dexa — метазон и тиазолидиндион троглитазон и впервые исследовали взаимодействие между этими препаратами в отношении множества метаболических параметров у людей.Наши результаты подтверждают предыдущие исследования, демонстрирующие, что глюкокортикоиды вызывают периферическую резистентность к инсулину (2-4) у здоровых людей, и они демонстрируют выраженный эффект стероидов на снижение максимальной чувствительности к инсулину, что не было ясно видно в предыдущем отчете. (2). Это отражает нарушение стимуляции захвата глюкозы в скелетных мышцах, которое, как было показано, является результатом дефектов в привлечении переносчиков глюкозы GLUT4 к поверхности клеток на животных (12,13) ​​и клеточных (14-17) моделях.Периферическая инсулинорезистентность вместе с инсулинорезистентностью на уровне печени (2-4) приводит к нарушению толерантности к глюкозе, что проявляется в текущих данных более высокими концентрациями глюкозы и инсулина после перорального введения глюкозы. Глюкокортикоиды могут также увеличивать липолиз и концентрацию циркулирующих FFA и предшествующее введение (2 мг два раза в день) (Dex), после лечения только троглитазоной (400 мг / день) в течение 4-6 недель (Trog) и после лечения троглитазоном. (400 мг / день) в течение 4-6 недель и дексаметазон (2 мг b.i.d.) за последние 4 дня (DexTrog). Данные представляют собой максимально стимулированный GDR во время клэмп-исследований с гиперинсулинемией-эугликемией и нормализованы по килограммам безжировой массы тела (LBM), оцененным DEXA. Данные представляют собой средние значения SE для 10 субъектов. * P <0,01 по сравнению с исходным уровнем; ** P <0,001 по сравнению с исходным уровнем. авторы связывают этот эффект с патогенезом индуцированной глюкокортикоидами инсулинорезистентности в скелетных мышцах (5,34,35). Наши данные подтверждают эту идею, поскольку дексаметазон снижает способность инсулина подавлять уровни FFA, и этот эффект хорошо коррелирует со снижением максимально стимулированного поглощения глюкозы.Хотя дексаметазон не влиял на уровни FFA натощак у наших субъектов, нарушение подавляемости FFA предсказуемо приведет к увеличению концентрации FFA после еды. Таким образом, хроническое повышение постпрандиальных профилей FFA может способствовать инсулинорезистентности в скелетных мышцах (27) (28) (29). Связь между способностью инсулина подавлять FFA и стимулировать усвоение глюкозы мышцами была ранее продемонстрирована Ривеном и его коллегами (45,46) на выборке здоровых ...

Контекст 7

… оценили эффекты глюкокортикоидов дексаметазона и тиазолидиндиона троглитазона и впервые исследовали взаимодействие между этими препаратами в отношении множества метаболических параметров у людей. Наши результаты подтверждают предыдущие исследования, демонстрирующие, что глюкокортикоиды вызывают периферическую резистентность к инсулину (2-4) у здоровых людей, и они демонстрируют выраженный эффект стероидов на снижение максимальной чувствительности к инсулину, что не было ясно видно в предыдущем отчете. (2).Это отражает нарушение стимуляции захвата глюкозы в скелетных мышцах, которое, как было показано, является результатом дефектов в привлечении переносчиков глюкозы GLUT4 к поверхности клеток на животных (12,13) ​​и клеточных (14-17) моделях. Периферическая инсулинорезистентность вместе с инсулинорезистентностью на уровне печени (2-4) приводит к нарушению толерантности к глюкозе, что проявляется в текущих данных более высокими концентрациями глюкозы и инсулина после перорального введения глюкозы. Глюкокортикоиды могут также увеличивать липолиз и концентрацию циркулирующих FFA, а предыдущие (2 мг b.id) введение (Dex), после лечения только троглитазоном (400 мг / день) в течение 4-6 недель (Trog) и после лечения троглитазоном (400 мг / день) в течение 4-6 недель и дексаметазоном (2 мг два раза в день) за последние 4 дня (DexTrog). Данные представляют собой максимально стимулированный GDR во время клэмп-исследований с гиперинсулинемией-эугликемией и нормализованы по килограммам безжировой массы тела (LBM), оцененным DEXA. Данные представляют собой средние значения SE для 10 субъектов. * P <0,01 по сравнению с исходным уровнем; ** P <0,001 по сравнению с исходным уровнем. авторы связывают этот эффект с патогенезом индуцированной глюкокортикоидами инсулинорезистентности в скелетных мышцах (5,34,35).Наши данные подтверждают эту идею, поскольку дексаметазон снижает способность инсулина подавлять уровни FFA, и этот эффект хорошо коррелирует со снижением максимально стимулированного поглощения глюкозы. Хотя дексаметазон не влиял на уровни FFA натощак у наших субъектов, нарушение подавляемости FFA предсказуемо приведет к увеличению концентрации FFA после еды. Таким образом, хроническое повышение постпрандиальных профилей FFA может способствовать инсулинорезистентности в скелетных мышцах (27) (28) (29).Связь между способностью инсулина подавлять СЖК и стимулировать усвоение глюкозы мышцами была ранее продемонстрирована Ривеном и его коллегами (45,46) на выборке здоровых ...

C ———- ————————————————— ———- C Подгонка магнитной восприимчивости как функции температуры. C Ссылка: Г. В. Р. Чандрамули, К. Балагопалакришна. С. М. В. Раджасекхаран и П. Т. Манохаран C Computers Chem.20 (1996), 353-358. C Дата: сентябрь 1994 г., изменено в апреле 2002 г. C C ————————————————- ——————— C C Максимальное количество точек данных — 300 C Максимальное количество параметров минимизации — 12 C Максимальное количество членов в функции — 10 C Файл входных данных содержит C NP — количество баллов C T (I) — Температура, XM (I) — Наблюдаемая восприимчивость на ион C при T (I), включая диамагнитную поправку.Бегу от 1 до НП. C Остальные данные выводятся на терминал. C g-значение может быть фиксированным или рассматриваться как параметр минимизации. C Рассматриваются следующие функции: C 1. Кюри 4. Блини-Бауэр 7. Модифицированный BF 10. DIM. C 2. Нулевое поле 5. Намагниченность 8. Димеры C 3. Изинг 6. Боннер-Фишер 9. Другое C Другое: 1. Определяется пользователем 3. Тетрамер Хэтфилда (S = 1/2) C 2. Тример (S = 1/2) 4.Димерные цепи Кана C Пользовательская функция требует модификации подпрограмм. C USRPAR и USER. См. Комментарии в этих подпрограммах. C Междимерные взаимодействия учитываются молекулярным полем прибл. C член Вейсса необязательно включен для члена Кюри. C Можно комбинировать вышеуказанные функции (термины), чтобы выбрать выражение C для расчета восприимчивости. Меню предварительно выбранных Комбинации C представлены для выбора желаемого выражения.C Также предусмотрены другие комбинации (вариант 99). C C Сначала отображается график наблюдаемых и расчетных данных. C Пользователь может выбрать один из трех вариантов сюжета — C SPLOT для графического экрана. PPLOT — для текстового экрана C HPPLOT для создания выходного файла с инструкциями HP по печати C Начальные параметры выражения могут быть интерактивно C изменен. C Выходные данные записываются в SUSCEP.OUT, только если были выполнены итерации.C Файл SUSCEP.HP создается, только выбрана опция HPPLOT. C EXPLMT — машинно-зависимый предел экспоненты по основанию e (= 2,714) C C Соглашение о значении J: синглет-триплетный разрыв рассматривается как 2J. C Единица измерения J — волновое число. Вычисляет восприимчивость на магнитный ион. C C ————————————————- ——————— Программа C использует графические процедуры Microsoft 5.0 для экранной графики. C Протестировано с компилятором Microsoft FORTRAN вер.5.0 C Подпрограмма NELMIN выбирается из СТАТИСТИЧЕСКИХ АЛГОРИТМОВ 47. C В него были внесены незначительные изменения для извлечения оценок ошибок. C Оценка погрешности выполняется в соответствии с C Грегори Р. Филлипс и Эдвард М. Эйринг C Анал. Chem. 1988, vol 60, pp 738-741 и ссылки в них. C Стандартный GJRD был взят из НАУЧНОГО ПОДПРОЧНОГО ПАКЕТА. C ————————————————- ——————— ВКЛЮЧИТЕ «FGRAPH.FI ‘ C ————————————————- ——————— НЕПРИЯТНОЕ РЕАЛЬНОЕ * 8 (A-H, O-Z) ХАРАКТЕР ПАР * 8, ИНФАЙЛ * 20, ОПЦИЯ * 1 ДВОЙНАЯ ТОЧНОСТЬ МИН (12), TC (300), XIC (300), XMO (300) РАЗМЕР СВПАР (12), СВПЭП (12), СИГМА (12) ОБЩИЙ / BLK1 / T (300), XM (300), DA, NP, IENTRY ОБЩИЙ / BLK2 / HCBIK, B2NBIK, G, G2B2NK, HCNB2, HCNGB2 ОБЩИЙ / BLK2A / BK, BN, GBK, GBN, TIP ОБЩИЙ / BLK3 / S (10), Z (10), NS (10) ОБЩИЙ / ЧЕР4 / ПУСК (12), ШАГ (12), ПАР (12) ОБЩИЙ / BLK5 / IN, IO ОБЩИЙ / BLK6 / ITYPE, NFN, NF (10), NPAR ОБЩИЙ / LMTS / EXPLMT, NPRLMT, NFNLMT C Раздел ввода (также смотрите подпрограмму GETPAR) IN = 5 IO = 6 EXPLMT = DLOG (1.0D307) NPRLMT = 12 NFNLMT = 10 C Наблюдаемые точки данных считываются из INFILE ЗАПИСАТЬ (IO, 5) 5 FORMAT (1h2, /// T27, ‘УСТАНОВКА ДАННЫХ ОБ УСТОЙЧИВОСТИ’ / 1 T22, авторы G.V.R.C, C.B.G., M.V.R., & P.T.M ‘/ 2 T22, «Ссылка: Компьютеры и химия» /// 3 ‘Имя файла экспериментальных данных?’) ЧИТАТЬ (IN, 10) INFILE 10 ФОРМАТ (A) ОТКРЫТЬ (БЛОК = 1, ФАЙЛ = ИНФАЙЛ, СОСТОЯНИЕ = ‘СТАРЫЙ’) ПРОЧИТАТЬ (1, *) NP ЧИТАТЬ (1, *) (T (I), XMO (I), I = 1, NP) ПИШИТЕ (IO, *) «Нет.наблюдаемых точек данных: ‘, NP ЗАПИСАТЬ (IO, *) 1 ‘Температура (K) X (M) Температура (K) X (M) Температура (K) X (M)’ ЗАПИСАТЬ (IO, 11) (T (I), XMO (I), I = 1, NP) 11 ФОРМАТ (3 (F7.2, F9.6,5X)) ЗАКРЫТЬ (1) WRITE (IO, *) ‘Если указанные выше значения X (M) НЕ относятся к магнитному иону,’ WRITE (IO, *) ‘умножьте их на коэффициент, чтобы получить значение X (M) на ион’ WRITE (IO, *) ‘См. Следующее меню. Нажмите клавишу Enter … ‘ ПРОЧИТАЙТЕ (IN, 10) ВАРИАНТ C Настройки по умолчанию DIA = 0.D0 TIP0 = 0.D0 ФАКТОР = 1.D0 NTH = 0 NEQN = 24 G = 0.D0 12 NPAR = 0 C Настоящее главное меню ЗАПИСАТЬ (IO, 20) DIA, TIP0 20 FORMAT (‘A. Диамагнитное значение, вычитаемое из наблюдаемого’ 2 ‘data =’, 1PE10.3 / ‘B. Значение TIP, которое нужно добавить к вычисленному’, 3 ‘Susceptibility =’, 1PE10.3) ЗАПИСАТЬ (IO, 21) FACTOR, NEQN 21 FORMAT (‘C. Умножьте наблюдаемые данные на множитель =’, F5.2, / 5 ‘D. Покажите настоящие ценности’ / 6 ‘E. Номер выражения =’, I2, / 7 ‘(Двузначный номер MN (01-99). Первая цифра M выбирает’, 8 ‘тип выражения.’ / 9 ‘Вторая цифра N выбирает функцию F (N). ‘, 9 ‘MN = 99 для нового выражения)’) ЗАПИСАТЬ (IO, 25) 25 FORMAT (1X, 75 (‘-‘) / ‘M Выражения для X (cal)’ 1, T35, ‘| N F (N) ‘, T55,’ Consts ‘, T62,’ Переменные ‘/ 2 1Х, 75 (‘-‘) / 3 ‘0 Простой F (N)’, T35, ‘| 1.Кюри ‘, T55,’ S ‘/ 4 ‘1 F (N) + DIM’, T35, ‘| 2. Нулевое поле, T55, S, T62, D / 5 ‘2 x F (N) + (1-x) Кюри’, T35, ‘| 3. Изинг », T62,« J-ISNG »/ 6 ‘3 x F (N) + (1-x) Z-F’, T35, ‘| 4. Блини-Бауэр, T62, J-BB / 7 ‘4 x F (N) + (1-x) Кюри + DIM’, T35, 8 ‘| 5. Намагничивание, T55, B, T62, J-MAGN) ЗАПИСАТЬ (IO, 30) G 30 ФОРМАТ (‘5 x F (N) + (1-x) Z-F + DIM’, T35, 1 ‘| 6. Боннер-Фишер, T62, «J-BF» / 2 ‘9 Ничего из вышеперечисленного’, T35, 3 ‘| 7.Модифицированные BF ‘, T62,’ J-MBF, Alpha ‘/ 4 ‘x — это COEFF1 (‘, E12.5) 100 F = CON * S * (S + 1) 200 ВОЗВРАТ КОНЕЦ C Ising Eq. ВСПОМОГАТЕЛЬНЫЙ ВЫХОД (A, T, F) НЕПРИЯТНОЕ РЕАЛЬНОЕ * 8 (A-H, O-Z) ОБЩИЙ / BLK2 / HCBIK, B2NBIK, G, G2B2NK, HCNB2, HCNGB2 ОБЩИЙ / BLK5 / IN, IO ОБЩИЙ / LMTS / EXPLMT, NPRLMT, NFNLMT AJ = A * HCBIK / T ЕСЛИ (AJ.GT.EXPLMT) ЗАПИСАТЬ (IO, *) 1 ‘Предупреждение об экспоненте в ISING: AJ =’, AJ X = DEXP (AJ) F = G2B2NK / (2.D0 * Т) * Х / (1 + Х) ВОЗВРАЩЕНИЕ КОНЕЦ Выражение лица C Блини Бауэра ВСПОМОГАТЕЛЬНЫЙ BB (AJ, T, F) НЕПРИЯТНОЕ РЕАЛЬНОЕ * 8 (A-H, O-Z) ОБЩИЙ / BLK2 / HCBIK, B2NBIK, G, G2B2NK, HCNB2, HCNGB2 ОБЩИЙ / BLK5 / IN, IO ОБЩИЙ / LMTS / EXPLMT, NPRLMT, NFNLMT X = 2.D0 * HCBIK * AJ / T ЕСЛИ (DABS (X) .GT.EXPLMT) + WRITE (IO, *) ‘Предупреждение об экспоненте в BB -> X =’, X TOP = 2,0 * G2B2NK / T БОТ = 3.D0 + DEXP (-X) F = TOP / BOT ВОЗВРАЩЕНИЕ КОНЕЦ C Обобщенная намагниченность ур. ВСПОМОГАТЕЛЬНЫЙ МАГН (AJ, B, T, F) НЕПРИЯТНОЕ РЕАЛЬНОЕ * 8 (A-H, O-Z) ОБЩИЙ / BLK2 / HCBIK, B2NBIK, G, G2B2NK, HCNB2, HCNGB2 ОБЩИЙ / BLK2A / BK, BN, GBK, GBN, TIP ОБЩИЙ / BLK5 / IN, IO ОБЩИЙ / LMTS / EXPLMT, NPRLMT, NFNLMT X = -2.D0 * HCBIK * AJ / T ЕСЛИ (DABS (X) .GT.EXPLMT) + WRITE (IO, *) ‘Предупреждение о экспоненте в MAGN -> X =’, X СРОК = GBK * B / T TOP = GBN * SINH (СРОК) БОТ = DEXP (X) +2.D0 * COSH (СРОК) + 1.D0 F = TOP / BOT / B ВОЗВРАЩЕНИЕ КОНЕЦ C Боннер-Фишер уравнение. ПОДРОБНЕЕ BF (AJ, T, F) НЕПРИЯТНОЕ РЕАЛЬНОЕ * 8 (A-H, O-Z) ОБЩИЙ / BLK2 / HCBIK, B2NBIK, G, G2B2NK, HCNB2, HCNGB2 ОБЩИЙ / BLK5 / IN, IO X = HCBIK * DABS (AJ) / T Х2 = Х * Х Х3 = Х2 * Х ANUM = 0,25 + 0,149445 * X + 0,30094 * X2 DENO = 1 + 1,9862 * X + 0,68854 * X2 + 6,0626 * X3 F = G2B2NK / T * ANUM / DENO ВОЗВРАЩЕНИЕ КОНЕЦ c Модифицированное уравнение Боннера-Фишера c (эта процедура изменена в апреле 2002-мвр) c Ссылка: W.E. Hatfield, J. Appl. Phys. 52 (1981), 1985 и c «Молекулярный магнетизм», О. Кан …, стр.263-265. ВСПОМОГАТЕЛЬНЫЙ MBF (AJ, AL, T, F) НЕПРИЯТНОЕ РЕАЛЬНОЕ * 8 (A-H, O-Z) ОБЩИЙ / BLK2 / HCBIK, B2NBIK, G, G2B2NK, HCNB2, HCNGB2 ОБЩИЙ / BLK5 / IN, IO AL2 = AL * AL AL3 = AL2 * AL AL4 = AL3 * AL ЕСЛИ (AL.LE.0.4) ТО A1 = 0,25 В = -. 12587 + .22752 * AL C = .019111-.13307 * AL + .509 * AL2-1,3167 * AL3 + 1.0081 * AL4 D = 0,10772 + 1,4192 * AL E = -. 0028521-.42346 * AL + 2,1953 * AL2-.82412 * AL3 F = .37754-6.702201E-2 * AL + 6.9805 * AL2-21.678 * AL3 + 15.838 * AL4 ЕЩЕ A1 = 0,25 В = -. 13695 + .26387 * AL C = 0,017025–0,12668 * AL + 0,49113 * AL2-1.1977 * AL3 + .87257 * AL4 D = 0,070509 + 1,3042 * AL E = -. 0035767-.40837 * AL + 3,4862 * AL2-.73888 * AL3 F = 0,36184-0,065528 * AL + 6,65875 * AL2-20.945 * AL3 + 15,42504 * AL4 ENDIF X = HCBIK * (- AJ) / Т Х2 = Х * Х Х3 = Х2 * Х АНУМ = A1 + B * X + C * X2 ДЕНО = 1 + D * X + E * X2 + F * X3 F = G2B2NK / T * ANUM / DENO ВОЗВРАЩЕНИЕ КОНЕЦ C Гейзенберг-интрадимерное взаимодействие для различных спиновых систем ВСПОМОГАТЕЛЬНЫЙ ДИМЕР (S1, S2, AJ, T, F) НЕПРИЯТНОЕ РЕАЛЬНОЕ * 8 (A-H, O-Z) ОБЩИЙ / BLK2 / HCBIK, B2NBIK, G, G2B2NK, HCNB2, HCNGB2 ОБЩИЙ / BLK5 / IN, IO ОБЩИЙ / LMTS / EXPLMT, NPRLMT, NFNLMT РАЗМЕР I (10), E (10) AJBIKT = HCBIK * AJ / T I1 = 2 * S1 + 0.01 I2 = 2 * S2 + 0,01 ITOT = I1 + I2 + 1 IMIN = IABS (I1-I2) +1 К = 0 DO 40 LL = IMIN, ITOT, 2 L = LL-1 К = К + 1 I (K) = L ЕСЛИ (LL.EQ.IMIN) IEK = 0 ЕСЛИ (LL.GT.IMIN) IEK = IEK + L СРОК = IEK * AJBIKT ЕСЛИ (DABS (TERM) .GT. EXPLMT) 1 WRITE (IO, *) ‘Предупреждение о экспоненте в DIMER ->, TERM =’, TERM E (K) = DEXP (СРОК) * LL 40 ПРОДОЛЖИТЬ TOP = 0 BOT = 0 L = K ДО 50 К = 1, L ВЕРХ = ВЕРХ + I (K) * (I (K) +2) * E (K) 50 БОТ = БОТ + E (K) F = G2B2NK / 12.D0 / T * TOP / BOT ВОЗВРАЩЕНИЕ КОНЕЦ C ДРУГИЕ функции C Trimer eq. S1 = S2 = S3 = 1/2, J12 = J23 = J31 = J C Ссылка: Mossoba M. M. et al. Варенье. Chem. Soc. (1980) 102,6864. ПОДРОБНЫЙ ТРИМЕР (A, IPAR, T, F) НЕПРИЯТНОЕ РЕАЛЬНОЕ * 8 (A-H, O-Z) ОБЩИЙ / BLK2 / HCBIK, B2NBIK, G, G2B2NK, HCNB2, HCNGB2 ОБЩИЙ / BLK5 / IN, IO ОБЩИЙ / LMTS / EXPLMT, NPRLMT, NFNLMT РАЗМЕР A (1) IPAR = IPAR + 1 AJ3 = 3.D0 * A (ИПАР) * HCBIK / T ЕСЛИ (AJ3.GT.EXPLMT) ЗАПИСАТЬ (IO, *) 1 ‘Предупреждение об экспоненте в TRIMER: AJ * 3 =’, AJ3 X = DEXP (AJ3) F = G2B2NK / (4.D0 * T) * (1.D0 + 5.D0 * X) / (1.D0 + X) ВОЗВРАЩЕНИЕ КОНЕЦ C Четырехспиновая связь Хэтфилда в кубе (S1 = S2 = S3 = S4 = 1/2) C Ссылка: Hall, J. W. et al. Неорг. Chem. (1977) 16, 1572. ВСПОМОГАТЕЛЬНЫЙ HATFLD (A, IPAR, T, F) НЕПРИЯТНОЕ РЕАЛЬНОЕ * 8 (A-H, O-Z) ОБЩИЙ / BLK2 / HCBIK, B2NBIK, G, G2B2NK, HCNB2, HCNGB2 ОБЩИЙ / BLK5 / IN, IO ОБЩИЙ / LMTS / EXPLMT, NPRLMT, NFNLMT РАЗМЕР A (1) IPAR = IPAR + 1 AJ1 = A (IPAR) * HCBIK / T IPAR = IPAR + 1 AJ2 = A (IPAR) * HCBIK / T IPAR = IPAR + 1 AJ3 = A (IPAR) * HCBIK / T IPAR = IPAR + 1 AJ4 = A (IPAR) * HCBIK / T AJ1PJ2 = AJ1 + AJ2 AJ1MJ2 = AJ1-AJ2 AJ3PJ4 = AJ3 + AJ4 AJ3MJ4 = AJ3-AJ4 TT1 = AJ1PJ2 + 2.D0 * AJ3PJ4 Q = — (TT1) / 2 Т1 = — (AJ1PJ2-2.D0 * AJ3PJ4) / 2 TT2 = 2.D0 * DSQRT (AJ1MJ2 * AJ1MJ2 + AJ3MJ4 * AJ3MJ4) Т2 = (AJ1PJ2 + TT2) /2.D0 Т3 = (AJ1PJ2-TT2) /2.D0 TT3 = 2.D0 * DSQRT ((AJ1PJ2-AJ3PJ4) ** 2 + 3.D0 * AJ3MJ4 * AJ3MJ4) S1 = (TT1 + TT3) /2.D0 S2 = (TT1-TT3) /2.D0 ЕСЛИ (DABS (Q) .GT.EXPLMT) 1 WRITE (IO, *) ‘Предупреждение об экспоненте в HATFLD -> Q =’, Q Q = DEXP (-Q) T1 = DEXP (-T1) T2 = DEXP (-T2) T3 = DEXP (-T3) S1 = DEXP (-S1) S2 = DEXP (-S2) ТОП = 5.D0 * Q + T1 + T2 + T3 БОТ = 5.D0 * Q + 3.D0 * (T1 + T2 + T3) + S1 + S2 КОН = 2.D0 * G2B2NK / T F = CON * TOP / BOT ВОЗВРАЩЕНИЕ КОНЕЦ C O. Kahn’s eqn. для цепочки димеров C Ссылка: Lorett, F. et al. Неорг. Chem. (1993) 32, 27. ПОДРОБНЕЕ KAHN (A, IPAR, T, F, ITERM) НЕПРИЯТНОЕ РЕАЛЬНОЕ * 8 (A-H, O-Z) ОБЩИЙ / BLK2 / HCBIK, B2NBIK, G, G2B2NK, HCNB2, HCNGB2 ОБЩИЙ / BLK3 / S (10), Z (10), NS (10) ОБЩИЙ / BLK5 / IN, IO ОБЩИЙ / LMTS / EXPLMT, NPRLMT, NFNLMT РАЗМЕР A (1) IPAR = IPAR + 1 AJ = A (IPAR) IPAR = IPAR + 1 AJ2 = A (IPAR) * HCBIK / T ВЫЗОВ ДАЙМЕР (S (ITERM), Z (ITERM), AJ, T, F) SSP1 = F / G2B2NK * 3.D0 * T X = AJ2 * SSP1 U = 1.D0 / DTANH (X) -1.D0 / X F = F * (1.D0 + U) / (1.D0-U) ВОЗВРАЩЕНИЕ КОНЕЦ C Исходные определения для ДРУГИХ уравнений ВСПОМОГАТЕЛЬНЫЙ OTHPAR (NF, I) НЕПРИЯТНОЕ РЕАЛЬНОЕ * 8 (A-H, O-Z) ОБЩИЙ / BLK3 / S (10), Z (10), NS (10) ОБЩИЙ / BLK5 / IN, IO ОБЩИЙ / BLK9 / NOTH 1 ЗАПИСЬ (IO, 5) 5 FORMAT (‘Выберите одну из следующих функций из других функций’ / 1 ‘1.Пользовательская функция ‘/ 1 ‘2. S = 1/2 Трехспиновая система’ / 2 ‘3. S = 1/2 Система Хэтфилда’ / 3 ‘4. Выражение Хана’) НАПИСАТЬ (IO, *) » WRITE (IO, *) ‘Введите номер опции:’ ПРОЧИТАЙТЕ (IN, *) NOTH GOTO (10,20,30,40), NOTH GOTO 1 10 WRITE (IO, *) ‘Пользовательская функция’ ПОЗВОНИТЬ USRPAR (NF, I) GOTO 200 20 WRITE (IO, *) ‘Три вращения S1 = S2 = S3 = 1/2 системы с одним J:’ ВЫЗОВ GETPAR (‘J-Trimr’) GOTO 200 30 WRITE (IO, *) 4-спиновую (S = 1/2) систему «Хэтфилда»: ВЫЗОВ GETPAR (‘J1-Tetra’) ВЫЗОВ GETPAR (‘J2-Tetra’) ВЫЗОВ GETPAR (‘J3-Tetra’) ВЫЗОВ GETPAR (‘J4-Tetra’) GOTO 200 40 WRITE (IO, *) ‘Выражение Кана для цепочки димеров:’ WRITE (IO, *) ‘Введите S1, S2 димера’ ПРОЧИТАЙТЕ (IN, *) S (I), Z (I) ВЫЗОВ GETPAR (‘J-Intra’) ПОЗВОНИТЬ GETPAR (‘J-Inter’) 200 ВОЗВРАТ КОНЕЦ C Выбор функции для ДРУГИХ уравнений.C ITERM необходим, если значение NF (ITERM) используется в любом C функций, используемых в этой программе. ПОДРОБНЕЕ ДРУГИЕ (A, IPAR, T, F, ITERM) НЕПРИЯТНОЕ РЕАЛЬНОЕ * 8 (A-H, O-Z) РАЗМЕР A (1) ОБЩИЙ / BLK5 / IN, IO ОБЩИЙ / BLK6 / ITYPE, NFN, NF (10), NPAR ОБЩИЙ / BLK9 / NOTH GOTO (10,20,30,40), NOTH 10 ПОЗВОНИТЬ ПОЛЬЗОВАТЕЛЮ (A, IPAR, T, F) GOTO 100 20 ТРИМЕР ВЫЗОВА (A, IPAR, T, F) F = F / 3 GOTO 100 30 ЗВОНОК HATFLD (A, IPAR, T, F) F = F / 4 GOTO 100 40 ЗВОНИТЕ КАНУ (A, IPAR, T, F, ITERM) F = F / 2 100 ВОЗВРАТ КОНЕЦ C Обе подпрограммы USRPAR и USER должны быть C изменен для изменения пользовательской функции.C USRPAR определяет имена и начальные значения минимизации Параметры C. Это выполняется оператором CALL GETPAR один раз для C каждым параметром. C C Спин 1-1 димера с изотропным D. Определите параметры. ВСПОМОГАТЕЛЬНЫЙ USRPAR (NF, I) НЕПРИЯТНОЕ РЕАЛЬНОЕ * 8 (A-H, O-Z) ОБЩИЙ / BLK3 / S (10), Z (10), NS (10) ОБЩИЙ / BLK5 / IN, IO C C ———- Измените следующую часть процедуры, чтобы изменить —— C ———- к новому уравнению ———————————— ——— C WRITE (IO, *) ‘S1 = 1, S2 = 1 димерное уравнение.включая D ‘ ПОЗВОНИТЬ GETPAR (‘J-1-1’) ПОЗВОНИТЬ GETPAR (‘D-1-1’) C C ————————————————- ——————— ВОЗВРАЩЕНИЕ КОНЕЦ C Димер Spin 1-1 с изотропным D. Оценка функции ПОДПРОГРАММНЫЙ ПОЛЬЗОВАТЕЛЬ (A, IPAR, T, F) НЕПРИЯТНОЕ РЕАЛЬНОЕ * 8 (A-H, O-Z) ОБЩИЙ / BLK2 / HCBIK, B2NBIK, G, G2B2NK, HCNB2, HCNGB2 ОБЩИЙ / BLK5 / IN, IO ОБЩИЙ / LMTS / EXPLMT, NPRLMT, NFNLMT РАЗМЕР A (1) C ————————————————- ——————— C Эта процедура определяет заданную пользователем функцию, которая может быть C рассматривается как термин в выражении восприимчивости.C A — массив параметров. C IPAR — это последний использованный номер параметра. Если эта рутина C определяет любые новые параметры, IPAR должен быть соответственно C увеличился перед ВОЗВРАЩЕНИЕМ к процедуре вызова. C T — Температура (K) C F — Восприимчивость, рассчитанная по этой методике C HCBIK = hc / k; B2NBIK = beta * beta * N / k С G2B2NK = g * g * beta * beta * N / k C HCNB2 = hc / (N * бета * бета) C HCNGB2 = hc / (N * g * бета * бета) C BK = бета * k BN = бета * N GBK = g * BK GBN = g * BN C Если g — параметр минимизации, g = A (1), иначе g = G.C Условие G SJ = ‘, SJ E1 = DEXP (TJ — 2.D0 * D / 3.D0) E2 = DEXP (TJ + D / 3.D0) E3 = DEXP (SJ — 2.D0 * D) E4 = DEXP (SJ — D) E5 = DEXP (SJ + 2.D0 * D) ВЕРХ = 2.D0 * E1 + 4.D0 * E2 + 6.D0 * E3 + 12.D0 * (E4 + E5) БОТ = 1.D0 + E1 + E3 + 2.D0 * (E2 + E4 + E5) C ————————————————- ——————— C F возвращается к процедуре вызова C ————————————————- ——————— F = CON * TOP / BOT / 2 ВОЗВРАЩЕНИЕ КОНЕЦ C EQN записывает выражение в файл SUSCEP.ВНЕ ПОДПРОГРАММА EQN НЕПРИЯТНОЕ РЕАЛЬНОЕ * 8 (A-H, O-Z) ПЕРСОНАЖ P * 3, PP (10) * 3 ОБЩИЙ / BLK2 / HCBIK, B2NBIK, G, G2B2NK, HCNB2, HCNGB2 ОБЩИЙ / BLK3 / S (10), Z (10), NS (10) ОБЩИЙ / BLK6 / ITYPE, NFN, NF (10), NPAR ОБЩИЙ / BLK9 / NOTH ДАННЫЕ PP / ‘A. ‘,’ Б. ‘,’ C. ‘,’ D. ‘,’ Э. ‘,’ F. ‘ 1, Г. ‘,’ЧАС. ‘,’Я. ‘,’ Дж. ‘/ WRITE (2, *) ‘Тип уравнения, используемый для подбора данных:’ ЕСЛИ (ITYPE.EQ.1) WRITE (2, *) ‘A + B + …’ ЕСЛИ (ITYPE.EQ.2) WRITE (2, *) ‘x (A) + (1-x) (B) + C + …’ ЕСЛИ (ITYPE.EQ.3) WRITE (2, *) 1 ‘y [x (A) + (1-x) (B)] + (1-y) (C) + D + …’ IF (ITYPE.EQ.1) WRITE (2, *) ‘Где,’ IF (ITYPE.EQ.2) WRITE (2, *) ‘Где x — коэффициент, а’ IF (ITYPE.EQ.3) WRITE (2, *) ‘Где x и y — коэффициенты, а’ ДО 100 I = 1, NFN P = PP (I) N = NF (I) ЕСЛИ (N.EQ.1) WRITE (2,5) P, ‘член Кюри’, S (I) 5 FORMAT (1X, 2A, ‘S =’, F3.1) IF (N.EQ.2) WRITE (2,5) P, ‘Термин нулевого поля’, S (I) IF (N.EQ.3) WRITE (2, *) P, ‘Термин Изинга’ ЕСЛИ (N.EQ.4) WRITE (2, *) P, ‘член Блини-Бауэра’ IF (N.EQ.5) WRITE (2, *) P, ‘Термин намагничивания’ ЕСЛИ (N.EQ.6) WRITE (2, *) P, ‘термин Боннера Фишера’ IF (N.EQ.7) WRITE (2, *) P, ‘Модифицированный член Боннера Фишера’ ЕСЛИ (N.EQ.8) WRITE (2,20) P, ‘Интрадимерный член’, S (I), Z (I) 20 ФОРМАТ (1X, 2A, ‘S1 =’, F3.1, ‘S2 =’, F3.1) 30 IF (N.NE.9) GOTO 40 ЕСЛИ (НЕТ.EQ.1) WRITE (2, *) P, ‘Пользовательский термин’ ЕСЛИ (NOTH.EQ.2) WRITE (2, *) P, ‘S = 1/2 тримерный член’ ЕСЛИ (NOTH.EQ.3) WRITE (2, *) P, S = 1/2 тетрамерного члена ‘Хэтфилда’ ЕСЛИ (NOTH.EQ.4) ЗАПИШИТЕ (2,20) P, «член димерной цепи Кана», S (I), Z (I) ЕСЛИ (NOTH.EQ.5) НАПИШИТЕ (2, *) P, термин Heucher ‘ 40 IF (N.EQ.10) WRITE (2, *) P, ‘DIM term’ ЕСЛИ (NS (I) .NE.1) GOTO 100 IF (N.EQ.1) WRITE (2, *) ‘Применена поправка Вейсса (тета)’ IF (N.GT.1) WRITE (2, *) ‘Применена коррекция соседства’ 100 ПРОДОЛЖИТЬ ЕСЛИ (DABS (G).GT. 1.D-6) WRITE (2, *) ‘фиксированное значение g’ IF (.NOT. (DABS (G) .GT. 1.D-6)) WRITE (2, *) ‘значение g изменилось’ ВОЗВРАЩЕНИЕ КОНЕЦ C Установка значений координаты X для графика ВСПОМОГАТЕЛЬНЫЕ XDATA (T, NP, TC, NC, NPLOP, IPLOP) НЕПРИЯТНОЕ РЕАЛЬНОЕ * 8 (A-H, O-Z) REAL * 4 XOBS, XCALC, YOBS, YCALC, XMIN, XMAX, YMIN, YMAX РАЗМЕР T (1), TC (1) ОБЩИЙ / BLK7 / XOBS (300), XCALC (300), YOBS (300), YCALC (300) ОБЩИЙ / BLK8 / XMIN, XMAX, YMIN, YMAX NC = 300 XMIN = T (1) XMAX = Т (1) ДО 10 I = 2, НП ЕСЛИ (XMAX.LT. Т (I)) XMAX = Т (I) ЕСЛИ (XMIN .GT. T (I)) XMIN = T (I) 10 ПРОДОЛЖИТЬ ЕСЛИ (IPLOP.EQ.2) GOTO 40 XINT = (XMAX-XMIN) / (NC-1) DO 30 I = 1, НЗ 30 TC (I) = (I-1) * XINT + XMIN GOTO 60 40 NC = NP DO 50 I = 1, НЗ 50 TC (I) = T (I) 60 ЕСЛИ (NPLOP.EQ.5) GOTO 100 ДО 70 I = 1, НП 70 XOBS (I) = Т (I) DO 80 I = 1, НЗ 80 XCALC (I) = TC (I) GOTO 150 100 ДО 110 I = 1, НП 110 XOBS (I) = 1.D0 / T (I) DO 120 I = 1, НЗ 120 XCALC (I) = 1.D0 / TC (I) TEMP = XMIN XMIN = 1.D0 / XMAX XMAX = 1.D0 / TEMP 150 ВОЗВРАТ КОНЕЦ C Установка значений координаты Y для графика ПОДПРОГРАММА YDATA (T, XI, NP, TC, XIC, NC, NPLOP) НЕПРИЯТНОЕ РЕАЛЬНОЕ * 8 (A-H, O-Z) REAL * 4 XOBS, XCALC, YOBS, YCALC, XMIN, XMAX, YMIN, YMAX ОБЩИЙ / BLK7 / XOBS (300), XCALC (300), YOBS (300), YCALC (300) ОБЩИЙ / BLK8 / XMIN, XMAX, YMIN, YMAX РАЗМЕР T (1), XI (1), TC (1), XIC (1) GOTO (10,30,50,70,10), NPLOP 10 ДО 20 I = 1, НП 20 лет (I) = XI (I) DO 25 I = 1, НЗ 25 YCALC (I) = XIC (I) GOTO 100 30 Т7 = 7.997D0 ДО 40 I = 1, НП 40 лет (I) = DSQRT (T7 * XI (I) * T (I)) DO 45 I = 1, NC 45 YCALC (I) = DSQRT (T7 * XIC (I) * TC (I)) GOTO 100 50 ДО 60 I = 1, НП 60 лет (I) = XI (I) * T (I) DO 65 I = 1, НЗ 65 YCALC (I) = XIC (I) * TC (I) GOTO 100 70 ДО 80 I = 1, НП 80 лет (I) = 1.D0 / XI (I) DO 85 I = 1, НЗ 85 YCALC (I) = 1.D0 / XIC (I) 100 YMIN = YOBS (1) YMAX = YOBS (1) ДО 110 I = 2, НП ЕСЛИ (YMAX.LT.YOBS (I)) YMAX = YOBS (I) ЕСЛИ (YMIN.GT.YOBS (I)) YMIN = YOBS (I) 110 ПРОДОЛЖИТЬ DO 120 I = 1, НЗ ЕСЛИ (YMAX.LT.YCALC (I)) YMAX = YCALC (I) ЕСЛИ (YMIN.GT.YCALC (I)) YMIN = YCALC (I) 120 ПРОДОЛЖИТЬ ВОЗВРАЩЕНИЕ КОНЕЦ Процедура вызова C для графика графика ВСПОМОГАТЕЛЬНЫЙ СПЛОТ (NP, NC) ОБЩИЙ / BLK7 / XOBS (300), XCALC (300), YOBS (300), YCALC (300) ОБЩИЙ / BLK8 / XMIN, XMAX, YMIN, YMAX ОБЩИЙ / BLK5 / IN, IO XDIV = (XMAX-XMIN) / 5.YDIV = (YMAX-YMIN) / 5. ВЫЗОВ XYPLOT (XOBS, YOBS, NP, XCALC, YCALC, NC, 1, 1 XMIN, XMAX, XDIV, YMIN, YMAX, YDIV) ВОЗВРАЩЕНИЕ КОНЕЦ C Графическая программа построения графиков для совместимых с IBMPC с использованием C Microsoft Fortran V 5.0 ПОДПРОГРАММА XYPLOT (XO, YO, M, XC, YC, N, L1, 1 XMIN, XMAX, XDIV, YMIN, YMAX, YDIV) C XC — РАСЧЕТ КООРДИНАТЫ C YC — РАСЧЕТ Y-КООРДИНАТ C XO — ЭКСПЕРИМЕНТАЛЬНАЯ ЧАСТЬ X-КООРДИНАТ C YO — ЭКСПЕРИМЕНТАЛЬНАЯ ЧАСТЬ Y-КООРДИНАТ C N — КОЛИЧЕСТВО РАСЧЕТНЫХ ТОЧЕК C M — КОЛИЧЕСТВО ТОЧЕК ЭКСПТАЛЬНЫХ ДАННЫХ C L1, L2 — ТИПЫ ЛИНИЙ ДЛЯ КРИВЫХ YO и YC ТИПЫ ЛИНИИ C: 1 — Сплошная 2 — Пунктирная 3 — Пунктирная C 4 — пунктирная линия 5 — длинная тире 6 — двойная точка C (X0, Y0) ИСТОЧНИК.(XL, YL) РАЗМЕРЫ В КОДАХ ПЛОТТЕРА. C ОЖИДАЕТ, ЧТО ФАЙЛ ШРИФТА ROMAN.FON ИМЕЕТСЯ В ТЕКУЩЕМ КАТАЛОГЕ. РАЗМЕР XC (1), YC (1), XO (1), YO (1) C ————————————————- —— ВКЛЮЧИТЬ ‘FGRAPH.FD’ C ————————————————- ——- ЦЕЛОЕ * 2 FX, FY, DUMMY, LS (6), X1, Y1, X2, Y2, MX, MY, NMX, NMY СИМВОЛ STR * 5, STR1 * 2 ЗАПИСЬ / XYCOORD / S ЗАПИСЬ / ВИДЕОКОНФИГ / SCR FX (X) = (X-XMIN) * XS + X0 FY (Y) = NMY — ((Y-YMIN) * YS + Y0) ДАННЫЕ LS / # FFFF, #AAAA, #CCCC, #EBEB, # 7878, # 8888 / C Инициализировать графический режим ВЫЗОВ GETVIDEOCONFIG (SCR) НАД = SCR.АДАПТЕР C CGA ВЫСОКОЕ РАЗРЕШЕНИЕ BW ЕСЛИ (NAD.EQ. # 0002) РЕЖИМ = 6 C ГЕРКУЛЕС ЕСЛИ (NAD.EQ. # 0020) РЕЖИМ = 8 C EGA МОНО И ЦВЕТ ЕСЛИ ((NAD.EQ. # 0004 .OR. NAD.EQ. # 0044) 1 .И. SCR.MONITOR .EQ. # 0001) РЕЖИМ = 15 ЕСЛИ ((NAD.EQ. # 0004 .OR. NAD.EQ. # 0044) 1 .И. SCR.MONITOR .NE. # 0001) РЕЖИМ = 16 C VGA, MCGA ЕСЛИ (NAD.EQ. # 0008.OR.NAD.EQ. # 0048.OR.NAD.EQ. # 0010) MODE = 17 DUMMY = SETVIDEOMODE (INT2 (РЕЖИМ)) C Определите размеры графика ВЫЗОВ GETVIDEOCONFIG (SCR) NMX = SCR.NUMXPIXELS-1 NMY = SCR.NUMYPIXELS-1 X0 = NMX * 0,16 Y0 = NMY * .20 XL = NMX * .80 YL = NMY * 0,75 C Масштабные коэффициенты X и Y XS = XL / (XMAX-XMIN) YS = YL / (YMAX-YMIN) C Нарисуйте первую кривую ЕСЛИ (L1.LT.1 .OR. L1.GT.6) L1 = 1 ВЫЗОВ УСТАНОВИТЬ СТИЛЬ (LS (L1)) ВЫЗОВ MOVETO (FX (XC (1)), FY (YC (1)), S) ДО 20 I = 2, N 20 DUMMY = LINETO (FX (XC (I)), FY (YC (I))) C Нарисуйте вторую кривую (отметьте x) XCR = 2 ДО 30 I = 1, М ВЫЗОВ MOVETO (FX (XO (I)) — XCR, FY (YO (I)) — XCR, S) DUMMY = LINETO (FX (XO (I)) + XCR, FY (YO (I)) + XCR) ВЫЗОВ MOVETO (FX (XO (I)) — XCR, FY (YO (I)) + XCR, S) 30 DUMMY = LINETO (FX (XO (I)) + XCR, FY (YO (I)) — XCR) C Нарисуйте рамку NH = 12 NW = 10 DUMMY = REGISTERFONTS (‘*.фон ‘) DUMMY = SETFONT («t’ROMAN ‘h22 w10″) MX = 5 * NW / 2 MY = NH / 2 СТИЛЬ ВЫЗОВА 40 (LS (1)) ВЫЗОВ MOVETO (FX (XMIN), FY (YMIN), S) DUMMY = LINETO (FX (XMIN), FY (YMAX)) DUMMY = LINETO (FX (XMAX), FY (YMAX)) DUMMY = LINETO (FX (XMAX), FY (YMIN)) DUMMY = LINETO (FX (XMIN), FY (YMIN)) C Нарисуйте отметки и значения по оси X L = ALOG10 (XMAX-XMIN) ST = XMIN Y1 = NMY- (Y0 + NMY * 0,025) Y2 = NMY- (Y0-NMY * 0,025) 50 ЗВОНОК MOVETO (FX (ST), FY (YMIN), S) DUMMY = LINETO (FX (ST), Y1) А = ST / 10.0 ** л ЗАПИСАТЬ (STR, ‘(F5.2)’) A X2 = FX (ST) -MX ЕСЛИ (X2.GT.INT2 (NMX * .95)) X2 = NMX * .95 ВЫЗОВ MOVETO (X2, Y2, S) ВЫЗОВ OUTGTEXT (STR) ST = ST + XDIV ЕСЛИ (ST.LE.XMAX) ПЕРЕЙТИ 50 ЕСЛИ (L.EQ.0) GOTO 55 ВЫЗОВ MOVETO (INT2 (X0 + XL-10 * NW), INT2 (Y0 + YL + 2,5 * NH), S) ВЫЗОВ OUTGTEXT (‘X10’) ВЫЗОВ MOVETO (INT2 (X0 + XL-7 * NW), INT2 (Y0 + YL + 2.0 * NH), S) ЕСЛИ (L.GT.0 .AND. L.LE.9) ЗАПИШИТЕ (STR1, ‘(I1)’) L ЕСЛИ (L.LT.0 .OR. L.GT.9) ЗАПИСАТЬ (STR1, ‘(I2)’) L 54 ТЕКСТ ВЫЗОВА (STR1) C Нарисуйте отметки и значения на оси Y 55 L = ALOG10 (YMAX-YMIN) ST = YMIN X1 = X0 + 0,008 * NMX X2 = X0-5,5 * NW 60 ЗВОНОК MOVETO (FX (XMIN), FY (ST), S) DUMMY = LINETO (X1; FY (ST)) A = ST / 10,0 ** L ЗАПИСАТЬ (STR, ‘(F5.2)’) A Y2 = FY (ST) -MY ЕСЛИ (Y2.LT.INT2 (0)) Y2 = 0 ВЫЗОВ MOVETO (X2, Y2, S) ВЫЗОВ OUTGTEXT (STR) ST = ST + YDIV ЕСЛИ (ST.LE.YMAX) ПЕРЕЙТИ 60 ЕСЛИ (Л.EQ.0) GOTO 70 ВЫЗОВ MOVETO (INT2 (X0-9 * NW), INT2 (NMY- (Y0 + YL-1.5 * NH)), S) ВЫЗОВ OUTGTEXT (‘X10’) ЕСЛИ (L.GT.9 .OR. L.LT.0) ЗАПИШИТЕ (STR1, ‘(I2)’) L ЕСЛИ (L.GT.0 .AND. L.LT.9) ЗАПИШИТЕ (STR1, ‘(I1)’) L ВЫЗОВ MOVETO (INT2 (X0-6 * NW), INT2 (NMY- (Y0 + YL-NH)), S) ВЫЗОВ OUTGTEXT (STR1) C Дождитесь возврата инструкции пользователя 70 ПРОЧИТАТЬ (*, *) ЗВОНИТЕ UNREGISTERFONTS () DUMMY = SETVIDEOMODE (INT2 (-1)) ВОЗВРАЩЕНИЕ КОНЕЦ C Программа печати на принтере ВСПОМОГАТЕЛЬНЫЙ PPLOT (NP) ХАРАКТЕР * 1 A (75), ПУСТО, ЗВЕЗДА, РАСЧЕТ ОБЩИЙ / BLK7 / X (300), XCALC (300), Y (300), Z (300) ОБЩИЙ / BLK8 / XMIN, XMAX, YMIN, YMAX ОБЩИЙ / BLK5 / IN, IO ДАННЫЕ ПУСТО, ЗВЕЗДА, РАСЧЕТ / », ‘*’, ‘c’ / NX = 74 NY = 20 XS = NX / (XMAX-XMIN) YS = NY / (YMAX-YMIN) YINT = 1.D0 / YS T2 = YMAX Т1 = Т2 NX1 = NX + 1 115 Т1 = Т1-ИНТ DO 120 Дж = 1, NX1 120 А (Дж) = ПУСТО NX1 = 0 DO 140 I = 1, НП ЕСЛИ (Y (I) .LT.T1 .OR. Y (I) .GT. T2) GOTO 130 J = (X (I) -XMIN) * XS + 1,5 ЕСЛИ (J.GT.NX1) NX1 = J A (J) = ЗВЕЗДА 130 IF (Z (I) .LT. T1 .OR. Z (I) .GT. T2) GOTO 140 J = (X (I) -XMIN) * XS + 1,5 ЕСЛИ (J.GT.NX1) NX1 = J A (J) = РАСЧЕТ 140 ПРОДОЛЖИТЬ ЗАПИСАТЬ (IO, 150) (A (J), J = 1, NX1) 150 ФОРМАТ (4X, 75A1) Т2 = Т1 ЕСЛИ (T1.GT. YMIN) GOTO 115 ПРОЧИТАЙТЕ (ВХОД, 160) A (1) 160 ФОРМАТ (A) 200 ВОЗВРАТ КОНЕЦ C Создает файл SUSCEP.HP для печати на плоттере HP7475A. ПОДПРОГРАММА HPPLOT (NP, NC, NPLOP) ЦЕЛОЕ X0, Y0, XX, YY, XSHFT, YSHFT РЕАЛЬНЫЙ ДЮЙМ ПЕРСОНАЖ Ch4 * 1 ОБЩИЙ / BLK7 / XOBS (300), XCALC (300), YOBS (300), YCALC (300) ОБЩИЙ / BLK8 / XMIN, XMAX, YMIN, YMAX ОБЩИЙ / BLK5 / IN, IO IX (X) = (X-XMN) * XS + X0 IY (Y) = (Y-YMN) * YS + Y0 Ch4 = СИМВОЛ (3) ДЮЙМ = 402 * 2.54 X0 = 2,0 * ДЮЙМ Y0 = 1,5 * ДЮЙМ XL = 8,5 * ДЮЙМ YL = 6,0 * ДЮЙМ XD = (XMAX-XMIN) / 10. YD = (YMAX-YMIN) / 10. WRITE (IO, *) ‘Значения из данных:’ WRITE (IO, 5) ‘X-min =’, XMIN, ‘X-max =’, XMAX, ‘X-divn =’, XD 5 ФОРМАТ (3 (A, F10.4)) WRITE (IO, *) ‘Введите значения для графика’ ЧИТАТЬ (IN, *) XMN, XMX, XD WRITE (IO, 5) ‘Y-min =’, YMIN, ‘Y-max =’, YMAX, ‘Y-divn =’, YD WRITE (IO, *) ‘Введите значения для графика’ ЧИТАТЬ (IN, *) YMN, YMX, YD XS = XL / (XMX-XMN) YS = YL / (YMX-YMN) C Нарисуйте кривую ОТКРЫТЬ (UNIT = 3, FILE = ‘SUSCEP.HP ‘, СТАТУС =’ НЕИЗВЕСТНО ‘) WRITE (IO, *) ‘Имя файла HPPLOT — SUSCEP.HP’ ЗАПИСАТЬ (3,10) C Инициализировать плоттер 10 FORMAT (1X, ‘IN; SP1; PU;’) C Постройте расчетную кривую DO 20 I = 1, НЗ 20 ЗАПИСАТЬ (3,21) IX (XCALC (I)), IY (YCALC (I)) 21 ФОРМАТ (1X, ‘PA’, I5, ‘,’, I5, ‘PD’) ЗАПИСАТЬ (3,22) 22 FORMAT (1X, ‘PU;’) C Постройте треугольники для экспериментальных точек DO 25 I = 1, NP 25 ЗАПИСЬ (3,26) IX (XOBS (I)), IY (YOBS (I)) 26 FORMAT (1X, ‘PA’, I5, ‘,’, I5, ‘CI60,10WG60,0,360,10;’) C Нарисуйте рамку ЗАПИСАТЬ (3,21) IX (XMN), IY (YMN) ЗАПИСАТЬ (3,21) IX (XMX), IY (YMN) ЗАПИСАТЬ (3,21) IX (XMX), IY (YMX) ЗАПИСАТЬ (3,21) IX (XMN), IY (YMX) ЗАПИСАТЬ (3,21) IX (XMN), IY (YMN) ЗАПИСАТЬ (3,22) C Запись X-меток ЗАПИСАТЬ (3,27) 27 FORMAT (1X, ‘SI0.2,0,3; ‘) XTB = ЦЕЛОЕ (XMN / XD) * XD ЕСЛИ (XTB .LT. XMN) XTB = XTB + XD NTX = (XMX-XTB) / XD +1 DO 50 I = 1, NTX XX = (XTB-XMN) * XS + X0 ЗАПИСАТЬ (3,30) XX, Y0 30 FORMAT (1X, ‘PA’, I5, ‘,’, I5, ‘XT;’) IXTB = XTB + 0,5 ЕСЛИ (NPLOP.EQ.5) GOTO 40 ЕСЛИ (IXTB.LT.10) ЗАПИСАТЬ (3,31) IXTB, Ch4 ЕСЛИ (IXTB.GE.10 .AND. IXTB.LT.100) ЗАПИСАТЬ (3,32) IXTB, Ch4 ЕСЛИ (IXTB.GE.100) ЗАПИСАТЬ (3,33) IXTB, Ch4 31 ФОРМАТ (1X, ‘CP-.5, -1,5; LB ‘, I1, A) 32 ФОРМАТ (1X, ‘CP-1., — 1.5; LB’, I2, A) 33 ФОРМАТ (1X, ‘CP-1.5, -1,5; LB’, I3, A) GOTO 50 40 ЗАПИСЬ (3,41) XTB, Ch4 41 ФОРМАТ (1X, ‘CP-1.5, -1.5; LB’, F3.2, A) 50 XTB = XTB + XD C Напишите Y-галочки YTB ​​= INT (YMN / YD) * YD ЕСЛИ (YTB.LT.YMN) YTB = YTB + YD NTY = (YMX-YTB) / YD + 1 DO 100 I = 1, NTY ГГ = (YTB-YMN) * YS + Y0 ЗАПИСАТЬ (3,60) X0, YY 60 FORMAT (1X, ‘PA’, I5, ‘,’, I5, ‘YT;’) GOTO (80,70,70,90,80), NPLOP 70 ЗАПИСЬ (3,71) YTB, Ch4 71 ФОРМАТ (1X, ‘CP-3.5, -. 5, LB ‘, F3.1, A) GOTO 100 80 IYTB = YTB * 1000 ЗАПИСАТЬ (3,81) IYTB, Ch4 81 ФОРМАТ (1X, ‘CP-3.5, -. 5, LB’, I3, A) GOTO 100 90 IYEXP = 0 ЕСЛИ (YMX.LE.0) GOTO 91 IYEXP = LOG10 (YMX) 91 YEXP = 10. ** IYEXP YTBE = YTB / YEXP ЗАПИСАТЬ (3,71) YTBE, Ch4 100 YTB = YTB + YD C Название оси X XSHFT = X0 + XL / 2. ЗАПИСАТЬ (3,110) XSHFT, Y0 110 FORMAT (1X, ‘PA’, I5, ‘,’, I5, ‘;’) ЕСЛИ (NPLOP.LT.5) ЗАПИСАТЬ (3,111) Ch4 ЕСЛИ (NPLOP.EQ.5) ЗАПИСАТЬ (3,112) Ch4, Ch4, Ch4 111 ФОРМАТ (1X, ‘CP-3, -3; LBT (K)’, A) 112 ФОРМАТ (1X, ‘CP-4, -3; LB1 / T (K’, A, ‘CP0, .5; SI.18, .25; LB-1’, A, 1 ‘SI.2, .3; CP0, -. 5; LB) ‘, A) НАПИСАТЬ (3,195) C название оси Y YSHFT = Y0 + YL / 2 ЗАПИСАТЬ (3,130) X0, YSHFT 130 FORMAT (1X, ‘PA’, I5, ‘,’, I5, ‘SI.2, .3; DI0,1;’) GOTO (140,150,160,170,140), NPLOP 140 ЗАПИСАТЬ (3,141) 141 ФОРМАТ (1X, ‘CP-6, -3;’) ЗАПИСАТЬ (3,180) ЗАПИСАТЬ (3,142) Ch4, Ch4, Ch4, Ch4, Ch4, Ch4, Ch4 142 FORMAT (1X, ‘LB (x10’, A, ‘CP0 ,.5; SI.18, .25; LB3 ‘, A, 1 ‘SI.2, .3; CP0, -. 5; LBcm’, A, ‘CP0, .5; SI.18, .25; LB3’, A, / 2 ‘SI.2, .3; CP0, -. 5; LBmol’, A, ‘CP0, .5; SI.18, .25; LB-1’, 3 A, ‘SI.2, .3; CP0, -. 5LB)’, A) GOTO 200 150 ЗАПИСАТЬ (3,151) 151 ФОРМАТ (1X, ‘CP-4, -3;’) C MU ЗАПИСАТЬ (3,185) ЗАПИСАТЬ (3,152) Ch4 152 ФОРМАТ (1X, ‘LB (B.M.)’, A) GOTO 200 160 ЗАПИСАТЬ (3,161) 161 ФОРМАТ (1X, ‘CP-10, -3;’) ЗАПИСАТЬ (3,180) ЗАПИСАТЬ (3162) Ch4, Ch4, Ch4, Ch4, Ch4, Ch4, Ch4 162 FORMAT (1X, ‘LB T (x10’, A, ‘CP0 ,.5; SI.18, .25; LB3 ‘, A, 1 ‘SI.2, .3; CP0, -. 5; LBcm’, A, ‘CP0, .5;’ / ‘SI.18, .25; LB3’, A, 2 ‘SI.2, .3; CP0, -. 5;’ / ‘LBmol’, A, ‘CP0, .5; SI.18, .25; LB-1’, 3 A, ‘SI.2, .3; CP0, -. 5LBK)’, A) GOTO 200 170 ЗАПИСАТЬ (3,171) 171 ФОРМАТ (1X, ‘CP-8, -3;’) ЗАПИСАТЬ (3,180) ЗАПИСАТЬ (3,172) Ch4, IYEXP, Ch4, Ch4, Ch4, Ch4 172 ФОРМАТ (1X, ‘CP-3,0LB1 / (x10’, A, ‘CP0, .5; SI.18, .25; LB’, I1, A, 1 ‘SI.2, .3; CP0, -. 5; LB см’, A, ‘CP0, .5;’ / ‘SI.18, .25; LB-3’, A, 2 ‘SI.2, .3; CP0, -. 5; LBmol) ‘, A) GOTO 200 C Написать Чи ФОРМАТ 180 (1X, ‘PR-1206,0;’ / 1 ‘ПР-213,24ПД; ПР-1,2ПР1,4ПР4,4ПР6, -2ПР3, -4ПР1, -6;’, / 1 ‘ПР-1, -6ПР-4, -4ПР-7,0ПР-9,4ПР-4,8ПР-1,12ПР1,12ПР4,12;’, / 2 ‘ПР9,16ПР15,16ПР78,44ПР4,8ПР0,4ПР-1,4ПР-3, -8ПР-3,0;’, / 3 ‘ПР-6,4ПР-1,4ПР-1,8ПР1,8ПР7,4ПР6, -2ПР6, -6ПР1, -8;’, / 4 ‘ПР1, -12ПР0, -8ПР-3, -12ПР-87, -52ПР-12, -12ПР-9, -16;’, / 5 ‘ПР-3, -8ПР0, -8ПР3, -4ПР4,0ПУ;’, 6 ‘ПР103,4ПД; ПР-9, -12ПР-108,144ПР9,12ПР108, -144ПУ;’) C Написать мю 185 ФОРМАТ (1X, ‘PR-1206,0;’ / 1 ‘ПР0,0ПД; ПР-120,40,0,20,42, -8,12,4,4,4,1,12;’, / 1 ‘ПР-1,8, -4,8, -7,8, -46,24,0,20,51, -20,4,2,1,6;’ / 2 ‘ПР-1,8, -4,8,4,0,4, -2,4, -6,1, -8, -1, -8, -4, -8;’, / 3 ‘ПР-6, -6,6, -6,7, -12,1, -16,0, -8, -3, -12, -4, -4;’, / 4 ‘ПР-4, -4, -24,6,90, -40,0, -20ПУ;’) C Стрелка вверх ФОРМАТ 190 (1X, ‘PR-1206,0;’, ‘PR-51,0PD; PR0,120PR18, -20PR-27,60’ / 1 ‘ПР-27, -60ПР18,20ПР0, -120ПР18,0ПУ;’) C Стрелка вправо 195 ФОРМАТ (1X, ‘PR100,50PD; PR90,0PR-15, -24PR45,36;’ / 1 ‘ПР-45,36ПР15, -24ПР-90,0ПР0, -24ПУ;’) 200 НАПИСАТЬ (3,210) 210 ФОРМАТ (1X, ‘PU; SP0;’) ЗАКРЫТЬ (UNIT = 3, STATUS = ‘KEEP’) ВОЗВРАЩЕНИЕ КОНЕЦ C C Ссылка: Оценка ошибок с использованием последовательного симплекс-метода в C Нелинейный анализ данных методом наименьших квадратов C Грегори Р.Филлипс и Эдвард М. Айринг C Анал. Chem. 1988, том 60, стр 738-741 C и ссылки в нем ВСПОМОГАТЕЛЬНАЯ ОШИБКА (XMIN, SIGMA, N, NP, IER) C C XMIN — это конвергентный параметр, установленный симплексным методом C P конечный симплекс; N — Нет параметров C NP — Количество наблюдаемых точек C IER = 1, если процедура терпит неудачу НЕПРИЯТНОЕ РЕАЛЬНОЕ * 8 (A-H, O-Z) ДВОЙНАЯ ТОЧНОСТЬ XMIN (12), MIN (12), SIGMA (12) РАЗМЕР B (12,12), BI (12,12), Q (12,12), A (12) 1, PBAR (12), Y (13), Y0I (12), BIA (12), BBA (12) 2, PP (20,21), YY (13) ОБЩИЙ / ERRCAL / P (20,21) ДАННЫЕ НУЛЬ, ПОЛОВИНА, ДВА, THR, ЧЕТЫРЕ / 0.0D0,0.5D0,2.0D0,3.0D0,4.0D0 / C DO 10 I = 1, N 10 СИГМА (I) = НУЛЬ ЕСЛИ (NP.LE.N) GOTO 550 IER = 0 КРУГЛЫЙ = 1.D-12 NN = N + 1 NM1 = N-1 DNN = FLOAT (NN) C DO 30 I = 1, NN DO 20 Дж = 1, Н 20 МИН (Дж) = P (Дж, I) 30 Y (I) = FN (МИН) C 100 КРУГЛЫЙ = КРУГЛЫЙ * 10.D0 ЕСЛИ (ROUND.GT.1.0D-5) GOTO 550 ИТЭР = 0 IEXPN = 1 C РАССЧИТАТЬ V0 110 МОТ = 1 YLO = Y (1) DO 200 I = 1, NN ЕСЛИ (Y (I).GE.YLO) GOTO 200 YLO = Y (I) МОТ = I 200 ПРОДОЛЖИТЬ ЕСЛИ (IEXPN.EQ.0) GOTO 310 C ПРОВЕРЬТЕ, НЕОБХОДИМО РАСШИРЕНИЕ ДО 220 I = 1, NN ЕСЛИ (I.EQ.ILO) GOTO 220 ЕСЛИ ((Y (I) -YLO) / YLO .LT. ROUND) GOTO 230 220 ПРОДОЛЖИТЬ GOTO 280 РАСШИРЕНИЕ C 230 ДО 250 I = 1, Н Z = НУЛЬ DO 240 Дж = 1, НН 240 Z = Z + P (I, J) 250 ПБАР (I) = Z / NN DO 270 I = 1, NN DO 260 Дж = 1, Н P (J, I) = PBAR (J) + ДВА * (P (J, I) -PBAR (J)) 260 МИН (Дж) = P (Дж, I) 270 Y (I) = FN (МИН) GOTO 110 C СОХРАНИТЬ СИМПЛЕКС 280 ДО 290 I = 1, Н YY (I) = Y (I) DO 290 Дж = 1, Н 290 PP (J, I) = P (J, I) РАСЧЕТ НАЧИНАЕТСЯ 300 IEXPN = 0 GOTO 110 C ПЕРЕМЕСТИТЕ КОЛОННУ ILO НА КОЛОННУ NN. 310 DO 320 Дж = 1, Н 320 МИН (Дж) = P (Дж, МОТ) Y0 = Y (МОТ) DO 330 I = МОТ, N I1 = I + 1 Y (I) = Y (I1) DO 330 Дж = 1, Н 330 P (J, I) = P (J, I1) DO 340 Дж = 1, Н 340 P (J, NN) = MIN (Дж) Y (NN) = Y0 C СОЗДАЙТЕ МАТРИЦЫ A, B, Q DO 360 I = 1, N DO 350 Дж = 1, Н Q (J, I) = (P (J, I) -P (J, NN)) 350 МИН (Дж) = (P (J, I) + P (J, NN)) * ПОЛОВИНА Y0I (I) = FN (МИН) A (I) = ЧЕТЫРЕ * Y0I (I) -Y (I) -THR * Y0 360 Б (I, I) = ДВА * (Y (I) + Y0-ДВА * Y0I (I)) DO 380 I = 1, NM1 I1 = I + 1 DO 380 Дж = I1, Н DO 370 K = 1, N 370 МИН (K) = (P (K, I) + P (K, J)) * ПОЛОВИНА YIJ = FN (MIN) B (I, J) = ДВА * (YIJ + Y0-Y0I (I) -Y0I (J)) 380 Б (J, I) = B (I, J) ВЫЗОВ ИНВЕРТИРУЙТЕ (B, BI, N, EPS) ЕСЛИ (EPS.LT.ZERO) GOTO 100 DO 400 I = 1, N BIA (I) = НУЛЬ ДО 400 К = 1, Н 400 BIA (I) = BIA (I) + BI (I, K) * A (K) DO 410 I = 1, N BBA (I) = НУЛЬ DO 410 К = 1, Н 410 BBA (I) = BBA (I) + BI (I, K) * BIA (K) АББА = НОЛЬ ДО 420 К = 1, Н 420 ABBA = ABBA + A (K) * BBA (K) ЕСЛИ (ABBA.GT. 0.25D0) ИТЭР = ИТЭР + 1 ЕСЛИ (ITER.LE.4) GOTO 450 C СЛИШКОМ МНОГО ИТЕРАЦИЙ. РАСШИРИТЬ SIMPLEX. DO 430 I = 1, NN Y (I) = YY (I) DO 430 Дж = 1, Н 430 P (J, I) = PP (J, I) GOTO 100 C ПРАВИЛЬНО THETA0 450 ДО 470 Дж = 1, Н QBAJ = НУЛЬ DO 460 К = 1, Н 460 QBAJ = QBAJ + Q (J, K) * BIA (K) P (J, NN) = P (J, NN) -QBAJ * HALF 470 МИН (J) = P (J, NN) Y (NN) = FN (MIN) ЕСЛИ (ABBA.GT. 0.25D0) GOTO 300 C РАСЧЕТ СИГМЫ 500 S2 = Y (NN) / (NP-N) DO 520 I = 1, Н DO 520 Дж = 1, Н B (I, J) = НУЛЬ DO 520 К = 1, Н 520 Б (I, J) = B (I, J) + BI (I, K) * Q (J, K) DO 540 I = 1, Н СИГМА (I) = НОЛЬ DO 530 К = 1, Н 530 СИГМА (I) = СИГМА (I) + Q (I, K) * B (K, I) ЕСЛИ (SIGMA (I) .LT.ZERO) GOTO 100 540 СИГМА (I) = DSQRT (СИГМА (I) * S2) DO 545 I = 1, N 545 XMIN (I) = MIN (I) ВОЗВРАЩЕНИЕ 550 IER = 1 ВОЗВРАЩЕНИЕ КОНЕЦ C Вызывает подпрограмму SSP GJRD для инверсии.C Критерий сходимости = (Наибольший элемент B) * 1.e-6 ВСПОМОГАТЕЛЬНАЯ ИНВЕРТАЦИЯ (B, BI, N, EPS) НЕПРИЯТНОЕ РЕАЛЬНОЕ * 8 (A-H, O-Z) РАЗМЕР B (12,12), BI (12,12), ТЕМП (144) IO = 6 EPS = 0.D0 К = 0 DO 5 J = 1, N ДО 5 I = 1, N ЕСЛИ (EPS.LT.DABS (B (I, J))) EPS = DABS (B (I, J)) К = К + 1 5 ТЕМП (K) = B (I, J) EPS = EPS * 1.D-6 ПОЗВОНИТЬ GJRD (TEMP, N, EPS) C IF (EPS.EQ.-1.0D0) WRITE (IO, *) ‘Ошибка инверсии матрицы’ К = 0 DO 10 Дж = 1, Н DO 10 I = 1, N К = К + 1 10 BI (I, J) = ТЕМП (K) ВОЗВРАЩЕНИЕ КОНЕЦ C IBM 1130 — НАУЧНЫЙ ПОДПРОГРАММНЫЙ ПАКЕТ ВСПОМОГАТЕЛЬНЫЙ GJRD (A, N, EPS) C Следующие строки изменены для двойной точности C все операторы печати подавлены C НЕПРИЯТНОЕ РЕАЛЬНОЕ * 8 (A-H, O-Z) РАЗМЕР A (1), L (50), M (50) ABS (X) = DABS (X) C ————————————————- ————— C N — ПОРЯДОК МАТРИЦЫ A (I, K) ОБРАТИТЬСЯ, EPS ЯВЛЯЕТСЯ A C ДОПУСК ПРИНЯТИЯ СИНГУЛЯРНОСТИ ДАННОЙ МАТРИЦЫ C ЕСЛИ МАТРИЦА ОСОБЕННО, ТО EPS = -1 C ОПРЕДЕЛЕНИЕ ПОВОРОТНОГО ЭЛЕМЕНТА C IPR = 6 ПЧ (Н-50) 1,1,99 1 ПРОДОЛЖИТЬ C ПОИСК КРУПНЕЙШИХ ЭЛЕМЕНТОВ D = 1.0 NK = -N ДО 80 К = 1, Н НК = НК + Н L (К) = К М (К) = К КК = НК + К БИГА = А (КК) DO 20 J = K, N IZ = N * (J-1) ДО 20 I = K, N IJ = IZ + I 10 IF (АБС (BIGA) — АБС ​​(A (IJ))) 15,20,20 15 BIGA = A (IJ) L (K) = I M (K) = J 20 ПРОДОЛЖИТЬ C ОБМЕН РЯДАМИ J = L (К) ЕСЛИ (JK) 35,35,25 25 КИ = К-Н DO 30 I = 1, N KI = KI + N УДЕРЖАТЬ = -A (KI) JI = KI-K + J A (KI) = A (JI) 30 А (JI) = УДЕРЖАТЬ C ОБМЕННЫЕ КОЛОННЫ 35 I = M (К) IF (I-K) 45,45,38 38 JP = N * (I-1) DO 40 Дж = 1, Н JK = NK + J JI = JP + J УДЕРЖАТЬ = -A (JK) A (JK) = A (JI) 40 А (JI) = УДЕРЖАТЬ C РАЗДЕЛИТЬ КОЛОНКУ НА МИНУС ПОВОРОТА (ЗНАЧЕНИЕ ПОВОРОТНОГО ЭЛЕМЕНТА C СОДЕРЖИТСЯ В БИГЕ) 45 IF (ABS (BIGA) — EPS) 46,46,48 46 EPS = — 1.0 C НАПИСАТЬ (IPR, 47) K, BIGA 47 FORMAT (24H0MATRIX SINGULAR * STAGE, I6,4X5HPIVOTE20.10) ВОЗВРАЩЕНИЕ 48 ДО 55 I = 1, N IF (I-K) 50,55,50 50 ИК = НК + Я А (ИК) = А (ИК) / (- БИГА) 55 ПРОДОЛЖИТЬ C УМЕНЬШИТЬ МАТРИЦУ DO 65 I = 1, N IK = NK + I УДЕРЖАТЬ = A (IK) IJ = I-N DO 65 Дж = 1, Н IJ = IJ + N IF (I-K) 60,65,60 60 IF (J-K) 62,65,62 62 кДж = IJ-I + K A (IJ) = ДЕРЖАТЬ * A (KJ) + A (IJ) 65 ПРОДОЛЖИТЬ C РАЗДЕЛИТЬ РЯД ПО ПОВОРОТУ КДж = К-Н DO 75 Дж = 1, Н КДж = КДж + N ЕСЛИ (J-K) 70,75,70 70 А (кДж) = A (кДж) / BIGA 75 ПРОДОЛЖИТЬ C ПРОДУКТ ПОВОРОТОВ D = D * BIGA C ЗАМЕНИТЕ ПОВОРОТ НА ВЗАИМНО. А (КК) = 1.0 / БИГА 80 ПРОДОЛЖИТЬ C ОБМЕН ЗАКЛЮЧИТЕЛЬНЫМ РЯДОМ И КОЛОНКОЙ K = N 100 К = (К-1) IF (К) 150,150,105 105 I = L (К) IF (I-K) 120,120,108 108 JQ = N * (К-1) JR = N * (I-1) DO 110 Дж = 1, Н JK = JQ + J УДЕРЖАТЬ = A (JK) JI = JR + J A (JK) = — A (JI) 110 А (JI) = УДЕРЖАТЬ 120 Дж = М (К) ЕСЛИ (JK) 100,100,125 125 КИ = К-Н DO 130 I = 1, N KI = KI + N УДЕРЖАТЬ = A (KI) JI = KI-K + J A (KI) = — A (JI) 130 А (JI) = УДЕРЖАТЬ ПЕРЕЙТИ К 100 150 ВОЗВРАТ 99 ПРОДОЛЖИТЬ C ЗАПИСАТЬ (IPR, 101) N 101 ФОРМАТ (39H0 ***** ПОДПРОГРАММА GJRA ‘N = 50 N =, I10, /) ОСТАНАВЛИВАТЬСЯ КОНЕЦ C C АЛГОРИТМ C КАК 47 ПРИКЛАДНАЯ СТАТИСТИКА (J.R.STATIST.SOC C), C (1971) ТОМ 20, № 3 C ПОДПРОГРАММА NELMIN (N, START, MIN, YNEWLO, REQMIN, STEP, KONGVE, ICOUNT) ДВОЙНОЙ ТОЧНЫЙ СТАРТ (N), МИН (N), YNEWLO, REQMIN, STEP (N), 1 P (20,21), PSTAR (20), P2STAR (20), PBAR (20), Y (20), 2 DN, DNN, Z, SUM, SUMM, YLO, RCOEFF, YSTAR, ECOEFF, Y2STAR, CCOEFF, 3 CURMIN, DEL, FN ОБЩИЙ / ERRCAL / P C ДАННЫЕ RCOEFF / 1.D0 /, ECOEFF / 2.D0 /, CCOEFF / 5.D-1 / C КОЭФФИЦИЕНТЫ ОТРАЖЕНИЯ, РАСШИРЕНИЯ И СЖАТИЯ.C C ПРОВЕРКА ДЕЙСТВИЯ ВХОДНЫХ ПАРАМЕТРОВ. KCOUNT = ICOUNT ICOUNT = 0 ЕСЛИ (REQMIN.LE.0.D0) ICOUNT = ICOUNT-1 ЕСЛИ (N.GT.20) ICOUNT = ICOUNT-10 ЕСЛИ (KONGVE.LE.0) ICOUNT = ICOUNT-100 ЕСЛИ (ICOUNT.LT.0) ВОЗВРАТ C JCOUNT = КОНГВЕ DN = FLOAT (N) NN = N + 1 DNN = FLOAT (NN) DEL = 1.D0 C C ПОСТРОЕНИЕ НАЧАЛЬНОГО СИМПЛЕКСА C 1001 ДО 1 I = 1, N 1 P (I, NN) = НАЧАЛО (I) Z = FN (СТАРТ) Y (NN) = Z СУММ = Z СУММ = Z * Z DO 2 J = 1, N СТАРТ (J) = СТАРТ (J) + ШАГ (J) * УДАЛ. ДО 3 I = 1, N 3 P (I, J) = СТАРТ (I) Z = FN (СТАРТ) Y (J) = Z СУММ = СУММ + Z СУММ = СУММ + Z * Z 2 СТАРТ (J) = СТАРТ (J) -ШАГ (J) * УДАЛ. C C SIMPLEX КОНСТРУКЦИЯ ЗАВЕРШЕНА C C НАЙДИТЕ НАИБОЛЬШИЕ И НАИБОЛЬШИЕ ЗНАЧЕНИЯ Y.YNEWLO (= Y (IHI)) УКАЗЫВАЕТ C ВЕРТЕКС СИМПЛЕКСА ПОДЛЕЖИТ ЗАМЕНИТЕ C 1000 YLO = Y (1) YNEWLO = YLO МОТ = 1 IHI = 1 DO 5 I = 2, NN ЕСЛИ (Y (I) .GE.YLO) ПЕРЕХОДИТЕ К ЭТАПУ 4 YLO = Y (I) МОТ = I 4 ЕСЛИ (Y (I) .LE.YNEWLO) ПЕРЕХОДИТЕ К ЭТАПУ 5 YNEWLO = Y (I) IHI = I 5 ПРОДОЛЖИТЬ СУММ = СУММ-YNEWLO SUMM = SUMM-YNEWLO * YNEWLO C C ВЫЧИСЛЕНИЕ PBAR, ЦЕНТРОИДА ВЕРТИФИКАЦИИ SIMPLEX C ИСКЛЮЧАЯ ЗНАЧЕНИЕ Y YNEWLO DO 7 I = 1, N Z = 0.D0 DO 6 J = 1, NN 6 Z = Z + P (I, J) Z = Z-P (I, IHI) 7 PBAR (I) = Z / DN C C ОТРАЖЕНИЕ ЧЕРЕЗ ЦЕНТРОИД C ДО 8 I = 1, N 8 PSTAR (I) = (1.D0 + RCOEFF) * PBAR (I) -RCOEFF * P (I, IHI) YSTAR = FN (PSTAR) ICOUNT = ICOUNT + 1 ЕСЛИ (YSTAR.GE.YLO) ПЕРЕЙДИТЕ К 12 C C УСПЕШНЫХ ОТРАЖЕНИЙ, ТАК РАСШИРЕНИЕ C ДО 9 I = 1, N 9 P2STAR (I) = ECOEFF * PSTAR (I) + (1.D0-ECOEFF) * PBAR (I) Y2STAR = FN (P2STAR) ICOUNT = ICOUNT + 1 C C УДЕРЖИВАЙТЕ РАСТЯНИЕ ИЛИ СЖАТИЕ C ЕСЛИ (Y2STAR.GE.YLO) ПЕРЕЙДИТЕ К ЭТАПУ 19 10 ДО 11 I = 1, N 11 P (I, IHI) = P2STAR (I) Y (IHI) = Y2STAR ПЕРЕЙТИ К 900 C C НЕТ РАСШИРЕНИЯ C 12 L = 0 DO 13 I = 1, NN ЕСЛИ (Y (I) .GT.YSTAR) L = L + 1 13 ПРОДОЛЖИТЬ ЕСЛИ (L.GT.1) ПЕРЕЙДИТЕ К 19 ЕСЛИ (L.EQ.0) ПЕРЕЙДИТЕ К 15 C C СОКРАЩЕНИЕ ИЛИ ОТРАЖЕНИЕ ЦЕНТРОИДА C DO 14 I = 1, N 14 P (I, IHI) = PSTAR (I) Y (IHI) = YSTAR C C КОНТРАКЦИЯ НА СТОРОНЕ Y (IHI) ЦЕНТРОИДА C 15 ДО 16 I = 1, N 16 P2STAR (I) = CCOEFF * P (I, IHI) + (1.D0-CCOEFF) * PBAR (I) Y2STAR = FN (P2STAR) ICOUNT = ICOUNT + 1 ЕСЛИ (Y2STAR.LE.Y (IHI)) ПЕРЕЙДИТЕ К 10 C C КОНТРАКТ ВСЕГО SIMPLEX C СУММ = 0.D0 СУММ = 0.D0 DO 18 J = 1, NN ДО 17 I = 1, N P (I, J) = (P (I, J) + P (I, ILO)) * 0.5D0 17 МИН (I) = P (I, J) Y (J) = FN (MIN) СУММ = СУММ + Y (Дж) 18 СУММ = СУММ + Г (Дж) * Г (Дж) ICOUNT = ICOUNT + NN ПЕРЕЙДИТЕ К 901 C C СОХРАНЯЙТЕ ОТРАЖЕНИЕ C 19 ДО 20 I = 1, Н 20 P (I, IHI) = PSTAR (I) Y (IHI) = YSTAR 900 СУММ = СУММ + Y (IHI) СУММ = СУММ + Y (IHI) * Y (IHI) 901 JCOUNT = JCOUNT-1 ЕСЛИ (JCOUNT.NE.0) ПЕРЕЙДИТЕ К 1000 C C ПРОВЕРИТЬ, Достигнут ли минимум C ЕСЛИ (ICOUNT.GT.KCOUNT) ПЕРЕЙДИТЕ К 22 JCOUNT = КОНГВЕ CURMIN = (SUMM- (SUM * SUM) / DNN) / DN C C CURMIN ЯВЛЯЕТСЯ РАЗНИЦОМ N + 1 ЗНАЧЕНИЙ FN НА ВЕРСИЯХ C ЕСЛИ (CURMIN.GE.REQMIN) ПЕРЕЙДИТЕ К 1000 C C ФАКТОРНЫЕ ИСПЫТАНИЯ, ЧТОБЫ ПРОВЕРИТЬ, ЯВЛЯЕТСЯ YNEWLO МЕСТНЫМ МИНИМУМОМ C 22 ДО 23 I = 1, N 23 МИН (I) = P (I, IHI) YNEWLO = Y (IHI) ЕСЛИ (ICOUNT.GT.KCOUNT) GOTO 30 DO 24 I = 1, N DEL = ШАГ (I) * 1.D-3 МИН (I) = МИН (I) + DEL Z = FN (МИН) ЕСЛИ (Z.LT.YNEWLO) ПЕРЕЙДИТЕ К 25 MIN (I) = MIN (I) -DEL-DEL Z = FN (МИН) ЕСЛИ (Z.LT.YNEWLO) ПЕРЕЙДИТЕ К 25 24 МИН (I) = МИН (I) + DEL C ДОБАВЛЕННАЯ ЧАСТЬ ДЛЯ ОЦЕНКИ ВАРИАЦИИ ПЕРЕМЕННЫХ 30 ДО 32 I = 1, Н СТАРТ (I) = 0.D0 ДО 32 Дж = 1, НН ЕСЛИ (J.EQ.IHI) ПЕРЕЙТИ 32 Z = DABS (P (I, J) -P (I, IHI)) ЕСЛИ (СТАРТ (I) .LT.Z) СТАРТ (I) = Z 32 ПРОДОЛЖИТЬ C МАКСИМАЛЬНЫЕ РАЗЛИЧИЯ В ПАРАМЕТРАХ ЕСТЬ В НАЧАЛЕ (I) ВОЗВРАЩЕНИЕ C C ПРОЦЕДУРА ПЕРЕЗАГРУЗКИ C 25 ДО 26 I = 1, N 26 ПУСК (I) = МИН (I) DEL = 1.D-3 ICOUNT = ICOUNT + 10000 ПЕРЕЙДИТЕ К 1001 КОНЕЦ

Root APK для Samsung Galaxy Note8 Greatqltechn Sm N9500 ​​[последняя версия]

Samsung Galaxy Note8 Greatqltechn Sm N9500 ​​Root APK

Загрузка и установка корневого APK-файла для Samsung Galaxy Note8 Greatqltechn Sm N9500 ​​с версией Android: Marshmallow / Nougat / Pie / Q / Lollipop / Oreo / KitKat / Lollipop.

Выберите вашу версию android:

Android Gingerbread 2.3.3 — 2.3.7 2010 год

• Имя приложения: AndroidRoot_5.6.APK
• Языки: английский, многоязычный
• Размер файла: 7,49 МБ
• Последнее обновление: 28 июня 2021 г.
• Тип файла: application / vnd.android.package-archive
• Версия: 5.6
• Загрузки: 39159

Загрузить сейчас

Android Ice Cream Sandwich 4.0.3 — 4.0.4 2011 год

• Имя приложения: AndroidRoot_7.3.APK
• Языки: английский, многоязычный
• Размер файла: 3.58 МБ
• Последнее обновление: 28 апреля 2021 г.
• Тип файла: application / vnd.android.package-архив
• Версия: 7.3
• Загрузки: 83732

Загрузить сейчас

Android Jelly Bean 4.1.x 2012 год

• Имя приложения: AndroidRoot_2.4.APK
• Языки: английский, многоязычный
• Размер файла: 8.62 МБ
• Последнее обновление: 09 марта 2021 г.
• Тип файла: application / vnd.android.package-archive
• Версия: 2.4
• Загрузки: 99 996

Загрузить сейчас

Android Jelly Bean 4.2.x 2012 год

• Имя приложения: AndroidRoot_4.3.APK
• Языки: английский, многоязычный
• Размер файла: 2.9 МБ
• Последнее обновление: 18 июля 2021 г.
• Тип файла: application / vnd.android.package-archive
• Версия: 4.3
• Загрузок: 22274

Загрузить сейчас

Android Jelly Bean 4.3 2013 год

• Имя приложения: AndroidRoot_3.7.APK
• Языки: английский, многоязычный
• Размер файла: 2.55 МБ
• Последнее обновление: 13 июня 2021 г.
• Тип файла: application / vnd.android.package-archive
• Версия: 3.7
• Загрузки: 44 156

Загрузить сейчас

Android KitKat 4.4 2013 год

• Имя приложения: AndroidRoot_7.5.APK
• Языки: английский, многоязычный
• Размер файла: 6.06 MB
• Последнее обновление: 09 июн 2021
• Тип файла: application / vnd.android.package-архив
• Версия: 7.5
• Загрузки: 93703

Загрузить сейчас

Android Lollipop 5.0 2014 год

• Имя приложения: AndroidRoot_4.2.APK
• Языки: английский, многоязычный
• Размер файла: 8.76 МБ
• Последнее обновление: 22 апреля 2021 г.
• Тип файла: application / vnd.android.package-archive
• Версия: 4.2
• Загрузки: 59186

Загрузить сейчас

Android Lollipop 5.1 2015 год

• Имя приложения: AndroidRoot_9.8.APK
• Языки: английский, многоязычный
• Размер файла: 2.81 MB
• Последнее обновление: 31 мая 2021 г.
• Тип файла: application / vnd.android.package-archive
• Версия: 9.8
• Загрузки: 59 775

Загрузить сейчас

Android Marshmallow 6.0 2015 год

• Имя приложения: AndroidRoot_4.6.APK
• Языки: английский, многоязычный
• Размер файла: 3.46 МБ
• Последнее обновление: 06 Апрель 2021
• Тип файла: application / vnd.android.package-archive
• Версия: 4.6
• Загрузки: