import numpy as np
from matplotlib.pylab import *
from math import *

n=100
p=0.5
# s(x)
x=np.linspace(-2,2,n)
y=np.zeros(n)

for i in range(n):
y[i]=x[i]**2

# (X(x),Y(x))
yy=np.zeros(n)
xx=np.zeros(n)

for i in range(n):
xx[i]=x[i]-p*2*x[i]/sqrt(1+4*x[i]**2)
yy[i]=x[i]**2+p/sqrt(1+4*x[i]**2)


#r(x)

def vinkel(h): #buelengden fra oppgave c
return 1.0/4.0*(2*h*sqrt(1+4*h**2)+asinh(2*h))/p

ry=np.zeros((n),float)
rx=np.zeros((n), float)
for i in range(n):
rx[i]=p*((2*x[i]/sqrt(1+4*x[i]**2))*cos(vinkel(x[i]))+(-1.0)/sqrt(1+4*x[i]**2)*sin(vinkel(x[i])))+xx(x[i])
ry[i]=p*((-1)*(2*x[i]/sqrt(1+4*x[i]**2))*sin(vinkel(x[i]))+(-1.0)/sqrt(1+4*x[i]**2)*sin(vinkel(x[i])))+yy(x[i])

plot(rx,ry)
plot(xx,yy)
plot(x,y)
show()

Модератор: Тема перенесена из форума "Программирование".