Задача: Создать компьютерную программу, которая сможет вычислить реальную зависимость напряжения от тока (Вольт-Амперную характеристику) для обычного резистора.
Обычно на уроках в школе предполагается, что резисторы подчиняются закону Ома
и ток зависит линейно от напряжения. Однако резистору во время долгой работы
свойственно нагреваться, а значит его сопротивление изменяется. Для большинства
проводников при не слишком больших температурах их удельное сопротивление зависит от температуры следующим образом:
ρ(T) = ρ0(1 + α(T − T0))
где ρ0 — удельное сопротивление при температуре T0, а α =
1\273 К−1 — температурная
постоянная, одинаковая для очень большого количества проводников. Процесс установления температуры можно в первом приближении считать достаточно быстрым, а
потери тепла — подчиняющихся закону Ньютона-Рихмана
Pпот. = β∆T.
Задание:
Постройте график зависимости сопротивления от протекающего тока

Мой код:
import math
import numpy as np
import matplotlib.pyplot as plt

for i in range (1, 50):
for t in 300, 400, 500, 600:
T0 = 200
T = T0 + t
if t == 300:
c = 443.6
if t == 400:
c = 484.1
if t == 500:
c = 524
if t == 600:
c = 591

pl = 8902
l = 10
s = 0.01
a = 14.01
q = a*(T-T0)
Q = c*s*l*pl*(T - T0) + q
α = 1/273
ρ0 = 8.700 * 10 - 8
p = ρ0*(1 + α*(T - T0))
R = p * l / s
t = Q / i**2 * R
i = math.sqrt(Q/R)
round(R, 1)
round(i, 1)
i = np.arange(0, 10.5, 0.5)
R= np.arange(0, 10.5, 0.5)
plt.scatter(i, R)
plt.show()
Что не так?

#include

#include

#include

#include

#include

#include

typedef struct

{

char imya[40];

char familiya[40];

char otchestvo[40];

int telefon;

int texpasport;

float pages ;

}gor;

void form(char mass[10], char *imya, char *familiya, char *otchestvo, int telefon, int texpasport)

{

file *f;

gor y;

int i;

if ((f=fopen(mass,"w"))==null) exit(1);

for (i=1; i< =3; i++)

{

cout< < "vvedite imya: ";

cin> > imya;

cout< < "vvedite familiya: ";

cin> > familiya;

cout< < "vvedite otchestvo: ";

cin> > otchestvo;

cout< < "vvedite telefon: ";

cin> > telefon;

cout< < "vvedite texpasport: ";

cin> > texpasport;

fwrite(& y,sizeof(gor),1,f);

if(ferror(f)! =null)

exit(2);

}

fclose(f);

}

void vivod(char mass[50])

{

file *f;

gor z;

if ((f=fopen(mass, "r"))==null) exit(3);

cout< < "fail: "<

while (! feof(f)) {

fread (& z,sizeof(gor),1,f);

if (! feof(f))

{

cout< < "imya: "<

cout< < "familiya: "<

cout< < "otchestvo: "<

cout< < "telefon: "<

cout< < "texpasport: "<

}

fclose(f);

}

}

void del(char fo[10], char fi[10])

{

file *f1,*f2;

gor d;

//char temp[40];

if ((f1=fopen(fi,"r"))==null) exit(4);

if ((f2=fopen(fo,"w"))==null) exit(5);

//cout < < "vvedite nazvanie na ydalenie - ";

//cin > > temp;

for (int i=0; i< 3; i++)

{

fread(& d,sizeof(gor),1,f1);

if (feof(f1)) break;

if (i> 2)

fwrite (& d,sizeof(gor),1,f2);

}

fclose(f2);

fclose(f1);

}

void dobav(char fi[10], char fo[10],char *imya, char *familiya, char *otchestvo, int telefon, int texpasport)

{

int q;

int i;

gor d;

file *f1,*f2;

char temp[40];

if ((f1=fopen(fi, "r"))==null) exit(4);

if ((f2=fopen(fo,"r+"))==null) exit(5);

cout < < "vvedite familiu, pered kotorym nuzno dobavit element: ";

cin > > temp;

for (i=0; i< 10; i++)

{

fread(& d,sizeof(gor),1,f2);

if (feof(f2)) break;

if (strcmp(d.familiya, temp)) break;

q=i;

}

for (i=0; i< 10; i++)

{

if (feof(f2)) break;

fwrite (& d,sizeof(gor),1,f2);

if (i==q-1)

{

cout< < "vvedite imya: ";

cin> > imya;

cout< < "vvedite familiya: ";

cin> > familiya;

cout< < "vvedite otchestvo: ";

cin> > otchestvo;

cout< < "vvedite telefon: ";

cin> > telefon;

cout< < "vvedite texpasport: ";

cin> > texpasport;

cout <

}

}

fclose(f2);

fclose(f1);

}

void main()

{

int n=10;

//cout < < "vvedite kol-vo zapisey: - ";

//cin> > n;

//cout <

form("f1.txt",n);

vivod("f1.txt");

del("f1.txt","f2.txt");

cout < < "posle udalenia: \n";

vivod("f2.txt");

getch();

n=1;

system("cls");

dobav ("f2.txt","f2.txt");

cout < < "posle dobavlenia: \n";

vivod ("f2.txt");

getch();

}

uses crt;

var i,sum,k: integer;

begin

      sum: =0;

      k: =0;

      writeln('вводите числа. 0 - конец последовательности');

      repeat

          read(i);

          if i mod 2=1 then begin

              sum: =sum+i;

              inc(k);

          end;

      until i=0;

      writeln;

      writeln('среднее  арифметическое нечетных чисел: ',sum/k: 5: 3);

end.