Контрольная работа: Решение задачи с помощью программ Mathcad и Matlab Моделирование движения

q=10^(-9);

qp=10^(-6);

r0=1;

e0=8.85*10^(-12);

e=1;

g=9.8;

rk=sqrt((qp*q)/4*pi*e0*e*m*g);

t1=[0:0.08:0.48];

r1=r0-(g*t1.^2)./2

subplot(2,2,1);plot(t1,r1)

grid on

xlabel('t')

ylabel('r1')

v1=r0./t1-g.*t1

subplot(2,2,2);plot(t1,v1)

grid on

xlabel('t')

ylabel('v1')

t2=[0.5:1:6.5];

r2=rk+((qp*q)/(4*pi*e0*e*m)).*(t2.^2)./2

subplot(2,2,3);plot(t2,r2)

grid on

xlabel('t')

ylabel('r2')

v2=((qp*q)/(4*pi*e0*e*m)).*t2

subplot(2,2,4);plot(t2,v2)

grid on

xlabel('t')

ylabel('v2')