Skip to content

Instantly share code, notes, and snippets.

@wactbprot
Created June 12, 2013 13:07
Show Gist options
  • Save wactbprot/5765071 to your computer and use it in GitHub Desktop.
Save wactbprot/5765071 to your computer and use it in GitHub Desktop.
play around --> persp()
remvalue(all);
v1 : %pi*r1^2*h1;
v2 : %pi*r2^2*h2;
v3 : %pi*r3^2*h3;
V : 2*v1 + v2 + 2*v3;
sr1 : diff(V, r1);
sr2 : diff(V, r2);
sr3 : diff(V, r3);
sh1 : diff(V, h1);
sh2 : diff(V, h2);
sh3 : diff(V, h3);
r1 : 34.0/2;
r2 : r1;
dd : 2.0;
r3 : r1 + dd; /* Dichtungsspalt Ventil-Volumen*/
h1 : 12.0; /* Tiefe zum Ventilsitz */
h2 : 900.0; /* Volumen zu 1l */
h3 : 0.1;
/*Unsicherheiten aus Messung*/
ur1 : 0.008;
ur2 : 0.0001;
ur3 : 0.001;
uh1 : 0.001;
uh2 : 0.0001;
uh3 : 0.02;
u: sqrt(
(sr1*ur1)^2 +
(sr2*ur2)^2 +
(sr3*ur3)^2 +
(sh1*uh1)^2 +
(sh2*uh2)^2 +
(sh3*uh3)^2 );
ev(float(V));
ev(float(u));
ev(float(u/V));
Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment