clear
clc
for i=0:8
    R(i+1)=0.06+i*0.005
end
for j=0:6
    C(j+1)=1200-100*j
end
q=0.32
for i=0:8
    for j=0:6
        Y(i+1, j+1) = (((9/16)*(q^2/(C(j+1)^2*R(i+1)^2)))^0.33)*1000
    end
end
for i=0:8
    for j=0:6
        B(i+1, j+1) = ((6*q*(R(i+1)/C(j+1)))^0.33)*1000
    end
end
for i=0:8
    for j=0:6
        σ(i+1, j+1) = C(j+1)*Y(i+1,j+1)*0.001
    end
end
for i=0:8
    for j=0:6
        τ(i+1,j+1) = 0.19*q*R(i+1)*10^4
    end
end
for i=1:5
    Nu(i) = i
end
v = 1.32
for i=1:5
    for j=0:6
        Mσ(j+1,i+1)=2*v*q*10^4*(Nu(i)/C(j+1))
    end
end
printf(" Исследование влияния параметров\n Варьируем радиус валов и жестковть покрытия\n")
printf("i:=0:8   ")
printf("R(i+1) = 0.06 + i*0.005")
printf("j:=0:8   ")
printf("C(j) = 1200 - 100*j\n\n")
printf("Влияния радиуса валов и жесткости покрытия на осадку покрытия, мм и ширину полости контакта\n")
printf("Y(i+1, j+1) = (((9/16)*(q^2/(C(j+1)^2*R(i+1)^2)))^0.33)*1000\n\n")
printf("  R(i)   C(j)\n")
for i=0..6
    printf(" %.3f   %.0f\n", R(i+1), C(i+1))
end
for i=7..8
    printf(" %.3f", R(i+1))
end
printf("\n")
printf("Y=\n")
for i=1..7
    printf("  %.3f   %.3f   %.3f   %.3f   %.3f   %.3f   %.3f\n", Y(i,1),Y(i,2),Y(i,3),Y(i,4),Y(i,5),Y(i,6),Y(i,7))
end
printf("B(i, j) = ((6*q*(R(i+1)/C(j+1)))^0.33)*1000\n\n")
printf("B=\n")
for i=1..7
    printf("  %.3f   %.3f   %.3f   %.3f   %.3f   %.3f   %.3f\n", B(i,1),B(i,2),B(i,3),B(i,4),B(i,5),B(i,6),B(i,7))
end
printf("\n")
printf("\n")
printf("  Влияние радиуса валов и жесткости резинового покрытия на нормальные и касательные напряжения в контакте\n")
printf("σ(i, j) = C(j+1)*Y(i+1,j+1)*0.001\n\n")
printf("σ=\n")
for i=1..7
    printf("  %.3f   %.3f   %.3f   %.3f   %.3f   %.3f   %.3f\n", σ(i,1),σ(i,2),σ(i,3),σ(i,4),σ(i,5),σ(i,6),σ(i,7))
end
printf("\n")
printf("\n")
printf("τ(i,j) = 0.19*q*R(i+1)*10^4\n")
printf("τ=\n")
for i=1..7
    printf("  %.3f   %.3f   %.3f   %.3f   %.3f   %.3f   %.3f\n", τ(i,1),τ(i,2),τ(i,3),τ(i,4),τ(i,5),τ(i,6),τ(i,7))
end
printf("\n")
printf("\n")
printf("  Влияние жесткости и вязкости покрытия на момент качения нормальных сил\n")
printf("i:=1:5   ")
printf("Nu(i)=i   ")
printf("Mσ(j,i-1)=2*v*q*10^4*(Nu(i)/C(j))   ")
printf("v:=1.32\n")
printf("\n")
printf("\n")
for i=1..7
    printf("  %.3f   %.3f   %.3f   %.3f   %.3f\n", Mσ(i,1),Mσ(i,2),Mσ(i,3),Mσ(i,4),Mσ(i,5))
end
scf(0)
plot(R,Y(:,1),'r',R,Y(:,3),'b',R,Y(:,5),"g",R,Y(:,7),'--')
xtitle('','R','Y')
legend('Y(:,0)','Y(:,2)','Y(:,4)','Y(:,6)')
scf(1)
plot(R,B(:,1),'r',R,B(:,3),'b',R,B(:,5),'g',R,B(:,7),'--')
xtitle('','R','B')
legend('B(:,0)','B(:,2)','B(:,4)','B(:,6)')
scf(2)
plot(R,σ(:,1),'r',R,σ(:,3),'b',R,σ(:,5),'g',R,σ(:,7),'--')
xtitle('','R','σ')
legend('σ(:,0)','σ(:,2)','σ(:,4)','σ(:,6)')
scf(3)
plot(R,τ(:,1),'r',R,τ(:,3),'b',R,τ(:,5),'g',R,τ(:,7),'--')
xtitle('','R','τ')
legend('τ(:,0)','τ(:,2)','τ(:,4)','τ(:,6)')
scf(4)
plot(R,Mσ(:,1),'r',R,Mσ(:,3),'b',R,Mσ(:,5),'g',R,Mσ(:,7),'--')
xtitle('','R','Mσ')
legend('Mσ(:,0)','Mσ(:,2)','Mσ(:,4)','Mσ(:,6)')