program lefficacie
    implicit none
    real(kind=8):: N1,N2,N3,N4,N5,D1,D2,Leff,dN,dD,dLeff
    real(kind=8), dimension(2) :: mcil,ma,Lasta,dA,h,R
    real(kind=8), dimension(2):: dM

    mcil=(/0.0755, 0.0001/)
    ma = (/0.0256, 0.0001/)
    Lasta = (/0.381, 0.002/)
    dA = (/0.166, 0.003/)
    h = (/0.01825, 0.00001/)
    R = (/0.01160, 0.00001/)
    dM(2) = 0.002

    print *, "Inserisci il valore di dM: "
    read *, dM(1)

    N1 = mcil(1)*dM(1)**2
    N2 = mcil(1)/12*h(1)**2
    N3 = mcil(1)/12*3*R(1)**2
    N4 = ma(1)*dA(1)**2
    N5 = ma(1)/12*Lasta(1)**2
    D1 = mcil(1)*dM(1)
    D2 = ma(1)*dA(1)
    Leff = (N1+N2+N3+N4+N5)/(D1+D2)

    dN = ((mcil(2)/mcil(1))+(2*dM(2)/dM(1)))*N1 + (mcil(2)/mcil(1)+(2*h(2)/h(1)))*N2 &
    +(mcil(2)/mcil(1)+2*R(2)/R(1))*N3 + (ma(2)/ma(1)+2*dA(2)/dA(1))*N4 + (ma(2)/ma(1)+2*Lasta(2)/Lasta(1))*N5

    dD = (mcil(2)/mcil(1)+dM(2)/dM(1))*D1 + (ma(2)/ma(1)+dA(2)/dA(1))*D2

    dLeff = (dN/(N1+N2+N3+N4+N5)+dD/(D1+D2))*Leff

    print *, Leff,"+- ",dLeff

end program