arxiv: v1 [math.ca] 30 Sep 2014

Transcription

1 SOLVING ABEL INTEGRAL EQUATIONS OF FIRST KIND VIA FRATIONAL ALULUS arxiv: v1 [math.a] 3 Sep 214 S. JAHANSHAHI, E. BABOLIAN, D. F. M. TORRES, AND A. R. VAHIDI Abstract. We give a new method for numerically solving Abel integral equations offirstkind. An estimation forthe errorisobtained. The method isbased on approximations of fractional integrals aputo derivatives. Using trapezoidal rule omputer Algebra System Maple, the exact approximation values of three Abel integral equations are found, illustrating the effectiveness of the proposed approach. 1. Introduction onsider the following generalized Abel integral equation of first kind: (1.1) f(x) = k(x, s)g(s) ds, < α < 1, x b, (x s) α where g is the unknown function to be found, f is a well behaved function, k is the kernel. This equation is one of the most famous equations that frequently appears in many physical engineering problems, like semi-conductors, heat conduction, metallurgy chemical reactions [1, 11]. In experimental physics, Abel s integral equation of first kind (1.1) finds applications in plasma diagnostics, physical electronics, nuclear physics, optics astrophysics [13, 14]. To determine the radial distribution of the radiation intensity of a cylinder discharge in plasma physics, for example, one needs to solve an integral equation (1.1) with α = 1 2. Another example of application appears when one describes velocity laws of stellar, then (1.1) is a fractional integral equation of order 1 α [24]. This problem is a generalization of the tautochrone problem of the calculus of variations, is related with the born of fractional mechanics [28]. The literature on integrals derivatives of fractional order is now vast evolving: see, e.g., [7, 8, 32 34]. The reader interested in the early literature, showing that Abel s integral equations may be solved with fractional calculus, is referred to [9]. For a concise recent discussion on the solutions of Abel s integral equations using fractional calculus see [16]. Many numerical methods for solving (1.1) have been developed over the past few years, such as product integration methods [2,3], collocation methods [5], fractional winds [13]. If k(x, s) = 1 Γ(1 α) 21 Mathematics Subject lassification. 26A33; 45E1. Key words phrases. Abel integral equation; singular integral equations; aputo fractional derivatives; fractional integrals. This is a preprint of a paper whose final definite form is in Journal of King Saud University (Science), ISSN Paper submitted 3/Aug/214; revised 25/Sept/214; accepted for publication 29/Sept/214. orresponding author: Delfim F. M. Torres (delfim@ua.pt), Department of Mathematics, University of Aveiro, Aveiro, Portugal. Tel: , Fax:

2 2 S. JAHANSHAHI, E. BABOLIAN, D. F. M. TORRES, AND A. R. VAHIDI multi step methods [18, 19, 23], backward Euler methods [2], methods based on wavelets [15, 29, 3]. Some semi analytic methods, like the Adomian decomposition method, are also available, which produce a series solution [4]. Unfortunately, 1 the Abel integral equation (1.1) is an ill-posed problem. For k(x, s) = Γ(1 α), Gorenflo [1] presented some numerical methods based on fractional calculus, e.g., using the Grunwald Letnikov difference approximation n ( ) α (1.2) D α f h α ( 1) r f(x rh). r r= If f is sufficiently smooth vanishes at x, then formula (1.2) has accuracy of order O(h 2 ), otherwise it has accuracy of order O(h). On the other h, Lubich [18, 19] introduced a fractional multi-step method for the Abel integral equation of first kind, Plato [23] considered fractional multi-step methods for weakly singular integral equations of first kind with a perturbed right-h side. Liu Tao solved the fractional integral equation, transforming it into an Abel integral equation of second kind [17]. A method based on hebyshev polynomials is given in [1]. Here we propose a method to solve an Abel integral equation of first kind based on a numerical approximation of fractional integrals aputo derivatives of a given function f belonging to n [a,b] (see Theorem 4.2). The structure of the paper is as follows. In Section 2 we recall the necessary definitions of fractional integrals derivatives explain some useful relations between them. Section 3 reviews some numerical approximations for fractional integrals derivatives. The original results are then given in Section 4, where we introduce our method to approximate the solution of the Abel equation at the given nodes we obtain an upper bound for the error. In Section 5 some examples are solved to illustrate the accuracy of the proposed method. 2. Definitions, Relations Properties of Fractional Operators Fractional calculus is a classical area with many good books available. We refer the reader to [2,24]. Definition 2.1. Let α > with n 1 < α n, n N, a < x < b. The left right Riemann Liouville fractional integrals of order α of a given function f are defined by ajxf(x) α = 1 (x t) α 1 f(t)dt Γ(α) xjb α f(x) = 1 (t x) α 1 f(t)dt, Γ(α) x respectively, where Γ is Euler s gamma function, that is, Γ(x) = a b t x 1 e t dt. Definition 2.2. The left right Riemann Liouville fractional derivatives of order α >, n 1 < α n, n N, are defined by d n adx α f(x) = 1 Γ(n α) dx n a (x t) n α 1 f(t)dt

3 SOLVING ABEL INTEGRAL EQUATIONS OF FIRST KIND 3 respectively. xdb α ( 1)n f(x) = Γ(n α) d n dx n b x (t x) n α 1 f(t)dt, Definition 2.3. The left right aputo fractional derivatives of order α >, n 1 < α n, n N, are defined by respectively. ad α xf(x) = 1 Γ(n α) xdbf(x) α = ( 1)n Γ(n α) a b x (x t) n α 1 f (n) (t)dt (t x) n α 1 f (n) (t)dt, Definition 2.4. Let α >. The Grunwald Letnikov fractional derivatives are defined by ( ) α D α f(x) = lim h α ( 1) r f(x rh) h r where D α f(x) = lim h α h [ α r ] r= [ α r r= ] f(x rh), = α(α+1)(α+2) (α+r 1). r! Remark 2.5. The aputo derivatives (Definition 2.3) have some advantages over the Riemann Liouville derivatives (Definition 2.2). The most well known is related with the Laplace transform method for solving fractional differential equations. The Laplace transform of a Riemann Liouville derivative leads to boundary conditions containing the limit values of the Riemann Liouville fractional derivative at the lower terminal x = a. In spite of the fact that such problems can be solved analytically, there is no physical interpretation for such type of boundary conditions. In contrast, the Laplace transform of a aputo derivative imposes boundary conditions involving integer-order derivatives at x = a, which usually are acceptable physical conditions. Another advantage is that the aputo derivative of a constant function is zero, whereas for the Riemann Liouville it is not. For details see [31]. Remark 2.6. The Grunwald Letnikov definition gives a generalization of the ordinary discretization formulas for derivatives with integer order. The series in Definition 2.4 converge absolutely uniformly for each α > for every bounded function f. The discrete approximations derived from the Grunwald Letnikov fractional derivatives, e.g., (1.2), present some limitations. First, they frequently originate unstable numerical methods henceforth usually a shifted Grunwald Letnikov formula is used instead. Another disadvantage is that the order of accuracy of such approximations is never higher than one. For details see [2].

4 4 S. JAHANSHAHI, E. BABOLIAN, D. F. M. TORRES, AND A. R. VAHIDI The following relations between aputo Riemann Liouville fractional derivatives hold [24]: n 1 adxf(x) α = a Dxf(x) α f (k) (a) Γ(k α+1) (x a)k α k= n 1 x Dα b f(x) = xdb α f(x) f (k) (b) Γ(k α+1) (b x)k α. Therefore, if f n [a,b] f (k) (a) =, k =,1,...,n 1, then if f (k) (b) =, k =,1,...,n 1, then k= a Dα x f = ad α x f, xd α bf = x D α bf. Other useful properties of fractional integrals derivatives are: all fractional operators are linear, that is, if L is an arbitrary fractional operator, then L(tf +sg) = tl(f)+sl(g) for all functions f,g n [a,b] or f,g L p (a,b) t,s R; if α,β >, then J α J β = J α+β, D α D β = D α+β ; if f L (a,b) or f n [a,b] α >, then (2.1) ad α x aj α xf = f xd α b xj α b f = f. On the other h, if f n [a,b] α >, then n 1 ajx α a Dα x f(x) = f(x) k= n 1 xjb α x Dα b f(x) = f(x) k= f (k) (a) (x a) k k! f (k) (b) (b x) k. k! From (2.1) it is seen that the aputo fractional derivatives are the inverse operators for the Riemann Liouville fractional integrals. This is the reason why we choose here to use aputo fractional derivatives for solving Abel integral equations. 3. Known Numerical Approximations to Fractional Operators Diethelm [6] (see also [8, pp ]) uses the product trapezoidal rule with respecttotheweightfunction(t k ) α 1 toapproximateriemann Liouvillefractional integrals. More precisely, the approximation tk t (t k u) α 1 f(u)du tk t (t k u) α 1 fk (u)du, where f k is the piecewise linear interpolator for f whose nodes are chosen at t j = jh, j =,1,...,n h = b a n, is considered. Odibat [21,22] uses a modified trapezoidal rule to approximate the fractional integral Jxf(x) α (Theorem 3.1) the aputo fractional derivative adxf(x) α (Theorem 3.2) of order α >.

5 SOLVING ABEL INTEGRAL EQUATIONS OF FIRST KIND 5 Theorem 3.1 (See [21,22]). Let b >, α >, suppose that the interval [,b] is subdivided into k subintervals [x j,x j+1 ], j =,...,k 1, of equal distances h = b k by using the nodes x j = jh, j =,1,...,k. Then the modified trapezoidal rule T(f,h,α) = + h α Γ(α+2) ((k 1)α+1 (k α 1)k α )f()+f(b)) h α Γ(α+2) k 1 ( (k j +1) α+1 2(k j) α+1 +(k j 1) α+1) f(x j ) j=1 is an approximation to the fractional integral J α x f x=b: Furthermore, if f 2 [,b], then J α x f x=b = T(f,h,α) E T (f,h,α). E T (f,h,α) c α f b α h 2 = O(h 2 ), where c α is a constant depending only on α. The following theorem gives an algorithm to approximate the aputo fractional derivative of an arbitrary order α >. Theorem 3.2 (See [21,22]). Let b >, α > with n 1 < α n, n N, suppose that the interval [,b] is subdivided into k subintervals [x j,x j+1 ], j =,...,k 1, of equal distances h = b k by using the nodes x j = jh, j =,1,...,k. Then the modified trapezoidal rule [ h n α ((k (f,h,α) = 1) n α+1 (k n+α 1)k n α) f (n) () Γ(n+2 α) ] k 1 +f (n) ( (b)+ (k j +1) n α+1 2(k j) n α+1 +(k j 1) n α+1) f (n) (x j ) j=1 is an approximation to the aputo fractional derivative Dα xf x=b : Furthermore, if f n+2 [,b], then D α xf x=b = (f,h,α) E (f,h,α). E (f,h,α) c n α f(n+2) b n α h 2 = O(h 2 ), where c n α is a constant depending only on α. In the next section we use Theorem 3.2 to find an approximation solution to a generalized Abel integral equation. The reader interested in other useful approximations for fractional operators is referred to [25 27] references therein. 4. Main Results onsider the following Abel integral equation of first kind: (4.1) f(x) = g(t) (x t) αdt, < α < 1, x b, where f 1 [a,b] is a given function satisfying f() = g is the unknown function we are looking for.

6 6 S. JAHANSHAHI, E. BABOLIAN, D. F. M. TORRES, AND A. R. VAHIDI Theorem 4.1. The solution to problem (4.1) is (4.2) g(x) = sin(απ) π f (t) (x t) 1 αdt. Proof. According to Definition 2.1, we can write (4.1) in the equivalent form Then, using (2.1), it follows that (4.3) Therefore, g(x) = f(x) = Γ(1 α) Jx 1 α g(x). D1 α x 1 x Γ(α)Γ(1 α) f(x) = Γ(1 α)g(x). f (t) sin(απ) x (x t) 1 αdt = π where we have used the identity π = sin(απ)γ(α)γ(1 α). f (t) (x t) 1 αdt, Our next theorem gives an algorithm to approximate the solution (4.2) to problem (4.1). Theorem 4.2. Let < x < b suppose that the interval [,x] is subdivided into k subintervals [t j,t j+1 ], j =,...,k 1, of equal length h = x k by using the nodes t j = jh, j =,...,k. An approximate solution g to the solution g of the Abel integral equation (4.1) is given by [ h α ((k (4.4) g(x) = 1) 1+α (k 1 α)k α) f ()+f (x) Γ(1 α)γ(2+α) ] k 1 ( + (k j +1) 1+α 2(k j) 1+α +(k j 1) 1+α) f (t j ). j=1 Moreover, if f 3 [,x], then g(x) = g(x) 1 Γ(1 α) E(x) with E(x) c α f x α h 2 = O(h 2 ), where c α is a constant depending only on α f = max t [,x] f (t). Proof. We want to approximate the aputo derivative in (4.3), i.e., to approximate 1 g(x) = Γ(1 α) D1 α x f(x) = sin(απ)γ(α) π D1 α x f(x). If the aputo fractional derivative of order 1 α for f, < α < 1, is calculated at collocation nodes t j, j =,...,k, then the result is a direct consequence of Theorem Illustrative Examples We exemplify the approximation method given by Theorem 4.2 with three Abel integral equations whose exact solutions are found from Theorem 4.1. The computations were done with the omputer Algebra System Maple. The complete code is provided in Appendix A.

7 SOLVING ABEL INTEGRAL EQUATIONS OF FIRST KIND 7 Example 5.1. onsider the Abel integral equation (5.1) e x 1 = g(t) (x t) 1/2dt. The exact solution to (5.1) is given by Theorem 4.1: (5.2) g(x) = ex π erf( x), where erf(x) is the error function, that is, erf(x) = 2 e t2 dt. π In Table 1 we present the approximate values for g(x i ), x i =.1,.2,.3, obtained from Theorem 4.2 with k = 1, 1, 1. Table 1. Approximated values g(x i ) (4.4) to the solution (5.2) of (5.1) obtained from Theorem 4.2 with k = 1, 1, 1. x i k = 1 k = 1 k = 1 Exact solution (5.2) Error k= Example 5.2. onsider the following Abel integral equation of first kind: (5.3) x = g(t) (x t) 4/5dt. The exact solution (4.2) to (5.3) takes the form (5.4) g(x) = 5 sin ( ) π 5 x 4/5. 4 π We can see the numerical approximations of g(x i ), x i =.4,.5,.6, obtained from Theorem 4.2 with k = 1, 1 in Table 2. Table 2. Approximated values g(x i ) (4.4) to the solution (5.4) of (5.3) obtained from Theorem 4.2 with k = 1, 1. x i k = 1 k = 1 Exact solution (5.4) Error k= Example 5.3. onsider now the Abel integral equation (5.5) x 7/6 = Theorem 4.2 asserts that (5.6) g(x) = π g(t) (x t) 1/3dt. t 1/6 (x t) 2/3dt.

8 8 S. JAHANSHAHI, E. BABOLIAN, D. F. M. TORRES, AND A. R. VAHIDI The numerical approximations of g(x i ), x i =.6,.7,.8, obtained from Theorem 4.2 with k = 1,1,1, are given in Table 3. Table 3. Approximated values g(x i ) (4.4) to the solution (5.6) of (5.5) obtained from Theorem 4.2 with k = 1,1,1. x i k = 1 k = 1 k = 1 Exact solution (5.6) Error k= Appendix A. Maple ode for Examples of Section 5 We provide here all the definitions computations done in Maple for the problems considered in Section 5. The definitions follow closely the notations introduced along the paper, should be clear even for readers not familiar with the omputer Algebra System Maple. > # Solution given by Theorem 4.1 > g := (f, alpha, x) -> sin(alpha*pi)*(int((diff(f(t), t))/(x-t)^(1-alpha), t =..x))/pi: > # Approximation given by Theorem 4.2 > gtilde := (f, alpha, h, k, x) -> h^alpha *(((k-1)^(1+alpha)-(k-1-alpha)*k^alpha) *(D(f))()+(D(f))(x)+sum(((k-j+1)^(1+alpha)-2*(k-j)^(1+alpha) +(k-j-1)^(1+alpha))*(d(f))(j*h), j = 1.. k-1)) /(GAMMA(1-alpha)*GAMMA(2+alpha)): > # Example 5.1 > f1 := x -> exp(x)-1: > g(f1, 1/2, x); erf( x)e x π > ExactValues1 := evalf([g(f1,1/2,.1), g(f1,1/2,.2), g(f1,1/2,.3)]); [ , , ] > ApproximateValues1 := k -> evalf([gtilde(f1, 1/2,.1/k, k,.1), gtilde(f1, 1/2,.2/k, k,.2), gtilde(f1, 1/2,.3/k, k,.3)]): > ApproximateValues1(1); > ApproximateValues1(1); > ApproximateValues1(1); > # Example 5.2 > f2 := x -> x: > g(f2, 4/5, x); [ , , ] [ , , ] [ , , ] 5 sin ( 1 π) x 4/5 5 4 π

9 SOLVING ABEL INTEGRAL EQUATIONS OF FIRST KIND 9 > ExactValues2 := evalf([g(f2,4/5,.4), g(f2, 4/5,.5), g(f2, 4/5,.6)]); [ , , ] > ApproximateValues2 := k -> evalf([gtilde(f2, 4/5,.4/k, k,.4), gtilde(f2, 4/5,.5/k, k,.5), gtilde(f2, 4/5,.6/k, k,.6)]): > ApproximateValues2(1); > ApproximateValues2(1); > # Example 5.3 [ , , ] [ , , ] > f3 := x -> x^(7/6): > g(f3, 1/3, x); x 7 t 1/ (x t) 2/3dt 2 π > ExactValues3 := evalf([g(f3,1/3,.6), g(f3,1/3,.7), g(f3,1/3,.8)]); [ , , ] > ApproximateValues3 := k -> evalf([gtilde(f3, 1/3,.6/k, k,.6), gtilde(f3, 1/3,.7/k, k,.7), gtilde(f3, 1/3,.8/k, k,.8)]): > ApproximateValues3(1); > ApproximateValues3(1); > ApproximateValues3(1); > ApproximateValues3(1); [ , , ] [ , , ] [ , , ] [ , , ] Acknowledgements This work is part of first author s PhD project. Partially supported by Azad University, Iran; IDMA FT, Portugal, within project PEst-OE/MAT/UI416/ 214. Jahanshahi was supported by a scholarship from the Ministry of Science, Research Technology of the Islamic Republic of Iran, to visit the University of Aveiro. The hospitality the excellent working conditions at the University of Aveiro are here gratefully acknowledged. The authors would like also to thank three anonymous referees for valuable comments. References [1] Z. Avazzadeh, B. Shafiee G. B. Loghmani, Fractional calculus of solving Abel s integral equations using hebyshev polynomials, Appl. Math. Sci. (Ruse) 5 (211), no , [2]. T. H. Baker, The numerical treatment of integral equations, larendon Press, Oxford, [3] P. Baratella A. P. Orsi, A new approach to the numerical solution of weakly singular Volterra integral equations, J. omput. Appl. Math. 163 (24), no. 2, [4] L. Bougoffa, A. Mennouni R.. Rach, Solving auchy integral equations of the first kind by the Adomian decomposition method, Appl. Math. omput. 219 (213), no. 9,

10 1 S. JAHANSHAHI, E. BABOLIAN, D. F. M. TORRES, AND A. R. VAHIDI [5] H. Brunner, ollocation methods for Volterra integral related functional differential equations, ambridge Monographs on Applied omputational Mathematics, 15, ambridge Univ. Press, ambridge, 24. [6] K. Diethelm, An algorithm for the numerical solution of differential equations of fractional order, Electron. Trans. Numer. Anal. 5 (1997), Mar., 1 6. [7] K. Diethelm, N. J. Ford, A. D. Freed Yu. Luchko, Algorithms for the fractional calculus: a selection of numerical methods, omput. Methods Appl. Mech. Engrg. 194 (25), no. 6-8, [8] K. Diethelm A. D. Freed, The FracPEE subroutine for the numerical solution of differential equations of fractional order. In: Proc. of Forschung und Wwissenschaftliches Rechnen: Beiträge Zum Heinz-Billing-Preis 1998, 22, pp [9] I. M. Gel f G. E. Shilov, Generalized functions. Vol. I, Translated by Eugene Saletan, Academic Press, New York, [1] R. Gorenflo, Abel integral equations with special emphasis on applications, Lectures in Mathematical Sciences, Vol. 13, The University of Tokyo, Graduate School of Mathematical Sciences, Tokyo, [11] R. Gorenflo S. Vessella, Abel integral equations, Lecture Notes in Mathematics, 1461, Springer, Berlin, [12] D. Kincaid W. heney, Numerical analysis, Brooks/ole, Pacific Grove, A, [13] O. Knill, R. Dgani M. Vogel, A new approach to Abel s integral operator its application to stellar winds, Astronom. Astrophys. 274 (1993), no. 3, [14] E. L. Kosarev, Applications of integral equations of the first kind in experiment physics, omput. Phys. ommun. 2 (198), [15] Ü. Lepik, Solving fractional integral equations by the Haar wavelet method, Appl. Math. omput. 214 (29), no. 2, [16] M. Li W. Zhao, Solving Abel s type integral equation with Mikusinski s operator of fractional order, Adv. Math. Phys. 213 (213), Art. ID 86984, 4 pp. [17] Y. Liu L. Tao, Mechanical quadrature methods their extrapolation for solving first kind Abel integral equations, J. omput. Appl. Math. 21 (27), no. 1, [18] h. Lubich, Fractional linear multistep methods for Abel-Volterra integral equations of the second kind, Math. omp. 45 (1985), no. 172, [19] h. Lubich, Discretized fractional calculus, SIAM J. Math. Anal. 17 (1986), no. 3, [2] A. B. Malinowska D. F. M. Torres, Introduction to the fractional calculus of variations, Imp. oll. Press, London, 212. [21] Z. Odibat, Approximations of fractional integrals aputo fractional derivatives, Appl. Math. omput. 178 (26), no. 2, [22] Z. M. Odibat, omputational algorithms for computing the fractional derivatives of functions, Math. omput. Simulation 79 (29), no. 7, [23] R. Plato, Fractional multistep methods for weakly singular Volterra integral equations of the first kind with perturbed data, Numer. Funct. Anal. Optim. 26 (25), no. 2, [24] I. Podlubny, Fractional differential equations, Mathematics in Science Engineering, 198, Academic Press, San Diego, A, [25] S. Pooseh, R. Almeida D. F. M. Torres, Approximation of fractional integrals by means of derivatives, omput. Math. Appl. 64 (212), no. 1, arxiv: [26] S. Pooseh, R. Almeida D. F. M. Torres, Numerical approximations of fractional derivatives with applications, Asian J. ontrol 15 (213), no. 3, arxiv: [27] S. Pooseh, R. Almeida D. F. M. Torres, Fractional order optimal control problems with free terminal time, J. Ind. Manag. Optim. 1 (214), no. 2, arxiv: [28] F. Riewe, Mechanics with fractional derivatives, Phys. Rev. E (3) 55 (1997), no. 3, part B, [29] H. Saeedi, M. Mohseni Moghadam, N. Mollahasani, G. N. huev, A AS wavelet method for solving nonlinear Fredholm integro-differential equations of fractional order, ommun. Nonlinear Sci. Numer. Simul. 16 (211), no. 3, [3] H. Saeedi, N. Mollahasani, M. Mohseni Moghadam G. N. huev, An operational Haar wavelet method for solving fractional Volterra integral equations, Int. J. Appl. Math. omput. Sci. 21 (211), no. 3, [31] E. Sousa, How to approximate the fractional derivative of order 1 < α 2, Internat. J. Bifur. haos Appl. Sci. Engrg. 22 (212), no. 4, 12575, 13 pp.

11 SOLVING ABEL INTEGRAL EQUATIONS OF FIRST KIND 11 [32] V. E. Tarasov, Review of some promising fractional physical models, Internat. J. Modern Phys. B 27 (213), no. 9, 1335, 32 pp. [33] J. A. Tenreiro Machado, V. Kiryakova F. Mainardi, Recent history of fractional calculus, ommun. Nonlinear Sci. Numer. Simul. 16 (211), no. 3, [34] J. Wang, Y. Zhou, W. Wei H. Xu, Nonlocal problems for fractional integrodifferential equations via fractional operators optimal controls, omput. Math. Appl. 62 (211), no. 3, Department of Mathematics, Science Research Branch, Islamic Azad University, Tehran, Iran address: s.jahanshahi@iausr.ac.ir Department of Mathematics, Science Research Branch, Islamic Azad University, Tehran, Iran address: Babolian@khu.ac.ir enter for Research Development in Mathematics Applications (IDMA), Department of Mathematics, University of Aveiro, Aveiro, Portugal address, orresponding author: delfim@ua.pt Department of Mathematics, Shahre Rey Branch, Islamic Azad University, Tehran, Iran address: alrevahidi@iausr.ac.ir