showing error like limit of integration must be scalars..please help

clc clear all

Pc=1.6*10^9 Pa=300*10^9 Ps=1*10^9 r=5 h=1 d0=0.4 d=1:10:1000 j= 2:2:3000 j1=2:2:3000 s=(d./r) w1=4.3 E0=8.85*10^-12 Eb=6.830890.*(w1.^(3/2)) A1=d.^2 B1=(d./d0).^(1/3) l=(0.5.*E0.*Eb.^2./B1.^2).*10^18 m=Pc n=(4.*Pa./(1+s).^4) x2=min(m,n) Es=min(x2,l)

Ea= (Es.*(d+r).^2)./(2.*r.^2)

u1= quad('0.5.*E0.*Ea.^2.*(d).^2.*y.^2./(d).^3.*(d+r).^2',r,r+d)