
    h                     P    S SK r S SKrS SKJr  S SKJr  S rS r	S rS r
S rS rg)	    N)p_rootsc           
          X!-
  U-  n[         R                  " [        US-   5       Vs/ s H	  oQXE-  -   PM     sn5      nUS-  U " US   5      U " US   5      -   [        SU " USS 5      -  5      -   -  $ s  snf N      r   nparrayrangesumfabNdxixjs          /home/user/lab7/myquad.pyqtrapnr      w    	QwB	51:.:aRT6:.	/Ba42a5!BrF)#C!BqH+$6677 /   A4c                 .   X!-
  U-  n[         R                  " [        US-   5       Vs/ s H	  oQXE-  -   PM     sn5      nUS-  U " US   5      U " US   5      -   [        SU " USSS2   5      -  5      -   [        SU " USSS2   5      -  5      -   -  $ s  snf )Nr      r   r      r   r	   )r   r   r   r   r   r   xis          r   qsimpnr      s    	QwB	51:.:aRT6:.	/BqD1RU8AbfI%c!Ab2ajM/&::3q2a1f:;OOPP /s   Bc           
          X!-
  U-  n[         R                  " [        US-   5       Vs/ s H	  oQXE-  -   PM     sn5      nUS-  U " US   5      U " US   5      -   [        SU " USS 5      -  5      -   -  $ s  snf r   r	   r   s          r   r   r      r   r   c           	         SnX!-
  nUS-  U " U5      U " U5      -   -  nSn[        U5       H  n[        R                  " [        SU-  5       V	s/ s H  oUS-  -   XY-  -   PM     sn	5      n
US-  US-  [        U " U
5      5      -  -   n[	        X5        US:  a  [        Xv-
  5      U:  a  Us  $ US-  nUnM     [	        S5        g s  sn	f )Nd   r   r   r   g      ?tol unreached)r   r
   r   r   printabs)r   r   r   tolmaxphI0Ipr   r   s              r   qtrapzr+      s    	DsAsAaD1IBA4[XX%1+6+Q1uQSy+67
Q$1S2Z
a
Q33qt9S=H	3  
/ 7s   C
c                     SnSn[        U5       H;  n[        XUSU-  5      n[        Xg5        US:  a  [        Xu-
  5      U:  a  Us  $ UnM=     [        S5        g )Nr!   r   r   r   r"   )r   r   r#   r$   )r   r   r   r%   r&   r(   r*   r)   s           r   qtrapz0r-   0   s]    	DB4[
Qq!t
a
Q33qt9S=H  
/    c           
          SnSn[        SU5       H]  n[        U5      u  pxX!-
  S-  [        X" X!-
  S-  U-  X-   S-  -   5      -  5      -  n	[        Xi5        [	        XY-
  5      U:  a  U	s  $ U	nM_     [        S5        g )Nr!   r   r   r   r"   )r   r   r   r#   r$   )
r   r   r   r%   maxnr(   nxwr)   s
             r   quadzr4   @   s    	DB1T]AJ3'#a13'!)QS!G+,,-
-a
rt9S=H  
/r.   )mathnumpyr
   matplotlib.pyplotpyplotpltscipy.specialr   r   r   r+   r-   r4    r.   r   <module>r<      s0       !8Q8"
 r.   