Архив: Перевести короткий код с phyton на c/c++
Бюджет
8$
/ 500
руб
в месяц
Создан: 7 лет назад
Закрыт
- Описание
- Нужно сегодня
Вот код
%matplotlib inline
import numpy as np
import math
import matplotlib.pyplot as plt
import pylab
from sympy import *
def f(x):
if v==1:
return math.sin(abs(x))
elif v==2:
return math.sin(x*math.pi)
elif v==3:
if (x>0):
return x*(-math.log(x))**3
else:
return 0
else:
if (x>0):
return x*math.log(x)
else:
return 0
def Rect(n):
x=np.linspace(a,b,n+1)
h=(b-a)/n
S=0
for i in range(n):
S+=f(x[i]+h/2)
S*=h
return S
def Trap(n):
x=np.linspace(a,b,n+1)
h=(b-a)/n
S=1/2*f(x[0])
for i in range(1,n,1):
S+=f(x[i])
S+=1/2*f(x[n])
S*=h
return S
def Simpson(n):
x=np.linspace(a,b,n+1)
h=(b-a)/n
S=0
for i in range(n):
S+=f(x[i])+4*f(x[i]+h/2)+f(x[i+1])
S=S*h/6
return S
def Gauss2(n):
S=0
x=np.linspace(a,b,n+1)
h=(b-a)/n
for i in range(n):
a1=x[i]
b1=x[i+1];
x1=(a1+b1)/2-(b1-a1)*math.sqrt(3)/6
x2=(a1+b1)/2+(b1-a1)*math.sqrt(3)/6
S+=f(x1)+f(x2)
S*=h/2
return S
def Gauss3(n):
S=0
x=np.linspace(a,b,n+1)
h=(b-a)/n;
for i in range(n):
a1=x[i];
b1=x[i+1];
x1=(a1+b1)/2-(b1-a1)/2*math.sqrt(3/5)
x2=(a1+b1)/2
x3=(a1+b1)/2+(b1-a1)/2*math.sqrt(3/5)
S+=5*f(x1)+8*f(x2)+5*f(x3)
S*=h/(2*9)
return S
v=1
n=100
if v==1:
a=-1; b=2;
I=2-math.cos(-1)-math.cos(2)
elif v==2:
a=-1; b=1;
I=0
elif v==3:
a=0; b=1;
I=3/8
else:
a=0; b=1;
I=-1/4
print(I)
print(Rect(n),"Rectangles")
print(Trap(n),"Trapeze")
print(Simpson(n),"Simpson")
print(Gauss2(n),"Gauss-2")
print(Gauss3(n),"Gauss-3")
# A=[2**i for i in range(1,30,1)]
A = np.arange(1, 100)
xR=[i for i in A]
xT=[i+1 for i in A]
xS=[2*i+1 for i in A]
xG2=[2*i for i in A]
xG3=[3*i for i in A]
yR=[Rect(i) for i in A]
yT=[Trap(i) for i in A]
yS=[Simpson(i) for i in A]
yG2=[Gauss2(i) for i in A]
yG3=[Gauss3(i) for i# A=[2**i for i in range(1,30,1)]
A = np.arange(1, 100)
xR=[i for i in A]
xT=[i+1 for i in A]
xS=[2*i+1 for i in A]
xG2=[2*i for i in A]
xG3=[3*i for i in A]
yR=[Rect(i) for i in A]
yT=[Trap(i) for i in A]
yS=[Simpson(i) for i in A]
yG2=[Gauss2(i) for i in A]
yG3=[Gauss3(i) for i in A]
plt.axhline(y=I,color='m',label="In tegral")
plt.plot(xR,yR,label="Rectangles")
plt.plot(xT,yT,label="Trapeze")
plt.plot(xS,yS,label="Simpson")
plt.plot(xG2,yG2,label="Gauss-2")
plt.plot(xG3,yG3,label="Gauss-3")
plt.grid()
plt.legend()
plt.xscale("log")
#print(yG2) in A]
plt.axhline(y=I,color='m',label="In tegral")
plt.plot(xR,yR,label="Rectangles")
plt.plot(xT,yT,label="Trapeze")
plt.plot(xS,yS,label="Simpson")
plt.plot(xG2,yG2,label="Gauss-2")
plt.plot(xG3,yG3,label="Gauss-3")
plt.grid()
plt.legend()
plt.xscale("log")
#print(yG2)
eR=[abs(yR[i]-I) for i in range(len(A))]
eT=[abs(yT[i]-I) for i in range(len(A))]
eS=[abs(yS[i]-I) for i in range(len(A))]
eG2=[abs(yG2[i]-I) for i in range(len(A))]
eG3=[abs(yG3[i]-I) for i in range(len(A))]
plt.plot(xR,eR,label="Error Rectangles")
plt.plot(xT,eT,label="Error Trapeze")
plt.plot(xS,eS,label="Error Simpson")
plt.plot(xG2,eG2,label="Error Gauss-2")
plt.plot(xG3,eG3,label="Error Gauss-3")
plt.yscale("log")
plt.xscale("log")
plt.grid()
plt.legend()
ER=[1/3*abs(Rect(i)-Rect(2*i)) for i in A]
ET=[1/3*abs(Trap(i)-Trap(2*i)) for i in A]
ES=[1/15*abs(Simpson(i)-Simpson(2*i)) for i in A]
EG2=[1/15*abs(Gauss2(i)-Gauss2(2*i)) for i in A]
EG3=[1/63*abs(Gauss3(i)-Gauss3(2*i)) for i in A]
plt.plot(xR,ER,label="Error Runge Rectangles")
plt.plot(xT,ET,label="Error Runge Trapeze")
plt.plot(xS,ES,label="Error Runge Simpson")
plt.plot(xG2,EG2,label="Error Runge Gauss-2")
plt.plot(xG3,EG3,label="Error Runge Gauss-3")
plt.yscale("log")
plt.xscale("log")
plt.grid()
plt.legend()
- Категория
Похожие проекты
0
3
Необходимо разработать систему, которая: Подключается к API Diagnocat и iDent (через токены доступа). Получает данные о пациенте, результатах диагностики, снимках и назначениях. Автоматически формирует структурированный отчет по шаблону (включая врачей, этапы лечения, стоимость и описание процедур).......
2 месяца назад
0
0
Необходимо связать БД 1C:ERP с БД Outlook, MS Access, Битрикс24. Например при вводе нового сотрудника в 1С:ERP он должен автоматически появляться в вышеперечисленных БД....
3 месяца назад
0
1
Требуется реализация отчёта на кастомизированной 1С:УПП 1. 2, ТЗ и контрольный пример во вложении, 80% кода должно быть покрыто модульными тестами....
3 месяца назад
$7
0
2
Мне нужно копировать мой текст и вставлять в одну группу ВК с открытой стеной для публикации рекламы. Задача скрипта/сервиса/программы - примерно каждые 10 секунд брать мой текст и вставлять в эту группу. Я пробовал через одну программу, которая записывает действия мышки и потом повторяет, но она как-то......
5 месяцев назад
0
0
Написать прошивку для A94B114, используя 2 подобные рабочие прошивки. Все файлы есть....
6 месяцев назад
$21
0
2
Нужно разархивировать файл архив. ujems, есть распаковщик, но выдает ошибку. Прислать пруфы(скрины)....
7 месяцев назад
0
1
Входные: Имеется скрипт по обработке email писем и добавление их в систему Redmine Проблема: Система обработки этих писем написана некорректная и негибкая Задача: Переписать логику обработки входящих писем, чтобы исключить их потерю при незначительных отклонениях от ожидаемой структуры....
7 месяцев назад
$37
0
1
- В обязанности сотрудника входит поддержание мобильного приложения системы автоматизации бизнеса. - Работа сдельная: задачи выдаются на оценку, далее после согласования передаются на реализацию. - Заключаем договор и NDA (соглашение о неразглашении)....
7 месяцев назад
$71
0
0
Информация о проекте: Проект включает в себя React Native приложение с Expo внутри которого находится webview с сайтом На устройствах с ОС IOS проблема с воспроизведением видео....
8 месяцев назад
0
0
Разрабатываю десктопные программы на языке Java, а так же Java + SQL....
8 месяцев назад
0
0
0
6
1. Чтобы можно было писать описание события (без заголовка) - превью в списке берется N количества букв с начала. (в списке писать дату события и сколько осталось до него) . Выставление даты и времени напоминания и ИНТЕРВАЛА ПОВТОРА ОТ 1 МИНУТЫ! Выставлять количество повторов, далее напоминания прекращаются.......
9 месяцев назад
0
4
Сайт arg2028t. beget. tech/. Нужно сделать вот такой калькулятор. airprint. by/outdoor/wide-format. и еще некоторые доработки....
9 месяцев назад
$214
0
1
0
7
Краткое описание: Мы ищем исполнителя для тестирования функционала мобильного приложения на платформе Android. Задача включает проверку основных функций приложения, таких как работа комнаты, передача видео, звука и чата, а также создание подробной отчетности о найденных ошибках. Требования к исполнителю:......
10 месяцев назад
0
12
$143
0
0
Oпpocы (1008pyб/чaс) Заходитe на сaйт: gonsù....
11 месяцев назад
0
4
Требуется создать скрипт для premiere pro. Данная суть скрипта: в выделеном бине нужно что бы автоматически проставлялся in out ровно по серидине файла так что бы выделенная область была длиной 3 секунды. И требуется кнопка что бы работало через окно expressions. Проблема текущего в том что он не работает......
11 месяцев назад
0
2
Нужно создать приложения под гемблинг/беттинг, на котлин. Приложение должно работать Firebase. Детали вышлю при обсуждение....
11 месяцев назад
$21
0
4
Нужно разархивировать файл archive. ujems, есть распаковщик, но выдает ошибку. Прислать пруфы ....
1 год назад
0
3
Разработать впн сервис ( Приложения на ios и android + сайт для оплаты подписки). Настройка Серверов. Оплата по договоренности. Работаем черед безопасную сделку....
1 год назад
Похожая удалённая работа
- Удаленная работа для веб-программиста
- Удаленная работа для разработчика баз данных
- Удаленная работа для прикладного программиста
- Удаленная работа для системного программиста
- Удаленная работа для разработчика игр
- Удаленная работа для разработчика мобильных приложений
- Удаленная работа для тестировщика ПО
- Удаленная работа для 1С программиста
- Удаленная работа для разработчика встраиваемых систем
- Удаленная работа для разработчика CRM и ERP
Перетащите