#!/usr/bin/env python3
# -*- coding: utf-8 -*-
"""
Created on Fri Oct 28 09:15:52 2016
@author: andreaskuk
"""
import numpy as np
import matplotlib.pyplot as plt
"""
function: Y=Y_0(cos(\theta)^(-f)*exp(b(1-1/cos(\theta))))
mit Y_0= 0.11, f=1.9, b=0.19 von 0°-90°
"""
y0, f ,b, = 0.11, 1.9,0.19
x= np.linspace(0,90,10000)
#y = y0*np.power(np.cos(x),-f)*np.exp(b*(1-1/np.cos(x)))
y = y0*np.cos(np.radians(x))**(-f)*np.exp(b*(1-1/np.cos(np.radians(x))))
plt.title('Plot Sputter Yield')
plt.xlabel('0°-90°')
plt.ylabel('y')
plt.grid(True)
plt.plot(x,y, label='$Y=$Y_0*(cos(\Theta))^(-f)*exp(b(1-(1/ cos(\Theta))))')
plt.legend(loc='upper left')
plt.show()
class Kukplot:
def __init__(self,y0,f,b):
#self.x=np.linespace(0,90,101)
self.y0=y0
self.f=f
self.b=b
# self.y=y0*np.cos(np.radians(x))**(-f)*np.exp(b*(1-1/np.cos(np.radians(x))))
def __call__(self,x):
self.x=x
y = self.y0*np.cos(np.radians(self.x))**(-self.f)*np.exp(self.b*(1-1/np.cos(np.radians(self.x))))
return y
x=np.linspace(0,90,101)
sjsj = Kukplot(0.11,1.9,0.19)
y= sjsj(np.linspace(0,90,101))
plt.title('Plot Sputter Yield')
plt.xlabel('0°-90°')
plt.ylabel('y')
plt.grid(True)
plt.plot(x,y, label='$Y=$Y_0*(cos(\Theta))^(-f)*exp(b(1-(1/ cos(\Theta))))')
plt.legend(loc='upper left')
plt.show()
